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Plate Motion Modelling 

Plate motion modelling has evolved from global models to the study 
of motions of small plates, and to the Investigation of possible plate 
deformation. Recent research under this grant has Involved: 

1. Verification of the relative motions of the 
Rivera plate with respect to its neighbors, l.e., 

North America, Cocc'S, and Pacific. 

2. Reassessment of the Interactions of the Caribbean 
and South American plates. 

3. Interaction of Indlan-Australlan plate (the Indian 
plate undergoes deformation) with the Eurasian plate 
In Southeast Asia, In an attempt to place bounds on 
Aslan deformation. 

4. Detailed study of the Pacific-North America plate boundary. 

The main technique used consists In assuming that the global plate 
motion model RM2 furnishes kinematic boundary conditions, and to use 
local geological and geophysical evidence for departure from Ideal rigid 
plate tectonics. The bulk of this research was presented at the IGC 
meetln<’ In Paris In July 1980, by Prof. Minster. 

Rivera plate motions 

The motions of the Rivera plate relative to Its neighbors (Pacific, 
North America, Cocos) have been determined by Inversion of data along 
the East Pacific Rise and the Rivera Fracture Zone (Minster and Jordan, 
1979). All data used are from the oceanic environment. The resulting 

set of poles and 20 error ellipses are shown on Figure 1, and calculated 
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relative notions at selected points on Figure 2. Of Interest is the 
nodel prediction of right lateral strike slip relative motion between 
Rivera and North America in a NH direction, at a rate of 1*3 cm/yr. The 
direction is compatible with offshore scarps near Tree Marlas islands* 
This means that 1) very significant changes in the nature of plate 
interactions, from subductlon (COCO-NOAM) to the South to strike-slip 
further north, and 2) this boundary may be the site of occasional 
strike-slip events* 

Figure 3 shows a focal mechanism for the 1966 December 15 event 
(mjj “ 5*1) in that region* Although the data are relatively few for 
such a small event, the solution is fairly well-constrained and is 
consistent with the predicted plate motions* Further confirmation will 
be obtained by studying the single largest earthquake in that region in 
the past 20 years (1976 February 9, mjj ■ 5.5). 

Figure 4 (from Lowman, 1980) illustrates the three plate boundaries 
discussed next* They are 1) the Caribbean-South America boundary, 
described as right-lateral strike slip on this figure, 2) Che 
Indian-Australlan boundary, described as left-lateral strike-slip on 
Figure 4, and 3) the Pacific-North America boundary, with complications 
in the West North American continent shown on Che map* 

Caribbean-South America Boundary 

The motions of the Caribbean plate were first described by Jordan 
(1975). We have reviewed these results in the light of Che new models. 
Figure 5 is a simplified tectonic map of the northern portion of South 
America, shown in an oblique mercator projection about the RM2 CARB-SOAM 
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pole 4 The boundary appears as a fairly broad zone of continental 
deformation* with conjugate strike-slip faulting* and a significant 
component of shortening between the Curacao Ridge and the frontal 
thrust* However* the amount of shortening predicted by the model Is 
somewhat less than that found by Jordan (1975)* Again* rapid (5.4 
cm/yr) left-lateral strike-slip motion Is predicted between Central 
America and the Panama Basin* pointing to the Panama escarpment as a 
possible seismic gap* 

The evidence for present-day plate deformation Is shown on Figure 
6a and 6b* This Is a contour map of the logarithm of cumulative moment 
release (per degree sequare) In the period 1960-1980* Only shallow (h < 
70 km) events with ro|, > 4 have been used and smoothing was performed 
with a Gaussian filter of aperture ~ 0*3"* Most of the activity Is 
confined between the major thrusts confirming the interpretation of a 
diffuse plate boundary* 

Direct evidence of subductlon of the Caribbean plate under the 
Curacao ridge Is found In the multichannel seismic date of Talwanl et 
al* (1978) (Figure 7)* The dip of seismic horizons thickening of 
sedimentary layers to the South indicate that slow subductlon Is 
probably currently taking place* 

Another confirmation that RM2 provides the correct boundary 
conditions Is shown on Figure 8 (from . ‘<1* 1978)* Near Caracas* 
Venezuela* N-NW treading faults offset the large right-lateral EW fault 
systems (San Sebastian* El Pilar) In a right-lateral direction* Yet the 
complex multiple event of 1967 July 29 (the largest in recent history) 
shows a left-lateral focal mechanism* This Is easily explained If the 


5 


regional naxlmum compreaalva atraaa la controlled by relative plate 
motions* Figure 8 shows that the compression axis for this focal 
mechanism is consistent with the RM2 prediction. 

Thus* RM2 does not change the basic conclusions of Jordan (1975) 
and is consistent with additional evidence published since that study. 

Indian Plate Deformation 

Figure 9 depicts the geometry of the Indlan-Australlan plate motion 
relative to Eurasia in an oblique mercator projection about the 
(INDI-EURA) pole. The Australian plate was allowed to have a motion 
Independent of the Indian plate so as to better satisfy data along the 
SE Indian Ridge (Minster and Jordan, 1978). The hypothetical boundary 
was arbitrarily assumed to coincide with the Ninetyeast ridge, in view 
of its high level of seismicity (Stein and Okal, 1978). Although the 
differential motion of the Indian and Australian portions is small 
(Figure 9) , it is calculated to be of the order of 1 cm/yr SE-NW 
compression. This is consistent with compression axes determined by 
Stein and Okal (1978) shown on Figure 10. 

In an attempt to refine the analysis, we Investigated the 
possbllity of interpreting the Java-Sumatra trench system as a 
convergent plate boundary between India (alternatively Australia) and 
Eurasia, keeping in mind possible significant relative motions between 
Southeast Asia and the Siberian platform. Figure 11 shows seismic slip 
vectors determined by Fitch (1972), plotted on an oblique mercator 
projection about the INDI-EURA pole. The observations are consistent 
with the general interpretation of an INDI-EURA convergent boundary, but 
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the scatter Is too large to detect snail discrepancies* Thus, although 
this permits us In principle to place bounds on the rates of plate 
deformation In the region, these bounds are not tight enough to be 
useful* 

Paclflc-Horth America 

Figures 12a and 12b show a contour map of log cumulative moment 
along the Western North American plate boundary* It Illustrates the 
diffuse and complex nature of the boundary, particularly in the Basin 
and Range province, and In Alaska* On the other hand, the plate margin 
appears to be simpler In the Gulf of California and along the Canadian 
coastline* 

A simplified Interpretive tectonic map of this area Is shown on 
Figure 13, again In an oblique mercator projection about the PCFC-NOAM 
pole* Because of the complications associated with Internal deformation 
of North America In the Basin and Range and Rio Grande Rift, the RM2 
data set along this boundary (shown on Figure 14) did not Include data 
Inbetween 32 'N and 50'N along the west coast of North America* Thus, it 
was hoped, no contamination of the model by complex tectonic 
deformations would occur* 

On the other hand, comparison of Figures 13 and 14 shows the 
following* 

1) Because of the paucity of tectonic information about northern 
Mexico, it Is not easy to assess possible contamination of data In the 
Gulf of California by deformation to the East In the southern extension 
of the Basin and Range and Rio Grande regions* 
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2) B«cau8« of obllqu* ■ubductlon along tha Quean Charlotta fault » 
earthquake slip vectors In that region may not be representative of true 
plate motion directions* 

3) The active Denali fault system Is Alaska may jeopardize the 
correct interpretation of slip vectors In the Eastern Aleutians and the 
Gulf of Alaska* 

A) Known effects of the downgoing slab have been shown to bias 
fault plane solutions In the western Aleutians (Engdahl et al* . 1977)* 
5) The plate boundary between Eurasia and North America Is very 
uncertain* and there Is no guarantee that subductlon In the Kuriles Is 
between Pacific and North America* 

In view of these many uncertainties, one can legitimately question 
Inferences based on the RM2 PCPC-NOAM rotation vector* To remove biases 
specifically associated with this plate boundary, we have performed a 
global plate motion Inversion after removing the entire data subset just 
described. Including the single rate datum between PCPC and NOAM at the 
mouth of the Gulf of California* 

The resulting pole Is shown on Figure 15 with the RM2 pole and 
their respective 2o error ellipses* It Is truly remarkable that the 
solutions are not resolvable from each other* The PCFC-NOAM motion 
determined from the rest of the world Is practically Identical to that 
determined from the plate boundary* 

Figure 16 shows that the predicted direction of relative motion In 
central California is distinguishable from the azimuth of the San 
Andreas fault system at the 2o level, whether or not PCFC-NOAM data are 
Included* This In turn means that the San Andreas direction is not 
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r 0 pr*s«ntAtlv« of relative oKitlon between PCFC and NOAM, but only of 
relative motion between the Pacific plate and the Great Valley region of 
California, or even between the Sallnaa block and the Great Valley* 

Geologically and geodetlcally determined ratee of allp along the 
San Andreas In central California yield a value of 3.7 cm/yr» If one 
accounts for auxiliary faults, this fate reaches up to 5.0 cm/yr* It la 
likely that the RM2 rate of 5*7 cm/yr Is an upper bound. Figure 17 
shows a plana geometry relative motion diagram between PCFC, NOAM and 
California (CALI) and Mojave (MOJA) blocks. 

CALI-NO^ is obtained by differencing the RM2 PCFC>N0AM speed and 
the slip vector along the San Andreas fault* Since there is good 
evience that this vector should have a non-aero component in a westerly 
direction (e.g., Zoback and Zoback, 1979) as shown on Figure 18, 
Interpretation of Figure 17 leads to an upper bound of 5*0 cm/yr for the 
San Andreas fault system, in agreement with geological data. If we 
assume that the Garlock fault is the site of 0.7 cm/yr of left-lateral 
displacement, then a lower bound of 1.3 cm/yr for the San Andreas fault 
system is obtained to Insure compression in the Transverse ranges* This 
lower bound is Increased to 2.0 cm/yr if 0.7 cm/yr thrust displacement 
is allowed across the Garlock-VRilte Wolf fault system. On the other 
hand, unless we have simultaneously a fairly high rate of slip on the 
San Andreas fault system, and a significant component of compression 
across the Garlock fault, the model predicts some amount of EW 
compression between the Mojave block and the stable part of NOAM, for 
which evidence Is lacking. 

The alternative is, of course, to assume that EW shortening takes 
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place b«ew««n the San Andraae fault and th« Pacific plate* Then the 
PCFC-c;/(M vector would trend In a more northerly direction than both the 
Sen Ancdreae fault and the PCPC-NOAM RM2 vector end EW Beeln end Range 
extenelon It easily accommodated. 

Thus, analyele of Inetentaneoue rii^Aatlve motlone on a reglnal 
baele. using the global model RM2 to provide klnematlcal boundary 
conditions leads us to the prediction of an EW component of shortening 
across the Salinas block* 

Independent confirmation Is provided by the orientations of 
compression axes of crustal events (Figure 19) and of fold axes (Figur? 
20) compiled respectively by Gawthrop (1978) and Burford (1967), which 
are consistent with a significant component of compression In a 
direction perpendicular to the San Andreas trend* 

Quantitative estimates require the determination of bounds on 
azimuth and direction of Basin and Range opening* This work la 
underway* 

Seismicity and Subductlon Processes 

We have studied correlations between coupling on the 
subductlon-zone thrust plane, parameterized by regional maximum 
earthquake magnitudes, with other parameters of the sbductlon zone* We 
found that convergence rate and lithospheric age correlate with 
coupling* It was also determined that penetration depth correlates with 
lithospheric age and contact width with convergence rate, results 
compatible with work done by Xsacks et al* (1968), Vlaar and Wortel 
(1976), and Wortel and Vlaar (1978). These correlations suggest that 
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large earthquake generation correlatea with Interface geometry or with 
age and rate in aoma other manner. Thta study was completed and the 
paper (Ruff, L., and H. Kamimori, Seismicity and the rubduction 
process) is now in press In the Rnyeici of, _the Earth and ?lanetary 
Interiors . 

Study of Slow Earthquakes 

Detailed analysis of IDA data was made to study a slow earthquake 
which occurred on December 16, 1976, in the Bouvet Islands region. 
Although the body>wave magnitude of this earthquake is very small (m(, •* 
5.4) it excited long-period Rayleigh waves with a very large amplitude. 
The seismic moment is determined to be 1.4 x 10^® dyne-cm, which is much 
larger than that for earthquakes with conq>arable body-wave magnitude. 
The mechanism of this event is consistent ..ith the geometry of the 
transform fault in the Bouvet Islands region, and suggcjta that the 
plate motion there is somewhat slower than that on other older transform 
faults. 

A similar analysis is currently being made for subduction-zone 
events. 

Study of Free Oscillations 

The earthquake that occurred in the Ecuador and Colombia border on 
December 12, 1979 is large enough to excite various long-period modes. 
Preliminary analysis identifies oS2, oS3, oS4 and oSo modes very 
clearly. We have determined the mechanism of this event, and are 
preparing the data sets to be used for more extensive analyses of the 
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free oaclllatlons» This type of study provides Information Intermediate 
In period to body waves and longer duration plate motions* 
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Figure 10 


Fig# 12 » Focal mechauismt magnitudes, and moments for major earthquakes near the 
Minetyeast Ridge* The mechanism of the 1928 event could not be realiably determined* 
Note that the 1939 and 1964 events on the ridge shov strike slip motion, the 1955 
event is a normal fault, and the 1974 e^vent is a thrust fault* The moments of the 
1928 and 1939 earthquakes are much greater than any of the more recent events* 
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Urn r«vl«w eh« progress in modeling Instsntaneous plate kinematics, 
with an emphasis on recently developed models of present-day plate 
motions derived by the syi«tematlc inversion of globally distributed 
data sets. The U-plate model RM2 (Minster and Jordan, 1978) is 
described, and the relative and absolute velocity flelda associated 
with this model are summarized in a series of maps. 
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Mod«Ung Ktiwaatlct 

In eh* •*v*r*l y**r* following Wllion** (1965a) outlln* of eh* basic 
klnanatlcal concapts, plaea eactonlcs was tranafomad from a bold now 
hypoehasls to an astabllshad Chtory. Many naw Idaas pronulgatad during 
this critical parlod wara synthaslzad In a book by Coulomb (1969» 1972); 
tan years latar these idaas continue to guide research In solid earth 
sclanc*. 

On* of the most spectacular succassas of the plat* theory was the 
early solution of the Instantaneous klnematlcal problem* namely: In the 

present geological epoch* where are the plate boundaries and what are the 
plate motions? Wilson (1965a) recognized that the three klnematlcal boun- 
dary types - divergent* convergent and transform - correspond to three 
classes of deformation zones - spreading oceanic ridges* active mountains 
and Island arcs* and major shear faults - "connected Into a continuous 
network of mobile belts about the earth which divide the surface Into 
several large rigid plates." McKenzie and Parker (1967) and Morgan (1968) 
Introduced the geocentric angular velocity vector as a description of 
Instantaneous plat* motion and began to determine these vectors by fitting 
sea-floor spreading rates* strikes of transform faults and azimuths of 
earthquake slip vectors. This pioneering work culminated In la Plchon's 
(1968) masterful synthesis of plate movements past and present* In which 
he produced a model for the Instantaneous relative velocities of six 
major plates: Africa* America* Antarctica* Eurasia* India and Pacific. 

Since 1968 the klnematlcal models have been considerably refined by 
numerous special studies of specific geographical regions. Among these 
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arc McKanzia and Sclatar'a (1971) traatlaa on plato tactonlcs In tha 
Indian Ocaan» Pieman and Talwanl'z (1972) analyslz of aaa-floor apraadlng 
in tha North Atlantic, tha atudy by Holnar (1975) of tha South 

Pacific, Hay* a (1974) aodallng of plata Intaractlona In tha East 
Pacific, and Jordan's (1975) study of tha Carlbbaan. 

Tha nunarlcal algorithms for confuting tha angular valocltias hava 
also baan Improvad. In tha original papers by McKanzla and Parkar (1967) 
and Morgan (1968) pola positions vara constructed by graphical methods. 
Thasa methods wars automated by La Plchon (1968) to handle large data 
sets; ha minimized tha root**maan-squara difference batvean tha Observed 
and computed azimuths of relative motion by searching over a grid of 
possible pola positions. McKenzie and Sclatar (1971) discussed various 
fitting criteria and devised fitting algorithms based on matrix solutions 
to tha minimization problem. As they noted, the problem Is complicated 
by the fact that, although the relative velocity vector at any point on 
a plate boundary Is a linear functional of the angular velocity vector, 
the unit vector describing the azimuth of relative motion Is not. In 
practice this nonlinearity poses no difficulties for most geometries and 
can be accommodated by the Iterative solution of suitably linearized 
equations. 

Thu Instantaneous kinematic problem for the global plate mosaic 
was attacked using iterative, nonlinear least-squares Inversion methods 
by Chase (1972) and Minster et (1974). Since all angular velocity 
vectors describing the relative motions among N plates are tmiquely 
specified by 3(M-1) parameters, to determine a global model from, say,. 


m 


A- 1- 5 


N data (racea and/or azimuths at arbitrary boundary points) requires the 
solution of an N K 3(M-1) system of equations, where N is typically an 
order of magnitude or more greater than M. In this formulation the con»~ 
patlbility conditions (e*g., triple- Junction closure) are automatically 
satisfied. Because the problem is over-determined an exact solution to 
this system generally does not exist, and the particular approximate 
solution chosen depends on the fitting criterion satisfied. Minster 
at al . (1974) advocated the minimization of a fitting function which 
is the sum of the squares of the differences between the observed and 
computed data, each normalized by an estimate of the standard observational 
error. If the data are assumed to have errors which are normally distri- 
buted, this choice yields the classical maximum likelihood estimate of the 
model perturbation, and the linear theory can be used to estimate uncer- 
tainties in the model Induced by uncertainties in the data. The theory 
also allows the computation of what Minster et, al. (1974) have termed 
"data Importances", quantities which are indicative of the information 
distribution among th^ data and which are useful in assessing a model's 
fit to the data set. 

Our experience in applying this algorithm to the plate kinematic 
problem has demonstrated to us its effectiveness. Although the algorithm 
Involves the linearization of a nonlinear problem, convergence to a good 

p i 

solution is generally rapid, and no difficulties with spuriOvts local | 

I 

minima have been encountered. The uncertainties in the model parameters 
derived from the linear theory have proven to be effective measures of 
the model errors caused by random observational errors, although, of 
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course* they do not account for data bias. 

TWO recent* Independent studies by Chase (1978) and Minster and 
Jordan (1978) have presented global kinematical models derived by the 
inversion of improved data sets. The wealth of new data in the ocean 
basins is impressive: every major oceanic ridge is now well 

populated with both azimuths and rates of relative motion (Figure 1)* 
and the uncertainties in the angular velocity vectors have been much re- 
duced. Although .some significant differences between Chase's solution 
and ours are revealed by a detailed comparison* most are ascribable to 
slight differences in data selection and interpretation. Clearly the 
models can be refined* but we anticipate that any future modifications 


will be minor. 


occurs, then the azimuth residuals (observed-computed) for transforms 
characterized by left-lateral faulting should be systematically more 
positive than those for right-lateral transforms. In fact, sintstral 
transform faults show slightly more negative residuals than dextral 
transform faults, opposite to that expected from the leaky transform 
hypothesis. Ihe difference Is not statistically significant, however, 
and the RM2 fit to the transform fault data appears to be unbiased. 
Iherefore, there Is no evidence In our results to Indicate that leaky 
transforming commonly occurs. 

A complete and convenient description of the Instantaneous relative 
motions among M plates is given by the geohedron, the set of 3(M-2) 
angular velocity vectors for plate pairs sharing common boundaries 
(McKenzie and Parker, 1974). Ttie 27 relative velocity vectors of the 
BM2 geohedron are listed In Table 1 In both spherical and cartesian co- 
ordinate systems. The poles and their associated 95Z confidence regions 
are depicted In Figures 2-4. Also shown are the best-fitting poles (BFP's) 
for Individual plate boundaries determined by separate inversion of data 
sets for specific plate pairs, wherever possible (Minster and Jordan, 1978, 
Table 3). The differences between BM2 and the BFP set result from the 
2M*5 vector compatibility conditions implicit In the simultaneous Inversion 
scheme. In most cases these differences are small, and BM2 provides a 
nearly optimal fit to localized data sid>sets. Significant differences 
do exist, however, for AFRC-NOAM, INDI-PCFC, INDl-ANTA and INDI-AFRC. 

The discrepancy for AFRC-MOAM involves RMZ's misfit to transform azimuths 
In the FAMOUS area, discussed above. The discrepancies between RM2 and 
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detailed history of plate motions as well as the spatial scale lengths 
over which the primary bathymerlc data are averaged. Complications arise 
In regions such as the central North Atlantic where an apparent reorien- 
tation of the motions has occurred in the last several million years 
(Macdonald, 1977; Fox et al. , 1979). Here two trends are superposed, 
an older trend of about S80 *’e represented by the mean azimuths of major 
transform faults (e.g.. Oceanographer and Atlantis) and a younger, nearly 
east-west trend of smaller, more recent faults (e.g.. Transforms A and B 
In the FAMOUS area and en echelon faults within the Oceanographer F.Z.). 

In this case we used the younger trends In the RM2 data set, although* 
surprisingly, they were not fit by the model; RM2 instead predicts azimuths 
compatible with the mean trends of the major transforms. (This misfit Is 
Indicative of Internal Inconsistencies In the data pertaining to the 
relative motions about the Azores triple junction not readily associated 
with the time-averaging problem; see Minster and Jordan, 1978, pp. 5341- 
5343, for further discussion.) Inmost cases, however, complications of 
this sort are not apparent, and the large-scale trends of transforms 
are probably representative of plate motions during the last several 
million years. This hypothesis Is supported by the good agreement, between 
transform fault azimuths and the slip vectors of transform fault earth- 
quakes, which truly do measure "Instantaneous" motions. 

The KM2 data set has also been examined for systematic discrepancies 
between observed transform fault azimuths and model-predicted 41rectlox». 
One possible source for such discrepancies could be widespread "leaky" 
transform faulting In the sense of Menard and. Atwater (1969). If this 
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Ralatlva Motions ; Model RM2 

Ihe model designated Relative Motion 2 <RM2) resulted from the simul- 
taneous Inversion of 110 spreading rates, 78 transform fault azlsuiths and 
142 earthquake slip vectors (Minster and Jordan, 1978) ; in isany cases the 
angular velocity vectors have been significantly Improved with respect to 
our previous model RMl (Minster et al. , 1974) • BM2 represents the instan- 
taneous relative motions among 11 major plates: Africa (AFRG), Antarctica 

(Am), Arabia (ARAB), Caribbean (GARB), Cocos (COCO), Eurasia (EURA), 

India (INDI) , Nazca (NAZC), North America (NOAM), Pacific (PCFC) and South 
America (SOAM). 

In compiling the RM2 data set, shown in Figure 1, the data from Intra- 
contlnental environments such as the Alplde belt and complex tectonic 
regions such as the western Pacific were generally excluded, because these 
areas exhibit conplexltleo not easily described by rigid plate kinematics, 
and the assumptions fundamental to a simple plate model may not apply. 

Hence the details of Mediterranean, Asian and Indonesian tectonics are not 
represented by this Instantaneous klnematical model. 

Of course, no klnematical model based on geological data is really 
"Instantaneous"; in particular, the spreading rate estimates and transform 
fault azimuths average in some way the motions over substantial time 
intervals. For a spreading rate datum the pertinent time Interval is 
simply the age of the magnetic anomaly used to obtain the rate. We generally 
employed anomalies 2 and 2', so that the mean averaging Interval for the 
RM2 rate data is less than 3 tty. The relevant averaging intervals for the 
transform fault azimuths are more uncertain, however, and depend upon the 
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thft BFP’s for the INDI poles are thou^t to be diagnostic of Internal 
deformation within the Indian plate (Minster and Jordan, 1978). although 
we concede that the prrblem could result from other, unrecognized sources 


of data bias. 
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Absolute Motions 

Nothing within the axiomatic framework of plate tectonics points to 
a natural* absolute reference frame In which to describe plate motions. 

Hie need to choose such a frame aidses only when one enlarges the class of 
observations to be explained beyond the manifestations of plate-plate 
Interactions or when one attempts to model the dynamics of plate motions. 

A particular frame of Interest in discussions of plate dynamics Is one 

I 

fixed with respect to the average position of the deep mantle* assused to 
be rigid or at least to have typical Internal motions much slower than 
the motions of the plates; we refer to this frame as the mean mesospheric 
frame. 

Ullson (1963, 1965b) proposed that various Island chains and aselsmlc 
ridges are generated by the motions of the crust over hotspots fixed in the 
mantle. Morgan (1971»1972a*b)expanded this concept and presented a model 
of 15 plates moving in a slowly changing frame defined by a score of hot- 
spots* each of which he associated with an upwelllng convective plxime. 

In a test of the Wllson-Morgan hypothesis Minster et al. (1974) used RMl 
to construct a model designated Absolute Motion 1 (AMI) by fitting the 
observed azimuths of twenty hotspot traces. They found no cause to reject 
the Wllson-Morgan hypothesis and concluded that these hotspots showed no 
significant relative motions In the last 10 My» although their ability to 
resolve such motions was admittedly limited. Other studies have concluded 
that^ averaged over lunger time Intervals (>40 Ify)* some hotspots have 
relative speeds of the order of 1 cm/yr (Morgan, 1972b; Burke ejt ^.* 1973; 
Molnar and Atwater* 1973; Molnar and Francheteau* 1975). Nevertheless* 
these relative speeds appear to be much smaller than the typical absolute 
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speeds of the faster plates, upon which the hotspot traces are best de- 
veloped, so that the concept of a hotspot reference frame remains useful 
for the short time scales used to derive instantaneous motions* Further- 
more, if the directions of hotspot migration are nearly random, this 
reference frame may be a very good approxiisatlon to the mean mesospheric 
frame. Chase (1978) and Minster and Jordan (1978) have derived absolute 
motion models based on hotspot data apd compatible with their revised 
estimates of relative motions. 

The models based on the Ifilson-Morgan hypothesis have been compared with 
other proposed criteria for fixing plate velocities in the mean ntesoapheric 
frame. One is to require that the lithosphere as a whole possess no net 
rotation with respect to this frame, a criterion discussed and applied by 
Minster e£ al. (1974), Lllboutxy (1974) and Solomon and Sleep (1974). The 
no-net-rotatlon model is one member of a more general class of mechanical 
models obtained by balancing a system of forces acting on the plates 
(Solomon and Sleep, 1974; Solomon ad., 1975; Gordon et al., 1978). 

Another class of absolute motion models is obtained by fixing a particular 
plate with respect to the mean mesospheric frame, such as the Antar:tlc 
(Knopoff and Leeds, 1972), African (Burke and Wilson, 1972) or Caribbean 
(Jordan, 1975) plates. TUllis (1972) suggested that downgoing slabs act 
as anchors in the mantle* and thus specify an absolute reference frame. 

This concept was discussed by Kaula (1975), who proposed a general boundary 
velocity minimization criterion, and Jordan (1975) , who used the' idea to 
rationalize his conclusion that the Caribbean plate is essentially fixed 
with respect to the mean mesospheric frame. 
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In Table 2 and Figuraa 5 and 6 wa present the absolute notion nodal 
designated AMl-2. This nodal was derived iron RM2 by fitting the aslnuths 
of nine hotspot traces and the rates iron five (Minster and Jordan, 1978) • 

The data set was restricted to Include only those constraints on hotspot 
nigratlon pertinent to the last 10 My. This tine span is really the 
ninlnun Interval for which good hotspot data can be obtained, although 
it exceeds by over a factor of three the nean averaging Interval for the 
relative notion data. Consequently, the data set is donlnated by infoma- 
tion iron the rapidly novlng plates in the Pacific Ocean; no Atlantic or 
Indian Ocean hotspot traces were employed. 

The AMl-2 velocity field depicted in Figures 5 and 6 shows a nuofter 
of interesting features which are common to most proposed absolute velocity 
models. The contrast between fast moving oceanic plates and slower moving 
continent bearing plates, noted by Minster et al . (1974), is striking, 
although the rapid motion of INDI, which has the same proportion of conti- 
nental area as ANTA, prevents an exact dlchoton^. The four fastest plates 
(PCFC, COCO, NAZe, INDI) are also the only plates being sid>ducted along a 
significant fraction of their boundaries, a correlation noted by Jordan and 
Minster (1974) and discussed in detail by Forsyth and Uyeda (1975). As 
Aggarwal (1978) has shown, an even better correlation between plate geo- 
metry and absolute velocities is obtained if the percentage length of ridge 
is Included in the analysis. There is a general increase of absolute velo- 
cities away from the spin axis, which is reflected in the distribution of 
global seismicity (Solomon et aJ.. , 1975). Systematic attempts to relate 
these correlations to plate mechanics have been made by Solomon et (1975), 
Forsyth and Uyeda (1975) and Chappie and Tullls (1977), but similar patterns 
do not seem to apply to Tertiary tl.nes (Solomon ^ al - 1977; Jurdy, 1978), so 
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caution roust be exercised in attributing to them any dynamical significance. 

Mien examining the absolute velocity field In Figures 5 end 6» It 
should be kept In mind that AMl-2 Is associated with much larger uncer> 
talntles than SM2, because the resolution of the hotspot data Is Intrinsi- 
cally limited and because caveats must be attached to the Wllson-Morgan 
hypothesis. The relative errors In the directions of absolute motion are 
generally small for points on the fast moving plates » but they can exceed 
lOOZ for points on the slowest plates, particularly EUAA and ANTA. In 
these cases, the directions are not usefully constrained by the data used 
to derive AMl-2. A more complete discussion of these uncertainties Is 
given by Minster and Jordan (1978). 

The large uncertainties associated with AMl-2 should also be kept 
in mind when comparing this model with other estimates of absolute motions. 
Models based on some of the alternate hypotheses discussed above can be 
constructed from Table 2 by vector addition. A detailed comparison 
(Minster and Jordan, 1978, pp. 5348-5349) yields the following conclu- 
sions: the no-net-rotatlon model and the model obtained by fixing AFRC 

do not provide adequate fits to the hotspot data, whereas any of the three 
plates GARB, EURA or ANTA could be stationary with respect to the hotspot 


frame. 
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Dlscusalon 

Plate tectonics is the unifying loodel which ties surface geology 
to planetary evolution. As noted by Deffeyes (1970) a useful and slnqple 
measure of global tectonic activity Is the rate of production and de- 
struction of lithosphere associated with plate motions. l!hls calculation 
requires only mlnlual Information about plate geometry (Deffeyes, 1970; 
Chase, 1972). Table 3 gives the area- budget for the RM2 plate set; In 
the computation we have assumed symmetric spreading at ridges, unilateral 
lithosphere destruction at trenches, and that the area loss Is shared 
equally by both plates In a continental collision. Clearly, the net growth 
of some plates (e.g. , ANTA, AFRC, SOAM AND NOAM) and the net shrinking of 
others (e.g., PCFC, INDI, NAZC, COCO) observed at the present time must 
be satisfied by any successful dynamical model, a point emphasized In the 
work of Hager and O'Connell (1978). 

There are conspicuous Instances where the ideal rigid plate model 
fails to describe adequately the complexities of tectonic Interactions. 
These may include continental collisions (e.g., Mblnar and Tapponnler , 
1975), zones of Incipient rifting (e.g., Courtlllot, 1979), and regions 
of deformation along continental margins (e.g., Jordan, 1975). In these 
complex areas the rigid plate model has only limited utility, but It 
does provide the kinematlcal boundary conditions which must be satisfied 
by local models of deformation. 
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Table 1; 


Table 2: 


TABLE CAPTIONS 


First plate named moves counterclockwise with respect to 
the second. Cartesian coordinates defined with respect 
to spin axis and Greenwich meridian. 

Plates move counterclockwise about poles. Cartesian 
coordinates' defined with respect to spin axis and 
Greenwich meridian. 
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TABLE 1; RM2 CEOHEDRON 
RELATIVE ROTATION VECTORS 


Plate Pair 

0 

•N 

♦ 

•e 

(It 

deg/My 

deg/Ify 

deg!^ 

deg^My 

AFRC-ANTA 

9.A5 

-41.70 

0.149 

0.1095 

-0.0976 

0.0244 

AFRC-EURA 

25.23 

-21.18 

0.104 

0.0874 

-0.0339 

0.0442 

AFRC-NOAM 

80.43 

56.36 

0.258 

0.0238 

0.0357 

0.2546 

AFRC-SOAM 

66.56 

-37.29 

0.356 

0.1126 

-0.0858 

0.3266 

ANTA-PCFC 

64.67 

-80.23 

0.964 

0.0700 

-0.4065 

0.8712 

ANTA>SOAM 

87.69 

75.21 

0.302 

0.0031 

0.0118 

0.3021 

ARAB-AFRC 

30.82 

6.43 

0.260 

0.2222 

0.0250 

0.1334 

ARAB-BURA 

29.82 

- 1.64 

0.357 

0.3097 

-0.0088 

0.1776 

CARB-SOAM 

73.51 

60.84 

0.202 

0.0280 

0.0502 

0.1940 

COCO-CARB 

23.60 

-115.55 

1.543 

-0.6098 

-1.2757 

0.6179 

COCO-NAZC 

5.63 

-124.40 

0.972 

-0.5463 

-0.7978 

0.0953 

COCO-NOAM 

29.80 

-121.28 

1.489 

-0.6707 

-1.1041 

0.7399 

COCO-PCFC 

38.72 

-107.39 

2.208 

-0.5150 

-1.6438 

1.3810 

EURA-NOAM 

65.85 

132.44 

0.231 

-0.0637 

0.0696 

0.2104 

EURA-PCFC 

60.64 

-78.92 

0.977 

0.0921 

-0.4701 

0.8515 

INDI-AFRC 

17.27 

46.02 

0.644 

0.4271 

0.4426 

0.1912 

INDI-ANTA 

18.67 

32.74 

0.673 

0.5366 

0.3450 

0.2136 

INDI-ARAB 

7.08 

63.86 

0.469 

0.2049 

0.4175 

0.0577 

INDI-EURA 

19.71 

38.46 

0.698 

0.5146 

0,4087 

0.2353 

INDI-PCFC 

60.71 

- 5.79 

1.246 

0.6066 

-0.0615 

1.0868 

NAZC-ANTA 

43.21 

-95.02 

0.605 

-0.0386 

-0.4396 

0.4145 

NAZC-CARB 

47.30 

-97.57 

0.711 

-0.0635 

-0.4779 

0.5226 

NAZC-PCFC 

56.64 

-87.88 

1.539 

0.0314 

-0.8460 

1.2857 

NAZC-SOAM 

59.08 

-94.75 

0.835 

-0.0355 

-0.4278 

0.7166 

NOAM-CARB 

-33.83 

-70.48 

0.219 

0.0609 

-0.1717 

-0.1220 

NOAM-PCFC 

48.77 

-73.91 

0.852 

0.1557 

-0.5398 

0.6411 

NOAM-SOAM 

25.57 

-53.82 

0.167 

0.0888 

-0.1215 

0.0720 


TABLE 2; MODEL AMl-2 


Plate 

AFRC 

ANTA 

ARAB 

CARB 

COCO 

EURA 

INDI 

NAZC 

NOAM 

PCFC 

SOAM 

Net 

Rotation 
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ABSOLUTE ROTATION VECTORS 


0 ^ W Wy W_ 

*N *E deg/Ny deg^ degAly deg/My 


18.76 

-21.76 

0.139 

21.85 

75.55 

0.054 

27.29 

- 3.94 

0.388 

-42,80 

66.75 

0.129 

21.89 

-115.71 

1.422 

0,70 

“23.19 

0.038 

19,23 

35,64 

0.716 

47,99 

-93.81 

0.585 

-58. 31 

-40.67 

0.247 

-61.66 

97.19 

0.967 

-82.28 

75.65 

0.285 

-51.09 

65.44 

0.251 


0.1221 

-0.0487 

0.0446 

0.0126 

0.0488 

0.0202 

0.3443 

-0.0237 

0.1781 

0.0375 

0.0872 

-0.0879 

-0.5724 

-1,1885 

0.5300 

0.0347 

-0.0148 

0.0005 

0.5492 

0.39 38 

0.2358 

-0.0261 

-0.3907 

0.4347 

0.0983 

-0.0844 

-0.2099 

-0.0574 

0.4553 

-0.8510 

0.0095 

0.0370 

-0.2819 

0.0655 

0.1434 

-0.1953 
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TABLE 3 


AREA BUDGET FOR RM2 


Plate 

Name 

Area Created 
(kmVyr) 

Area Destroyed 
(km2/yr) 

Net Growth 
(km2/yr) 

Relative Area 
Change(l0~^ yr"*) 

AFRC 

0.322 

0.025 

0.297 

3.76 

ANTA 

0.575 

0.026 

0.549 ' 

9.31 

ARAB 

0.024 

0.048 

-0.024 

-4.92 

GARB 

0.001 

0.009 

-0.008 

-2.08 

COCO 

0.174 

0.211 

-0.037 

-12.86 

EURA 

0.060 

0.117 

-0.057 

-0.83 

INDI 

0.339 

0.693 

-0.354 

-5.90 

NAZC 

0.439 

0.505 

-0.066 

-4.43 

NOAM 

0.107 

0.017 

0.090 

1.51 

PCFC 

0.690 

1.207 

-0.517 

-4.79 

SOAM 

0.157 

0.030 

0.127 

3.09 

TOTALS 

2.888 

2.888 

0. 
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FIGURE CAPTIONS 


Figure 1 ; Relative velocities from RM2 at selected points along 
plate boundaries. Mercator projection; latitude 
and longitude marks every 30*. Numbers are rates in 
cm/yr. The RM2 data Is used to outline major plate 
boundaries. Circles are rate data, lozenges and triangles 
are transform fault and slip vector azimuth data respectively. 
Velocity vectors are drawn symmetrically about point of 
calculation In direction of relative motion with length 
proportional to rate. 


Figure 2 : Poles for model ^2, with their 952 confidence ellipses. 

RMl and best fitting poles, where available, are shown 
for comparison. 

Figure 3 ; Poles for model RM2, with their 95% confidence ellipses. 

RR1. and best fitting poles, where available, are shown 
for comparison. 


Figure 4 ; Poles for model RM2, with their 95% confidence ellipses, 
RMl and best fitting poles, where available, are shown 
for comparison. 
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Figure 5 : Absolute motions at selected points from model 

Conventions are similar to figure !• Stars are absolute 
poles, about which plate rotation Is counterclockwise. 
Asterisks are antipoles. 

Figure 6 ; Absolute velocity field from AMl-2. Vectors drawn along 
small circles about absolute poles with exaggeration 
factor 2.22 x 10^, so that 1 cm/yr Is mapped Into 
2 geocentric degrees at 90* from pole. Orthographic 
projections, centered at a) 0"N, 0"E; b) 0“N, 90"E; 
c) 0*N, 180*E; d) 0*N, 90*W; e) North pole; 
f) South pole. 
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Abstract 

Dislocations contribute to both seismic wave attenuation and steady- 
state creep in the mantle. The two phenomena Involve quite different strains 
and characteristic times but they can both be understood with simple 
dislocation models. The roost satisfactory model for creep Involves the 
glide of dislocations across subgrains but rate limited by the climb of 
jogged dislocations In the walls of a polygonlzed network. The Jog 
formation energy contributes to the apparent activation energy for creep, 

making it substantially larger than the activation energy for self-diffusion. 

2 3 

The theory leads to either a a or a law for creep rate, depending on the 
length of the dislocations re.latlve to the equilibrium spacing of thermal 
jogs. 

Attenuation In the mantle at seismic frequencies is probably caused by 
the glide of dislocations la the subgrains. Kink and impurity drag can both 
contribute to the glide time constant. The kink-forma tlon, or Pelerls 
barrier, model for ii-slocatlon glide appears to be a low- temperature, high- 
frequency mechanism most appropriate for pure systems. A small amount of 
impurity drag brings the dislocation glide characteristic time into the 
seismic band at upper mantle temperatures. 
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The attenuation and creep behavior of the mantle are related through 
the dislocation structure. Discussion of the various possible mechanisms is 
facilitated by casting them and the geophysical data in terms of a pre- 
exponential characteristic time and an activation energy. The relaxation 
strength is an additional parameter that can be used to Identify the 
attenuation mechanism. Mobile dislocations In sid>grains» rather than 
cell-walls, have the appropriate characteristics to explain the damping 
of seismic waves In the upper mantle. The grain boundary peak may be 
responsible for attenuation In the lithosphere. 
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I. Introduction 

The high traptratur* cr««p of aoot crystalllno solid* Is controlled 

by dislocation cllnb or glide plus cliad> [Weercnsnn» 1968] . For a wide variety 

of Mtalllc, Ionic and polar crystals the activation energies for creep end 

self-diffusion are essentially Identical, indicating that the notions of 

dislocations are rate-limited by self-diffusion. Furthermore, In the stress 
-2 -5 

range of about 10 to 10 G, where G is the shear modulus, the stress 
exponent of the creep rate is about 3, consistent with dislocation climb or 
the glide of Jogged dislocations. Both of these mechanisms of steady-state 
creep require self-diffusion. The tUlrd-power creep law Is satisfied by olivine 
In the stress range of several hundred bars to several kilebars [Kohls tedt 
and Goetze, 1974 , Goetze and Kohlstedt, 1973, Kohlstedt et al . . 1976]. 

The laboratory creep law for olivine is comnonly assumed to hold for 
the much lower stresses and strain rates appropriate for the mantle and this 
has been the basis for the calculation of many viscosity profiles for the 
mantle. It has also been assumed that the activation energy of creep, 

125 kcal/mole, represents the dlffusional activation energy for the slowest 
moving species. It has recently been found, however, that the activation 
energies for diffusion of oxygen and silicon in olivine are only about 
90 kcal/mole [Reddy e^a^. , 1980, Jaoul e^al. , 1979]. In addition, at 
low-stresses the stress exponent may be closer to two than to three. 

The mechcnlsm of attenuation In crystalline solids also appears to be 
an activated process. Dislocation damping Is a hlgh-temperature mechanism 
that is rate-llmiwed by the diffusion of point defects. The glide of dis- 
locations does not necessarily require self-diffusion and can therefore 
operate on a miJich shorter time scale than cllnd). The purpose of this paper 
Is to Investigate the role of dislocations in both creep and attenuation. 
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II. a. Activated procaaaea 

Theraally activated relaxation proceaaes» Including thoae that are 
controlled by dl£fualon» are functions of a characteristic relaxation time that 
can be written 


T 


T 


O 



( 1 ) 


where T^, the pre-exponential » is the relaxation time at infinite temperature 
and E* Is the effective activation energy. For simple point defect mechanisms 
Is close to the Debye frequency, 'vio sec . For mechanisms Involving 
dislocations or dislocation-point defect interactions depends on other 
parameters, such as dislocation lengths. Burgers vectors, kink and log 
separations, Pelerls stress. Interstitial concentrations, etc. 

The activation energy, depending on the mechanism. Is a composite 
of activation energies relating to self-diffusion, kink or jog formation, 
Peierls energy, bonding energy of point defects to the dislocation, core- 
dlffiision and so on. E* for creep of metals Is often just the self-diffusion 
activation energy, 

Non-elastic processes In geophysics are commonly expressed In terms 
of viscosity, for long-term phenomena* and Q for oscillatory and short-term 
phenomena. It Is convenient to express these In terms of the characteristic 
time. For creep 


T - a/Ge - n/G (2) 

« 

where a, G, e and n are the devlatorlc stress, rigidity, strain rate and 
viscosity, respectively. Thus, the characteristic time can be simply computed 
from creep theories and experiments. 
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Rtlaxation theories, of attonuatlon can !>• described by 

q”^ ■ «t/(1 + W*T^) (3) 

where u Is Che frequency, Q**^ le Che scanderd iMesure of eccenuecion and 

Is Che peak value of eCCenueCion ec otr • 1. In Che ■ancle, where Q la low, we . 
Infer chac Che cheracCerlsClc relesteclon cine Is close Co l/w, Che reciprocal of 

the MsaureinenC frequency. values of Q IndlcaCe ChaC w » 1 /t or u « 1 /t. 

In Chase cases Q 'n* u or aT^, When Q la elowly varying we ere probably in Che aidsc 

of an ebsorpeion band whidi is convenienCly explained as Che reaulc of a dlscrlbuclot 

* 

of relaxaCion Claes [Liu ec al.. 1977; Anderson ec el. , 1977; MinaCer, 1980]. 

Since Q is small in Che higher CenperaCure regions of Che mancle ic 
is likely chac aC aeiamlc periods Che mancle la on Che low-Ceiq>eraCure side of an 
ebsorpeion band. Thus 

q“^ - 2Q^^/wt - (2Q“^/wt^) exp (-E*/RT) (4) 

in regions of rapidly varying aCCenuaclon. In chls paper we use Che seismic 
data CO esclmace Che relaxaCion clme in Che upper mancle and ics varlaClon wlch 
depCh. Derails of Che aCCenuaClon in che midsC of an absorpCion band, and Che 
effecC of a specCrum of relaxaCion Clmes are discussed in a companion paper 
(MlnsCer and Anderson, 1980). The High TemperaCure Background (HTB) is a domlnane 
mechanism of accenuation in crysCalllne solids aC high- CenperaCure and low- 
frequency. This usually satisfies a Q a. u or Q relation. The latter has 

been Interpreted by Anderson and Minster (1978) in terms of a distribution of 
relaxation tines. The former is valid for all relaxation mechanisms when the 
measuremer^t frequency Is sufficiently high* In the following"^ will be used 
for the pre-exponential creep, or Maxwell, relaxation tine and will be tised 
for the attenuation, or Q, characteristic tine. 
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n.b. Activation for oXi»i> iiid gli4t 

Tha craop of ductlla natorlala such aa natala la rata llml^'.ad by aalf- 
dlf fusion. This is tha cssa for althar pura dlffuslonal craap or dislocation 
cllab. In thasa clrcuastancas tha activation anargy fcr craap is tha sama as 
thst for salf-dlf fusion. In ollvlna tha activation anargy for craap is appraciably 
hlghar than that for salf-dlffuaion of 0^" or Sl^"^. This suggasts that craap 
nay ba controllad by kink or Jog formation (Guaguan, 1979]* 

Tha activation anargy for craap uhan dislocation cllnb is due to 
nuclaatlon and lateral drift of Jogs is (1/2) + 2Ej) or (1/2) (E^^ 

lAf A 

Eq) + 4Ej ) depending on whether tha dislocations are longer or shorter than 

A A ^ 

the equilibrium Jog spacing. Here E^^, E^^^ and Ej are the activation energies 
for self-diffusion, core-diffusion and Jog formation, respectively [Hlrth and 
Lotha» 1968]. 

if 

Egjj is 90 kcal/mole for both oigrgen and silicon diffusion [Reddy et, al. , 
1980, Jaoul et al . . 1979]. The Jog formation energy is somewhat greater 
than kink formation energy [Hlrth and Lothe, 1968]. The activation energy 
controlling high-temperature deformation of olivine has been estimated by 
Goetze and Kohlstedt [1973] to be 135 kcal/mole from the clini> of dislocation 
loops in olivine and 122-128 kcal/mole by Ashby and Verrall [1978] and Goetze 
[1978] from creep data on olivine. 
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In th« abs«nc« of Intorstltlal dofocts n«ar th« dlalocation tho 

* * 

activation anargy for gllda la E|^ If tha kink danalty la high and 2 ^ whan 
tha dlalocatlon langth la ahortar than tha aquUlhrlua dlatanca batwaan klnka. 

la tha kink formation anargy which Stockar and Aahby [1973] aatlmata to ba 
26 kcal/mola. Tha activation anargy for gllda la tharafora 52 kcal/nola If 
tha kink apaclng la larga or If doubla kink formation la raqulrad for glide. 

If point dafacta muat dlffuaa with t>«a dlalocatlon llna than thalr dlffualvlty 
may ba rate limiting. Tha activation anargy for Mg^Fa Intardlffualon In 
ollvlna la about 58 kcal/mola for 0.9 mola fraction Mg» dacraaalng to 50 kcal/mola 
for pura foratarlta [Mlaanar, 1974]. Dlalocatlon gllda In ollvlna will thara- 
fora have an activation anargy of tha ;'>rdar of 50-60 kcal/mola unlaa<^, there 
are alower moving apeclaa In tha vicinity of the dlalocatlon. Ollvlna In con- 
tact with pyroxene la likely to have exceaa alllcon Interatltlala aa a dominant 
point defect [Stocker, 1978 ]. The activation energy In thla caae may be due to 
drag of alllcon interstitials, about 90 kcal/mole. There Is also a small term 
representing the binding energy of Interstitials to the dislocation. This is 
ordinarily a few kcal/mole [Van Bueren, 1961] but data are lacking for silicates. 
The kink and jog formation energies in silicates are also highly uncertain and 
the above estimates for olivine are approximate. 

The activation energy for dislocation core diffusion should be similar 

to that for surface diffusion. This, for metals, is usually about 1/2 to 2/3 

of the lattice diffusion value. If this rule applies to olivine the core 

diffusion activation energy would be about 50 to 60 kcal/mole. The hot-pressing 

experiments of Schwenn and Goetze [1978] yield a value of 85 kcal/mole, close 

4+ 

to that for lattice diffusion of 0 and Si . The dislocation mobility experi- 
ments of Goetze and Kohlstedt [1974] also suggest that the activation energies 
for lattice and core diffusion are similar. The actual values determined in 
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ch«lr «xp«rliMnti may ba Intarpratad at tha aua of tha diffusion and Jog foriBa*=* 
tlon anarglaa. A high valua for tha cora and surfaca activation anarglas la 
axpactad from tha haat of vaporization* which la much hlghar for slllcatas 
than for matals. 

If wa adopt E* * E* > 90 kcal/mola and E* • 125 kcal/mola than tha 
Implied jog formation energy la 35 kcal/mola. This la much greater than la 
typical of metals. This la also expected since tha kink and Jog formation 

3 

energies are proportional to Gb which la much greater for silicates than for 
metals . 

Because of the high values for core diffusion and Jog formation we 
expect relatively high values for the creep activation energy for silicates 
and a non-correspondance with the self-dlffusional activation energy* except 
for Coble and Nabarro-Herrlng creep which will dominate for very small grain 
sizes and low-stresses. 

A summary of the physical properties of olivine Is assembled in 
Table 1. These will be used in subsequent calculations. 
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II, c, Gtophyelcul con«tralnts on th« activation param»t«r« 

SoiM of the ectlvetlon parameters for creep and attenuation can be 
estimated directly from the geophysical data. The characteristic time scales 
of telaicatlon processes In the mantle depend on temperature, through the 
activation energy, and the pressure, through the activation volume. Stress 
has an Indirect effect on rela'<atlon times since It controls such parameters 
as dislocation density and subgrsin size. The characteristic frequency of 
atomic processes and processes involving dislocation motions arc typically 
10~^^ to lO'^^ seconds. The characteristic relaxation time for creep in 
the upper numtle, the ratio of viscosity to rigidity. Is about 10^^ seconds. 

For temperatures of 1500-1600 R, appropriate for the upper mantle, the required 
activation energy Is 145-170 fccal/mole. This can be compared with the activation 
energy for creep of olivine 125-165 kcal/mole[Goetze, 1978, Jaoul et al . . 1979] 
and similar values Inferred for the climb of dislocations In olivine [Qoetze, 1978, 
Kohlstedt ct al . . 1976]. 

Judging from the attenuation of surface waves [Anderson and Archambeau, 

1964, Anderson and Hart, 1978] the characteristic time controlling the attenuation 

2 

of seismic waves In the upper mantle Is of order of 10 seconds. This Implies 
an activation energy In the range 90-100 kcal/mole for the above T^. This Is 

close to the activation energy found for diffusion of 0 and SI In olivine 
[Reddy et al. , 1980; Jaoul et al . , 1979]. A higher value for gives a 

it 

smaller value for E . 

The trade-off between and e’*' is shown In Figure 1. The cutves 
to the right cover a range of upper mantle temperatures and viscosities. The 
curves to the left represent the seismic band. Theories of creep and attenuation 

■k 

predict ^ 8 pairs > Those which fall in the areas shown are geophysically 
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plausible mechanisms. Low values of require low values for E , As we will 
show, the characteristic pre-exponentlal times are not the same for creep and 
attenuation. 

It Is not simple to estimate the pre-exponentlal characteristic time 

appropriate for attenuation In the mantle. Hlgh-temperature background damping 

-12 -14 

processes In metals give In the range of 10 - 10 seconds [Nowlck and 

Berry » 1972]. Lower temperature peaks give characteristic times one or two 

orders of magnitude higher [Wolrgard, 1976]. Theoretically, the times scale as 

the dlffuslvlty which Is typically three orders of magnitude lower in silicates 

than In metals. On the other hand, subgrain sizes and dislocation lengths in 

the mantle can be expected to be one or two orders of magnitude greater than 

2 

typical laboratory grain sizes In metals. Since T jt /d the characteristic 

—6 —8 

time for the mantle may therefore be of the order of 10 to 10 seconds. 

There Is very limited data on silicates. For sintered forsterlte with grain 
-4 -13 

size of 5 X 10 cm Jackson [1969J obtained 4 x 10 sec for a relatively low- 

temperature relaxation peak. His temperatures were not high enough to resolve 

the HTB peak but one can calculate from his data that « 10~^ seconds. The 

subgrain size In olivine at 10 bars, a reasonable mantle stress. Is of the 

order of 10 ^ cm [Durham et al . . 1977]. The mantle relaxation times can there- 

-8 

fore be estimated from the forsterlte data as about 10 sec. This Implies an 
activation energy of 58-74 kcal/mole If the upper mantle seismic absorption 
band is centered at 100 seconds. 

2 

Although T Is proportional to I for most attenuation and creep 
mechanisms there are other factors which affect the characteristic tlm.Hs. 
Steady-state creep is controlled by self-diffusion of the slowest moving 

it it 

species. Tills affects both and E . E for creep, which also Includes the 




energy of Jog fometion, le therefore greeter then E** for ettcnuetion which 
need not Involve 8elf->diffuelon. For e polygonlzed network the climb of dis- 
locetions in cell wells is rete-limiting but the cherecteristic time for creep 
Involves both the motion of the long mobile dlslocetions in the cells end the 
shorter dlslocetions in the wells. For ettenuetion, the seperete populetlon 
of dlslocetions would give rise to widely sepereted Internel friction peeks. 

A wide seperatlon of cherecteristic times cen '"««fore be obtelned from e 
single dlsloceted solid. 

Although there is e lerge uncertelnty in the sppropriete for 
ettenuetlon in the mentle the renge of inferred ectlvetlon energies is consistent 
with diffusion control. The lower end of the renge is consistent with ectlvetlon 
energies found at high temperatures, low frequencies and low strains in AI 2 OJ 
and Mg2S10^ [Jackson, 1969]. These in turn ere comparable to those appropriate 
for interstitial drag (Mg-Fe) or double-kink formation. 

e 

II. d. Estimates of t and E for olivine 

o 

The characteristic relaxation time for creep of olivine can be 
obtained from high- temperature' laboratory experiments. The relatively low 
stress, <2 kbar, creep data satisfies 

E - Ao“ exp^- E/Rt) (5) 

Using the Kohlstedt-Goetze data with n >■ 3 and E 125 kcal/mole and the 
relation 

■ oG/e (6) 

o 0 

••IX 

we obtain'?' ■ 3 x 10 to 3 x lO" sec for the stress range 1-10 bars. 

0 

2 *“9 ""XO 

The low-stress data has a trend. In this case 'r' is 10 to 10 sec. 

0 

A summary of 'r^ values consistent with the olivine creep data is given in 
Table 2. A range of activation energies is given. 
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From the rigidity and vlfcoalty the relaxation time of the upper mantle 
la 10^'^"^ seconds. Using an activation energy of 125 kcal/mole and a range of 
upper mantle temperatures from 1100 to 1300*C the Inferred value for t 

o 

Is 10~^ to 10~^^ seconds. This Is consistent with the laboratory creep data 
for differential stresses lu the range of a fraction of a bar to about 5 bars. 

The mantle values also overlap the lower values for ^ determined from 
dislocation mechanisms rate limited by silicon dlf fission. 

The effect of pressure can be estimated from the observation that both 
viscosity and Q do not vary by more than 2 orders of magnitude In the mantle 
[Anderson and Hart, 1978, Peltier, 1979]. This constrains the activation 

3 

volume to be between 4 and 9 cm /mole. Since the activation volume of oscygen, 
the largest major Ion In the mantle. Is about 11 cm' /mole and decreases with 
pressure and since the effective activation volume depends on all 
diffusing Ions for coupled diffusion, the result from the mantle Is consistent 
with dlffuslonal control for both dl.^locatlon creep and attenuation. 

A crude estimate of the activation energy controlling attenuation can 
be obtained from the observation that Q can vary by about an order of magnitude 
over relatively short distances both laterally and vertically In the upper 
mantle. Using 200” C as a reasonable difference In temperature the Inferred 
activation energy Is about 52 kcal/mole. This, In turn, . requires a of about 

10**'^ - 10~° seconds for the attenuation mechanism. 

These results support previous conclusions that both creep and attenuation 

in the mantle are activated processes that are controlled bjf diffusion. They 
are also consistent with numerous studies on a variety of materials that high 
temperature creep and Internal friction are controlled by the diffusive motion 
of dislocations. The large modulus defect Implied by the upper mantle low- 
velocity, low-Q zone is also consistent with a dislocation relaxation mechanism 
[Gueguen and Mercier, 1972; Anderson and Minster, 1980]. Point defect Internal 
friction peaks are generally relatively sharp, occur at low temperatures and 
do not exhibit the low Q's and large velocity dispersion that characterize the 
mantle and the hi^ temperature Internal friction peaks. 
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III, Calculation of creep of olivine 

The climb velocity of a dislocation, when both e#l£-dif fusion and jog 
foraatlon are Important is [Hlrth and Lothe, 1968] 

V - (AirDgjjOb^/kT) exp(Ac. „£.j)/2RT (7) 

We use “SD ■ “st * <“si>o ~Eg^/RT) since Sl^*** appears to be the slowest 
moving species. We assume that the material contains cells or subgrains of 
diameter L. The creep equation Is [Glttus, 1976] 

e - K^^on^^^/G^T (8) 

where a « Gb/L has been taken as the effective stress operating on the mobile 
dislocations of length L. The parameter K Is the ratio of cell diameter to 
the average dislocation length. Including mobile and cell^all dislocations. 
This parameter has not been measured In olivine but the ratio of average mobile 
dislocation length to cell diameter Is about 15 [Durham et al . . 1977]. This 
means that the Glttus parameter K > 15. Glttus also ast. jmes that each svd>graln 
contains only one mobile dislocation. Equation (8) still applies If there are 
more dislocations per subgrain but the K parameter would have a different 
value. The strain rate Increases with an Increase in mobile dislocation 
density but decreases with a decrease in length of the dislocations. K Is 
therefore a slowly varying function of the mobile dislocation density in the 
subgrains. AE in equation (7) :»s Egjj - E^p. 

The calculated creep curves are given In Fig. 2 using for silicon 
and K of 29 and 15. A K of 29 fits the data. There is a direct trade-off 


between and K. If oxygen proves to be the slowest moving species then, 
with current estimates of Its dlffuslvlty, the compatible K Is 8. 

The subgrain structure of stressed olivine, a characteriatic It shares 
with other crystalline materials, also affects the attenuation. If we assume 
that only the free dislocations contribute to the high temperature attenuation, 
their lengths are 


L - K(Gb/o) (9) 

The relaxation time then depends Inversely on stress squared and, with K ■ 10, 
the characteristic times are 10 greater than they would be for a uniform 
Frank network. Note that the characteristic time for creep Is proportional 
to i.e. very much shorter than the case for a uniform network. 


Steady-state creep data for olivine for stresses between several 
hundred bars and several kllobars is well fit by a dislocation climb model 
where the strain rate Is proportional to the third power of the stress [Goetze, 
1978]. At lower stresses the data deviates from this relationship and approaches 
a stress-squared dependency. The stress-squared dependency at low stresses Is 
suggested by data of Kohlstedt and Goetze [1974], Durham et al. [1979] and 
Rerckhomer et al [1979 ] . 

The dislocation climb velocity for the case where the equilibrium Jog 
spacing Is greater than the average dislocation length is [Hlrth and Lothe, 

1968], 

4ir Z Dgjj fi b a 

V « eiip (AE* - 4E*)/2RT (10) 

b kT in(z/b) 

where z is the mean free path of a core vacancy along the dislocation line 
(al.Aa), a is the height of a thermal jog at high temperature, Z is the dis- 
location length, the pre-exponential self-diffusion coefficient (approxi- 

mately equal to the dlffuslvity of the slowest moving species) , 0 is the atomic 
volume, b is the Burgers vector, a is the stress, AE is the difference in 
activation energy between self-diffusion and core diffusion. Setting 
a b 'V and writing e ■ pVb, we have 

e 4irDgjjba^/kTG (11) 

We have used the relation Gb/Z • a where G is the rigidity. Equation 
(11) has the required stress dependency and an effective activation energy which 
is substantially greater than the self-diffusion activation energy. 
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IV. a. Attanuatlon 

The actual phyalcal nechanlam of attenuation In the nantle la uncertain 
but Ic la likely to be a relaxation proceaa Involving a dlatrlbutlon of 
relaxation tliaea (e.g., Mlnater, 1980). Many of the attenuation nechanlano 
that have been Identified In aollda occur at relatively lov teag>eraturea and 
high frequenclea and can therefore be eliminated from conalderation. Iheae 
Include point defect and dialocatlon reaonanca atchaalama which typically give 
abaorptlon peaka at kilohertz and megahertz frequenclea at tenperaturea below 
about half the melting point. The ao«called grain boundary and cold>work peak 
and the **hl^ temperature background" occur at lower frequenclea and hlj^er 
temperatures. These mechanisms involve the stress Induced diffusion of 
dislocations. The Bordonl peak occurs at relatively low temperature In 
metals but may be a higher temperature peak In silicates. 

Even In the laboratory It Is often difficult to identify the mechanism 
of a given absorption peak. The effects of asiplltude. frequency, temperature, 
Irradiation, annealing, deforuLitlon and Impurity content must be studied 
before the mechanism can be Identified with certainty. This Information is 
not available for the mantle or even for the silicates which may be components 
of the mantle. Nevertheless, there is some Information which helps constrain 
the possible mechanism of attenuation in the mantle. 

1. The frequency dependence of Q Is weak over most of the seismic 
band. At frequencies greater than about 1 Hz Q appears to Increase linearly 
with frequency [Kanamorl and Anderson, 1977, Minster, 1978 , Slpkin and Jordan, 
1979] . This Is consistent with the behavior expected on the low-temperature 
side of a relaxation band. A weak frequency dependence is best accomplished 
by invoking a distribution of relaxation times. A distribution of dislocation 
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lengths » grain sizes and activaticn energies nay be involved. 

2. Althou^ it has not been specifically studied » there has been no 
evidence brought forward to suggest that aelomic attenuation is amplitude or 
stress dependent. Laboratory measurements of attenuation are independent of 
amplitude at strains less than 10**^. Strains associated with seismic waves 
are generally much less than this. 

3. The radial and lateral variations of Q are our best clues to 

the effects of temperature and pressure. The lower Q regions of the mantle 
seem to be in those areas where the temperatures are highest. This suggests 
that mo£t. of the upper mantle is on the low- temperature side of an absorption 
band or in the band itself. At a depth of 100 km the temperature of the 
continental lithosphere is about 200 K less than under oceans. Q is 
roughly 7 times larger under continent. This implies an activation energy 
of about 50 kcal/mole. 

4. The variation of Q with depth in the mantle covers a range of less 
than two orders of magnitude. This means that the effects of temperature and 
pressure are relatively modest or that they tend to compensate each other. 

5. Losses in shear are more important than losses in compression. 

This is consistent with stress Induced motion of defects rather than a thermo- 
elastic mechanism or other mechanisms involving bulk dissipation. 


IV. b. Dislocation danplng 


The general expression for the characteristic time of stress induced 
diffusion of a dislocation line is 

T - (RTL^/D^Gb^) exp(E*/RT) (12) 

In the bowed string approximation for climb L is the dislocation length, 

is the pre-exponential self-diffusion coefficient and E is the activation 
energy for self-diffusion. If the climb is rate-limited by core diffusion (CD) 
then [Wolrgard, 1976] 

l} - i^/Xb (13) 

\diere I Is the dislocation length, X b exp^Ej/RT^is the jog separation and 
D^ and E* refer to core diffusion. The effective activation energy is there- 
fore the sum of the core diffusion and Jog formation energies. If D^(CD) - 

10^ D (SD) and i " 10^ b the T for core-diffusion Is 10^ times that for self- 
0 o 

diffusion. At modest temperature, however, t(CD) < t(SD) If E^^ < E^^^. In 
general, relaxation peaks with high activation energies occur at higher 
temperature and lower frequency than those with lower activation 
energies. The presence of point defects, kinks or jogs along the dislocation 
line changes the characteristic time. 

An estimate of the relaxation time for climb can be obtained from the 
values In Table 1. We first assume that the dislocation length Is the same as 
the subgrain size, which In turn Is related to the tectonic stress (Durham et al. 
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Using and O bstwssn 1 and 10 bars wa obtain of 1-10 aaconda. Nalthar 

tha salf-dtffuslon nor tha cora-dlffuslon charactaristic tinas can ba brought 

into tha upper mantle seismic band with the appropriate activation energy 

(Table 1 and Figure 1) . The actual mobile dislocation lengths in olivine are 

.3 

about 15 times smaller than the subgrain size. This reduces to 4 x 10 to 

4 X 10**^ seconds. With an activation energy of 90 kcal/mole the relaxation 

8 6 

time at upper mantle temperatures is 10 to 10 seconds, still well outside the 
seismic band but straddling the Chandler period. Dislocation cllnd> in subgrains 
is therefore unlikely to contribute to seismic attenuation but may contribute 
to damping of the Chandler wobble. 

The dislocation lengths in cell, walls are an order of magnitude smaller 

than in the cells, nils makes a fariner ^’eduction in to about 10~^ to lO**^ sec. 

o 

This gives relaxation times Just outside the seismic band at mantle temperatures 

if 

and E ■ 90 kcal/mole but the relaxation strength is very small because of the 
limited bow-out possible for small dislocations. 

Therefore, it appears that dislocation climb can be ruled out as a 
mechanism for attenuation in the mantle. We therefore turn our attention to 
dislocation glide. 

In general dislocations glide much faster than they climb and the 
relaxation time is therefore much reduced. Glide is rate limited by lattice 
(Pelerls) stresses or by point defects. 

.wr a gliding dislocation, rate limited by the diffusion of interstitial 
defects [Schoeck,1963] 

- C^A^ (15) 

where C^ is the concentration of interstitials along the dislocation line at 
sufficiently high temperature 
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and it tha bulk concantratlon of Intaratltial point dafacta or 
lapurlcy atoaa and la ttia binding anargy of tha point dafacta to tha dla- 
location llna. Tha dlffualvlty la that approprlata for tha dlffualon of tha 
lapurlty or Intaratltial lona« 

Ualng - lO”^ cu^/aac and -v - lO"® va obtain of 10“® to 

-8 

10 aaconda for typical dla location langtha appropriate for nantla atreaaaa 

of 1 to 10 bara. Thaaa, conblned with activation anarglaa approprlata for 

dlffualon of catlona In alllcataa glva relaxation tlsaa In tha aalanlc band 

at upper nantla temparaturaa . Tha Inpurlty content rafara to that In the aub- 

gralna. Moat Impurltlea In tha nantla are probably at grain boundarlaa and 

therefore the grain Intarlora are relatively pure. Tha above aatinata for C la 

D 

therefore not unreaa enable. 

IV. c. Relaxation atrength 

-1 

The relaxation atrength» or naxlnun value for ^ la given by 

2 - (l/eVT 

for a collection of randomly oriented gliding dlalocatlona (Minster and Anderson, 

2 

1980). In the subgrains 1/1 for a rectanglular grid of dislocations 
and the relaxation strength is about 7.SZ or a Q of about 26. The same 
relationship holds for a Frank network In t.he cell walls but because of the 
small volume in the walls the average relaxation strength Is very small for 
the wall contribution to the attenuation. 

The relaxation strength Is proportional to the area swept out by the 
gliding dislocations and gives the difference between the hlgh~frequency or low- 
temperature modulus and the relaxed modulus. Gueguen and Mercler (1972) and 
Anderson and Minster (1980) pointed out that the velocity and attenuation In the 
upper mantle low-velocity zone could be explained by dislocation relaxation. 
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IV, d. Ijpmpltt of dll location wuchanttai of atfnuatlon 

The K0tt«r peak occura in col4-wovkad aatala containing 
Intatafcltlal lapurltlaa, soaa of which have concantratad along the dlilocatlona. 
It occura at high taaparatura and low fraquancy and usually glvaa a ^o 
10“^^ to 10~^^ aaconda for aatals with 0, N or H Intaratltlals . It la bs 'laved 
to ba cauaad by an Interaction betwaan dlilocatlona and Intaratltlal point 
defects In the neighborhood of the dislocation [Schoeck, 1963]. With reasonable 
choices of parameters It can explain attenuation In the mantle. 

The grain boundary peak (GBP) occurs In metals at about one-half 
the melting temperature at 1 Hz. Actually, a series of peaks are often observed 
between about 0.3 and 0.6 The activation anergles range from about 0.5- 
1.0 of the self-diffusion activation energy end are higher for the high 
temperature peaks . An Increase In solute concentration Increases the magnitude 
of the higher temperature peak [Nowlck and Berry, 1972 ]. These peaks are all 
superimposed on an even stronger high- temperature background which dominates 
at very high temperature. In synthetic forsterlte a strong peak occurs at 
about 0.25-0.5 Hz at lOOO^C [Jackson, 1968] with an activation energy of 

57 fccal/molel At upper montle temperatures this absorption would shift to 

-2 

about 10 seconds and therefore would not be laqiortant at seismic frequencies. 

In addition the peak Is much less pronounced In natural olivine and the HTB 
gives greater absorption than the GBP at temperatures greater than about 1500” C* 
The GBP may be responsible for att ;;. lation In the lithosphere. 

The Bordonl Peak was first Identified In deformed ("cold worked") 
metals af low temperature. It Is unlikely to be Important at mantle temperatures 
and seismic frequencies but it would be useful to know Its properties for 
mantle minerals since it contains information about the Pelerls energy. 
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Th« Hlrth and Lotha [1968] thaory for tha charactarlatlc ralaxatlon 
elm give a 


b^kT 

cxp(2E^/RT) 


(17) 


With nominal valuaa for ollvlna tha ralaxatlon tlm at mantla temperaturaa la 
10~^ seconds. Thus , this Is a high frequency, low tamparature mechanism. 

The controlling activation energy Is 2Ej^ or about 52 kcal/nola for 

ik 

olivine. This conblnatlon of and E make It unlikely that the Pelerls' 
energy alone controls attenuation In the mantle. It may be a responsible for 
the "grain boundary" peak In olivine for which Jackson [1968] obtained of 
10“^^ sec and E* of 57 kcal/roole. The value for Ej^, however, is uncertain. 

V. D iscussion 

Values of of order 10~' to 10 seconds for the upper mantle are 

9 10 

implied by relaxation times of 10 - 10 seconds at temperatures of 1400-1600 K 

and an activation energy of 125 ± 5 kcal/mole. These values are consistent with 

both laboratory creep data and creep in a polygonlze'? network model of olivine 

3 

if the stresses are less than about 10 bars. This applies to both the a and 
2 

a laws. If kllobar level stresses existed In the upper mantle then the in- 
ferred viscosity and relaxation tiroes would be at least 6 orders of magnitude 
lower than observed. On the other hand, there is no contradiction with 
kllobar level stresses being maintained In the lithosphere for 10^ years If 
temperatures are less than about 900 K. Kllobar level stresses also shift the 
absorption band to much higher frequencies. 
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The same dislocations contribute to phenomena with quite different 
time scales. Tlie climb of jogged dislocations Is rate-limited by self-dlf fusion 
r.;'d jog-nucleatlon. In silicates > and other materials with high Pelerls energy • 
the activation energy for creep can be appreciably greater than for self- 
diffusion. At finite temperature the creep rate will therefore be slower, 
and the characteristic time longer, than for materials rate-limited by self- 
diffusion alone. The effective activation energy for creep Is either or 
2Ej greater than for self-diffusion, depending on whether the dislocation 
length Is smaller or larger than the equilibrium spacing of thermal jogs. 

These two situations lead to a c or a~ creep law. If the activation energy 
for core-diffusion is less than for self- diffusion the effective, 
activation energy Is reduced by one-half the difference of the two. This is 
a relatively small effect and appears to L«s negligible for olivine. 

On the other hand, the glide of dislocations Is controlled by a much 
smaller activation energy and Is therefore much faster than cllinb at finite 
temperature. This is the case whether glide is controlled by kink nucleatlon 
or diffusion of Interstitials. Both creep and attenuation depend on dislocation 
length and therefore tectonic stress. They are both also exponentially de- 
pendent on temperature although the controlling activation energies are different. 
In principle, the viscosity of the mantle can be estimated from the Q. 





A- 11-24 


VI. Directions for futura raaearch 

There is now abundant data on the hli^-tenperature creep of olivine. 

The various creep theories can be tested with this data and relatively con- 
fident extrapolations can be made to lower attain rates. There la no In- 
consistency between the laboratory and geophysical data. The major taicertalnty 
Is the stress dependence and activation energy at low stresses and the micro- 
structure at low stresses. The Glttus theory satisfies the laboratory and geo- 
physical data. The assumption that there Is only one mobile dislocation per 
grain and that this traverses the grain by glide alone may need modification. 

The theory also -assumes that cell-walls of all orientation behave the same. 
Nevertheless, the Glttus theory seems closer to reality than theories Invol-vlng 
climb alone In Frank networks. 

The calculation of relaxation times and relaxation strengths for glide 
in the subgrains la stral^t forward and appears capable of explaining the 
seismic attenuation data. The laboratory data for checking the theory Is 
basically non-existent for silicates. Low frequency, small strain experiments 
bn deformed and annealed single and polycrystals at a variety of temperatures 
is required. The effect of impurity conte.it must also be studied In both 
olivine and peridotlte. In principal, seismic data can be used to Infer 
temperature, dislocation density, stress and subgrain ln|>urlty content. This 
paper has hopefully made the case for a new generation of geophysical experiments. 
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Table 1 

Material Properties for Olivine 


Property 

Syniiol 

Value 

Units 

Burger's vector 

b 

6 X 10"® 

cm 

Oxygen ion volume 

fl 

1 X 10"^® 

cm® 

Shear modulus 

G 

00 

K 

O 

dy/cm^ 

Silicon diffusivity^^^ 

pre-exponent lal 

“si 

1.5 X 10"® 

2 

cm /sec 

activation energy 

®Si 

90 

kced/mole 

(2) 

Oxygen diffusivity' ' 

pre-exponential 

D 

ox 

3.5 X 10"® 

2/ 

cm /sec 

activation energy 

E 

ox 

89 

kcal/mole 

Mg-Fe diffusivity^^^ 

pre-exponential 

°Mg 

3.4 X 10 ® 

2/ 

cm /sec 

activation energy 

^g 

47 

kcal/mole 

Sub grain size f 

Dislocation length 

K' 

15 

— 

Kink energy 


26 

kcal/mole 

(5) 

Jog energy 


35 

kcal/mole 


^^^Jaoul et al., 1979 
Reddy et al. , 1980 
^^^rfisener, 1974 

Stocker and Ashby, 1973 
^^^From crs?ep data 
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Table 2 

Summary of Experimental Values 
Consistent with Kohls tedt-Goetze Creep Data 



0^-law 

1 bar 

Stress 

10 bar 


120 

3 X 10~® sac. 

3 X 10”^° sec 

B 

(kcal/mole) 

125 

0.7 X 10“® 

0.7 X 10“^° 


135 

0.35 X lO"^ 

0.35 X lO”^^ 


2 

a -law 




120 

3 X lO"^^ 

-12 

3 X 10 


125 

0.7 X 10"^^ 

0.7 X lO"^^ 


135 

0.35 X 10~^^ 

0.35 X lO"^® 


125 kcal/mole Is experimental activation energy 
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Figure Caption 

Figure 1. Pre-exponentlal t vs, activation energy cond>lnatlons that 
satisfy the aelsnlc data (left) and rebound data (right) for a range 
of upper mantle temperatures. 

Figure 2. Creep in polygonlzed olivine for several values of K. data from 
Kohlstedt et al. (1976). 
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Abstract 

The dislocation mlcrostructure of mantle materials can account 
simultaneously for long term steady-’State creep, and for stress wave 
attenuation at seismic frequencies. The hypothesis that a single micro- 
structural model explains the rheology for characteristic times ranging 
from 1 to 10^® seconds can be used to restrict the class of permissible 
rheological models for the mantle. We review steady-state dislocation 
damping models in order of increasing complexity, and reject those which do 
not satisfy laboratory data or geophysical constraints. 

This elimination procedure leads •ir' to consider an organized micro- 
structure, in which most dislocations are found inside subgrain w,. 

The cells contain relatively few dislocation links. These are free to bow 
under small, i.e. seismic, stresses. The time constant of this mechanism 
is controlled either by the diffusion of kinks or of point defects bound to 
the dislocation line. The glide of intragrain dislocations explains the 
magnitude and frequency range of seismic attenuation. Steady state creep 
is governed by recovery through climb and annihilation in cell walls. Under 
conditions of jog undersaturation, climb is controlled by jog formation 

addition to self diffusion, and the model requires a higher creep activation energy 
than for self-diffusion, in agreement with observations on olivine. Quantitative 
agreement with laboratory data is achieved if the density of cell wall dislocations 
is one to two orders of magnitude higher than the density of Intracell dislocations. 
Self'»dlf fusion is probably controlled by silicon diffusion at low pressure and 
by oxygen diffusion at high pressure. The long term tectonic stress is the 
dominant factor which determines scale lengths; as a result, the total strength 
of the relaxation associated with bowing of intracell dislocation links is 
fixed by the geometry and is of the order of 10%. This limits the width of 


ill 


Che seismic absorption band Co 2-3 decades in frequency for each mantle mineral. 
The actual position of the seismic absorption band Is determined primarily as 
a result of a tradeoff between temperature, pressure and tectonic stress. 

This model provides a physical framework within which Q and viscosity 
are related via the dislocation microstructure. 


tv 


(ilossnry of Symbols 


a 

b 

3 . 


C 

d 

D, 

*’oQ 

‘'os 

D(T) 

y(<»>) 

A 

E 


m 


< E < E 


E 


M 


0 


AE* 

• 

r, e 

J(t) 

J 

u 

J 

r 

6J 


K, K , K 
in s 


exponent foi. In o. *ancy dependent 0 

llurRers vector 

geometrical constant 

Jog concentration 

normoli zntion constant 

climb distance for annihilation 

diffusivity and preexponential diffusion constant 

diffusion constant for attenuation 

diffusion constant for self-diffusion 

retardation spectrum 

attenuation coefficient 

relaxation strength 

activation energy and range 

activation energy for attenuation 

activation energy for self-diffusion 

.Jog formation energy 

kink formation energy 

difference between bulk and core diffusion activation energy 

strain and strain rate 

compliance response 

unrelaxed compliance 

re 1 a xo d comp I i ance 

c omp I i an (>e de f e c t 

ml crostructural scaling constants 


V 


k 

k(a») 

t 

L, L 
£ 

c 

A 

V 

A 

U 

n 

w 

0) 

n 


V 

I>(£), p(E*) 
P 

'l'(t) 

Q(u) 

R 


r 

w 



Boltr.mmiii eonatont 

complex wave number 

dislocation link length » and average 

Biibgrain alze^ and average 

critical length for multiplication 

scale length for dislocation drag mechanism 

scale length for vacancy diffusion 

mean free path of dislocations 

rigidity 

number of lattice spacings between cell wall dislocations 

freq iiency 

atomic volume 

vacancy volume 

probability distributions 

pressure 

creep function 

quality factor 

gas constant 

radius of curvature of dislocation line 

dislocation spacing in cell walls 

inner and outer cutoff radii for vacancy diffusion 

total (smeared) dislocation density 

density of nK>bile dislocations 

density of dislocations composing cell boundaries 


0 


s 


T 



T 


V. V 


D 


V (u>) 

P 

V , V 
u’ r 

* 

V 

z 


seismic (transient) stress 
tectonic (steady state) stress 
temperature 

time for multiplication 
retardation time 
preexponential retardation time 
limits of retardation spectrum 
Maxwell time 

dislocation migration (drift) velocity 
phase velocity 

unrelaxed and relaxed phase velocities 
activation volume 

mean free path of vacancy along dislocation core 


Subscript index 

NH Nabarro-Herring diffusion 

IG Intragraln recovery mechanism 

pu dislocation pile-up mechanism 

cw cell-wall recovery mechanism 

2 

sq o creep law 


INTRODUCTION 


It is probable that steady-itate creep in the mantle is controlled by 
dislocation climb (Weertman, 1970, 1972; Goetze and Brace, 1972; Green and 
Radcliffe, 1972; Goetze and Kohlstedt, 1973) and is therefore limited by the 
slow process of self-diffusion and, possibly, jog nucleation. On the other 
hand, cl f -'.location motion at short times and under very low stress is probably 
due to glide and is rate limited by faster dlffusional processes involving 
small point defects or kinks (e.g., Simpson and Sosin, 1972). This has led Gueguen 
and Herder (1973), and Anderson and Minster (1980a) to propose diffusion-controlled 
dislocation bowing as an attenuation mechanism for seismic waves in the mantle. 

There is therefore a possible close connection between seismic attenuation and 
viscosity via the dislocation structure of the mantle. 

In this paper we develop a simplified physical model which allows us to 
bridge the vastly different time scales and strains involved in steady state 
creep and anelastic behavior. The fundamental assumption is that the micro- 
structure of mantle grains is controlled by long-term processes and is there- 
fore statistically time invariant, and that mobile dislocation links are free 
to bow by glide under small transient stresses. The long-term, large strain 
deformation is controlled by climb, and therefore by self-diffusion and jog 
formation. The small anelastic strain associated with bowing of mobile dis- 
location links leads to stress relaxation and consequently to attenuation of 
stress waves. Although greatly simplified, the formalism is quite general 
and can accommodate a distribution of grain sizes, dislocation link lengths 
and activation energies. Tlte detailed analysis will be 

carried out for a physically plausible and geophysically interesting family 
of such distributions. The necessary scaling relations are given in section I. 




Vie then develop a formal theory of attenuation by diffuslon-controlled 
dislocation bowing » and provide the r^ecessary results to analyze absorption 
bands with a mild frequency dependence of Q. In section III» we present 
a sequence of creep models of Increasing complexity » and recast all results In 
terms of Maxwell characteristic times. This permits us to test the various 
models against quantitative observational constraints as well as to impose the 
necessary constraints for self consistency. This Is done In section IV 
for creep models as well as attenuation models. We show that quantitative 
agreement with laboratory data and geophysical observations can be achieved 
simultaneously for both aspects of the riieology If the microstructure presents 
a fairly high degree of organization. The preferred creep model is a generali- 
zation of Glttus' (1976a) model. The width of the absorption band Is probably 
associated with a distribution of activation energies and/or dislocation 
link lengths. 

I. SCALING LAWS 

Although there is abundant evidence that mantle minerals contain many 
dislocations (e.g,, Raleigh and Kirby, 1970; Green and Radcliffe, 1972; 

Goetze and Kohlstedt, 1973), the actual microstructure under mantle conditions 

is uncertain. Plausible assumptions can be made, however, since the bulk of 

2 

mantle material has remained under quasi-constant low stress (1-10 bars) 
over long times (>10 years)# In view of the high temperatures and pressures^ 
we can assume that the dislocation microstructure of mantle minerals ^ situ 
is that pertinent to high temperature steady state creep. 

The substructural characteristics associated with high temperature steady 
state creep have been observed for a wide variety of materials. Summaries of 
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fcha most significant properties show that they are remarkably slsdlar from 
material to material (e.g. , Muknerjee et al . , 1969; Poirier^ 1976; Takeuchi 
and Argon. 1976a; Durham et al>. 1977). By comparison with laboratory results, 
and in view of the supporting theories t>f Gittus (1976a. b). we shall assume 
that most dislocations in the mantle are associated with cell walls of sub- 
grains. and that the subgrains contain only a few dislocations, with no tangles 
nor pile-ups. We shall further assume that these few free dislocations may be 
strongly pinned at isolated points, and that the segments between pinners can 
bow-out under low stress. Under larger stress multiplication may take place, 
either by glide (Frank-Read sources) or climb (Bardeen-Herrlng sources) . 

The problem of superposing a low level transient (seismic) stress on a 
much larger steady (tectonic) stress is quite complex. Because of the very 
different time scales involved, the problems can be decoupled to first order 
and we shall assume that the transient is actually superposed on a static, 
"equilibrium'' configuration of internal stress. 

Further specification and simplification of the model can be achieved by 
use of empirical scaling laws; namely (e.g., Mukher.jee et al. 1969; Takeuchi 
and Argon, 1976a); 

1) The average subgrain size is Inversely proportion to the (tectonic) 
stress 

L/b « m/o^, (1) 

where b is the Burgers vector and p the rigidity. 

2) The dislocation density within subgrains is proportional to the square 
of stress 

b^Pjij « (o^/u)^ (2) 

Here we have Identified intragrain dislocations with the mobile dislocations, 
of density p^. Comparison of these two scaling laws leads to 
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b2p - K 2 (b/L)2 (3) 

III IP 

Tlie observations of Durham at al^. (1977) on laboratory deformed olivine 
lead to a proportionality constant K of the order of lOt although this value l8 

IP 

quite uncertain. Since the volume of subgrains is proportional to this 
implies that the length of dislocation line per cell is proportional to the 
cell dimension L. We may assume that this statistical scaling law holds 
locally in the sense that it holds in a neighborhood which is very large com- 
pared with the cell size but very small compared with the averaging length of 
a seismic wave. 

If p is the total dislocation density, then we may write 

p ■ p + p (4) 

where p^ is the density of dislocations composing cell walls. For materials 
which show a strong degree of organization, p^ may be over an order of magnitude 
greater than p^ (e.g. , Takeuchi and Argon, (1976a) , and thus Is a good approxi- 
mation to the "smeared" dislocation density, in the sense of Holt (1970) and 
Glttus (1976a, b). On the basis of theoretical and observational consideration, 
these authors argue that a relation similar to (3) holds for the smeared density, 
in the form 


b^p ~ (b/L)2 (5) 

Observations on metals point to values of K ranging from 10 to 20, but 

higher values may be expected theoretically for materials with a higher Pelerls 

energy (Glttus, 1976 a,b). Consistency between the three equations (3)-(5) 

requires that p and p scale with grain size, and hence with stress in a 
IP w 
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similar way. Further, is an order of magnitude larger than if K 
is only three times greater than K . 

Although roost of the preceding considerations are based on observational 
evidence on metals, the observations of Durham and Goetze (1977) and Durham 
et al . (1977) indicate olivine is no exception. Of course, detailed description of 
rhe mlcrostructure involves many complications, so that the present model must 
be taken as an idealization. Our purpose is to test some of its consequences. 

II . DIFFUSION-CONTROLLED DISLOCATIOH DAMPING 
1) Retardation spectrum 

The most commonly -* and roost successfully ~ used model of dislocation 
damping is the classical Koehler-Granato-LUcke model (Koehler, 1952; Granato 
and LUcke, ri^56). This is a string model for a dislocation strongly pinned at 
both ends, .wd we use it as a starting point. A recent treatment of the 
problem at low frequency, where Inertial effects can be neglected, is given by 
Simpson and Sosin (1972). 

We may write the force balance equation for the string model of the dis- 
location line in the form! 

Ob- Mb2/2R - VkT/AD - 0 (6) 

8 C 

Here Og is the applied (seismic) stress, R^ the radius of curvature, and V the 
migration velocity (fig. 1). X and D are a scale length and a diffusion 
coefficient which are characteristic of the particular viscous drag mechanism 
considered. Since we treat the problem as the superposition of a transient 
(seismic) stress on an equilibrium state, R^ and V take the values which 
correspjnd to equilibrium with o^. For small bow-out (6) reduces to a diffusion 
equation, and solutions have been derived by Friedel ^ al . (1955), Weertroann 
(1955), Schoeck (1963), and numerous subsequent authors. 
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For a link of longth It strongly pinnsd at both anda, tha analaatlc 
strain dua to bow~out Is givan by 

e(t) « c(«) (l~oxp(~t/T) 1 

which is tha classical rasult for a standard linear solid 
tardation time is 

T - af^icT/w^pb^AD (8) 


(7) 

where the re** 


If wa introduce the steady-state drift velocity of an infinite dislocation 

line under the applied stress o , 

& 


T 


i iL 


(9) 


which is also of the general form 


T ■ exp[EQ*/RT] 
. A(T) 


( 10 ) 


A 

Here la the activation energy for tlds attenuation mechanism, and D^q is 
the preexponential diffuslvlty factor for the controlling diffusing species. 

A distribution of retardation times can be achieved through a spectrum 
of preexponential factors and/or a spectrum of activation energies . It 
is well known (e.g.. Nowick and Berry, 1972) that a localized distribution of 
retardation times, with a spectrum decaying rapidly away from a dominant value, 
can be replaced for most purposes by an effective, single mechanism. In such 
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casM, tha aaaanCial diaractariatica of cha nadiuM do not dlffar draatlcally 
from thoaa of a atandard linaar aolid. 

It iR more uaaful, for our purpoaea, to allow for the poaalblllty that a 
broad apectrum of retardation timaa la available. Wa have auggaatad pravloualy 

(Anderaon and Mina t«r, 1960a » Hinster and Andaraon, 1980) a normalized danalty of tha 
form 

D(t) - H(t-t ) (11) 

Among the many possible physical models which entail the distribution (11) we 
shall consider two extreme cases: 

1) A distribution of lengths with probability density 

P(^) - £*“■*'* { t < (12) 

where is the normalization factor. 

2) A distribution of activation energies with probability density 

p(E*) - Cg exp (aE*) ; S* < E* < eJ (13) 

where Cg is again a normalization factor. 

The general case clearly consists of a simultaneous distribution of lengths 
and activation energies. In his discussion of the problem, Macdonald (1963) also 
considers juxtaposed intervals in t with different values of o. 

As m shall see shortly, the nonlocal character of the absorption band 
associated with the spectrum (11) is only significant for small values of a. 

The richest variety of rheological behavior is observed in the range 
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a € l-l.Xj} in th« remainder of this paper, we confine our attention to this 

range. Note, however, that the limit a o must be taken with some care, and 

results for this case will be given separately. 

Tlie cutoffs t and t„ are taken here as phenomenological parameters to 
m Pi 

be constrained by the physical model. Since the actual behavior of the spectrum 

outside the range It ,t„J is immaterial for our purposes, as long as D(t) decays 

m M 

rapidly away from this range, we have assumed sharp cutoffs. Localised spectra 
can be conveniently simulated by taking a narrow range, t^^ Tj^, irrespective 
of a. 

Let Jy be the unrelaxed compliance of the material, and 5J be the 
compliance defect. The compliance response is then 

J(t) - - Jy(l+Ai|/(t)] (lA) 

Here ijj(t) is the normalized creep function and A the relaxation strength. The 
medium is anelastic, that is, undergoes time dependent, but linear and recoverable 
strain response to an applied stress (Nowick and Berry, 1972; Minster, 1979). 

The complete anelastic behavior of Che material under low stress can be described 
in terms of A and ijj. 

2) Relaxation strength 

The compliance defect iSj depends on the total area swept by mobile dis- 
locations as they bow out, as t ■> «. The calculation is a classical one (e.g., 
Friedel al,. , 1955) and yields 

A - 0d T f* P(i) dl (15) 

I 'o 

Here is the density of mobile dislocations, p(£) the density function of the 
link length distribution, and J the average link length. 3^ is a geometrical 
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factor which depends on the orientation and geometrical arrangement of the 
dislocations. For edge dislocations which may only glide, and are all optimally 
oriented 1/6. If they arc randomly oriented drops to l/6i/5. If they are 
allowed to bow by climb as well as glide, 3^“ /7/6i/i5. 

actually, the averaging should only take place over the possible orienta- 
tions of the glide planes (Friedel et , 1955; Simpson and Sosln, 1972). For 
our purposes, approximate estimates are sufficient. Evaluation of (IS) leads 
to the following results: 

a) Case of a distribution of activation energies 

A ~ 3d ^m 

It may be noted tb.;t this result is the correct one whenever the distribution 
of lengths is highly localized. Whenever the separation of mobile dislocations 
is comparable to the length 1 , (e.g. , Frank network), the relaxation strength 
is about c% for glide, and 11% if climb Is permitted as well. 

b) Case of a broad distribution of lengths, 

We find in that case if a is sm;; i.l 


I 


3— 2ot n 
2- 2a m 


and we may distinguish three cases, according to the sign of a: 


(17) 


a < o 


A 


(1-g)^ 

a(3-2o)^ 




A 



I 


m 


d m 


(18) 


a 


o 


(19) 
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a > o 


A * 


(3-2a) 



cl m 


( 20 ) 


Comparison witli (12) and (17) indicates that the relaxation strength is dominated 
by the more numerous shorter links if a is small, and that it Increases rapidly 
with a. This is in agreement with expectations since an increase in a leads to 
relatively more abundant long dislocations. Note that in equations (16)-(20) 
we can use the scaling low (3) in the form 


0 




( 21 ) 


which illustrates more precisely tlie effect of the microstructure. If the link 
length scales linearly with subgrain size, then a spread in the distribution of 
lengths tends to increase A toward values greater than the 10!S! predicted by (16). 
If the link lengtii f does not increase with L, then the material tends to be 
underpopulated in mobile dislocations and the modulus defect is rather small. 

3) Creep functi on 

The normalized creep function is given by 


i/»(t) 




D(r)dT 


Wltl) the density (11), evaluation of the integral yields 


( 22 ) 


|^(t) 


at 


a 




a 

T 

m 



(23) 


where T is tiu* Incomplete gamma function. This expression possesses a proper 
limit as IX ^ o, which can be expressed in terms of exponential integral functions: 
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♦ ^(t) « I “ r~i — k~T U* ~ I ~ r ~l <2^‘) 

Vh, [ \ ^»/ \ 'm/J 

III the limit cC a narrow npectrum (i ~ t„), we get the case of a standard 

in n 

linear solid, which does not require elaboration. In the case of a broad 

spectrum (t « t„) , the creep history is best characterized in terms of three 
m H 

successive regimes. 

a) Short time t « t_ 

in 

For such short times, none of the dislocation segments has relaxed. 

We find 


a < o : i|»(t) ~ 

-g t 

l~a T 


(25) 


m 



a * o : i/t(t) ~ 


JL, 

T 

m 

(26) 


^ / r 



a > o : i|»(t) ~ 

g / Mn 

l-a Itm 

/ m 

(27) 


This linear time dependence implies Newtonian behavior for small time. It 
is of interest that the normalized creep Function i|( depends on the width of the 
spectrum in (26) and (27), hut does not in (25), for which the density of long 
retardation times Is much smaller. If we consider the combination Aij;, this 
remains true in the case of a distribution of activation energies (eq, 16). 

On the other hand, combination of eqs. (18-20), and (25-27) show that irres- 
pective of u, Aiji depends only on for small times if the spectrum is due to 
a length distribution. This late result could have been anticipated since only 
the shorter links can react on this time scale. 

b) Transitional Regime t « t « r.. 

During this interval, more and more dislocation links reach their final 
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equilibrium position. As a result we expect this regime to be one of strain 
hardening. The dominant terms are, according to the sign of a; 

( 28 ) 

(29) 

(30) 

These are reminiscent of various transient rheologies which have been observed 

on many materials. Once more, we note that Aij» is independent of the width of 

the spectrum if (18)-(20) are used, as opposed to (16). And again we note the 

influence of a on the rheology. Kor a negative, strain hardening is very rapid 

and the anelastic strain rapidly approaches an equilibrium value. For a positive, 

the higher density of slow retardation mechanisms is felt for times much longer 

than T and anelastic strain continues to grow as a power law. The intermediate 
m 

value « = o yields a logarithmic creep law. 

c) "L ong time” regime « t 

For such long times, all links have relaxed to their equilibrium, bowed 
position. The normalized creep function is exponentially close to unity in all 
cases. This stage Is reached when all dislocation segments have reached the 
same equilibrium radius of curvature. This radius may be smaller than 0.5 
however, for high enough applied stress. In that case, an equilibrium con- 
figuration is never reached because the longer links begin to multiply and 

this gives rise to unbounded plastic strain. We shall return to this point 


a < o : i|»(t) ~ 1 - r(l-a) 

“ - « = *<'=) ~ [<-■ * t 

where C is Euler's constant 



later. 
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4 ) Attenuation propertie s 

In the frequency Uomain, the anelastic properties described above give 
rise to an absorption band with fairly simple properties. The basic concepts 
are found In Liu e^ al. (1976), and have been recently reviewed by Kanamocl 
and Anderson (1977) and Minster (1979). Let be the unrelaxed high-frequency 
wave velocity; the complex wave number is then given by 


k2(u») 



1 + 2^ r^M _-jl! 

■^2 - / t*“®(l+iWT) 

M III J 

m 


(31) 


which may be represented In terms of hypergeometric functions. As sho\ a by 
O'Connell and Budiansky (1978) an unequivocal and convenient definition of the 
quality factor is 


Q(ui) »• - Re [k^(u))j / Till [k^(w)] 


(32) 


Although no simple analytical form can be derived for Q(w) , asymptotic expressions 
can be obtained fairly easily in the various regimes already mentioned. In the 
limit of a very narrow spectrum, the absorption bund reduces to the classical 
Debye peak. For a broad spectrum, we consider three frequency bands, 
a) Long periods wt « aiT.. 1 

II) 


a < o ; Q(w) ^ - 


l-fg HA 
a, A 


(“I (WT)^) 

\Tm/ ” 


a = o : Q(u)) 


HA , 

n “ 1 


<»T„) 


-1 


u > o ; Q(w) 


('oTm) 


(33) 


(34) 




1+0 

a 


i+A 

A 


(35) 
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Q(b») is found to increase linearly with period, the classical behavior of the 
long period side of an absorption band. For monochromatic waves at such 
periods, all dislocation links bow out to equilibrium and the medium exhibits 
relaxed elastic behavior, with phase velocity 


V V / /t+A 
r u 


(36) 


b) Absorption band « 1 << u»Tj^ 

Asymptotic approximations are more laborious to derive in this case. We 

find 


a < 0 


Q(m) ''' cotan y- 


1+A cos oif/2 / ^a 
-T ~W2 <“^m^ 


(37) 


a “ o : 


a > 0 ; 


. 1+A 2 . M 

m 


- inwT,. 
W M 


V . ^ aiT . 1 cos onf/2 , _ 

Q(o)) <v cotan - 5 - + j (u.t„) 


(38) 

(39) 


Depending on the sign of a, Q(d)) decreases or increases as w®. With a = o 
Q is quasi-independent of frequency; this assumption is often used in seismology 
(e.g. , Liu £t al. , 1976; Anderson et al. , 1976; Minster, 1978a, b). A mild fre- 
quency dependence such as predicted by (39) has also been suggested by various 
authors (e.g., Jeffreys, 1970; Jeffreys and Crampln, 1970; Macdonald, 1961, 1963; 
Anderson and Minster, 1980a, b; Strick, 1970). 

It is interesting to note that there is little generality to be gained 
by letting a take values outside the range 1-1,11. A power law dependence of 
Q(o)) on (D still holds, but Q(w) cannot vary faster than u» or w ^ within the 
framework of linear models such as the present one. Thus, a sharp drop of the 
relaxation spectrum is difficult to distinguish from a simple cut-off Insofar 


as the absorption band is concerned. 
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The physical dispersion Is obtained from the definition 


k(o)) » YTvi) " 

P 

where V (o)) is the phase velocity and v(a)) the attenuation coefficient. 

P 

Asymptotic results are conveniently derived only if Q(u) » l.e., in the 

low loss approximation. We have then 


a < o : 


a “ o : 


a > o ; 


V (w) ~ V 
P u 


1 A . 1 ^ on /%■ 

1 - j + 2 cotan -j Q 


^(w)! 


Vp(u)) ~ j^l + (ninujT^)"^ Q”^(u)j 

Vp(o)) ~ 1^1 - ~ cotan ^ Q”^(aj)j 


(4p 

(42) 

(43) 


Thus, for a 4 the phase velocity varies as a power of frequency, whereas 
a * o yields the classical logarithmic dispersion, 
c) HlRh Frequencies 1 « « on],) 

At high frequencies, none of the dislocation links reach their relaxed 
equilibrium configuration. Little anelastic strain takes place, and we have the 
following limiting behavior: 


a < o : 


a = o : 


a > o : 



(44) 

(45) 

(46) 


Once again, comparison with (17)-(20) shows that these limits are Independent of 
the distribution of lengths when such a distribution controls the absorption band. 
This is consistent with our discussion of i/>(t) for short times. On the other 
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hand, the width of the band enturu (45) and (4b) explicitly when one considers 
a distribution of activation enetRies instead. The linear increase of Q with 
frequency is characteristic of the hiRh frequency side of an obsorption band, 
Vp(u>) asymptotically approaches in that xlmit. 

These properties are summarir.ed on figure 2 which depicts three absorption 
bands with a ■ -0.25, 0, and 0.25 respectively, as well as the associated dis- 
persive characteristics of the medium. 

The so-called high temperature background Internal friction in crystalline 
solids possesses the general characteristics derived above for intermediate 
and high frequencies (e.g.. Minster, 1979). There is therefore a suggestion 
that laboratory frequencies are on the high frequency-low temperature side of 
an absorption band of the type discussed here. A similar relationship appears 
to hold in the upper mantle at seismic, frequencies (Anderson and Minster , 1980b) . 
Since an investigation of attenuation as a function of temperature depends 
critically on the activation energy, we defer discussion of these aspects until 
a specific physical model is discussed. 

II 1 . STEADY STATE RHEOLOGY 
1) An Idealized micros tructural mod el. 

For a given applied (tectonic) stress o,j, the equilibrium radius of 
curvature of the dislocation lines is obtained by solving (6) in the static 
limit 

= wb / 2o,j, (47) 


liinks with lengths f '• = 2R^ bow beyond a serai-circle and thereafter multiply, 
either through a Frank-Read mechanism, if glide only is permitted, or a Bardeen- 
Herrlng median ism for glide plus climb. A crude estimate of the time required 
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Cor multiplicacion ia obtained by conaidering the mid-point of the link, and 
assuming that the link Is a circular arc at all times during bowing out; then 

r ^/2 -1 

t^a) -/ V \x) dx (48) 

"^o 

where V(x) is obtained by salving (6) for this geometry. The solution takes 
the form 

t, 

til) « ^ f(£/f ) (49) 

M c 

where 


f(x) 


1+ Ln 






(50) 


Figure 3 is a graphical comparison of the relaxation time (9) and the multipli- 
cation time (49) under the same applied stress a •* Og » 0,^,. It illustrates 
the region of validity of the anelastic theory developed in the previous section. 
For sufficiently large stress, is shorter than some or all of the mobile 
links. These links multiply on a time scale and we have non-recoverable 

plastic strain. The function C(x) has a minimum for x » 1.66; thus links with 
Z = 1.66 Z^ are the first to multiply. It is reasonable to expect that they 
will dominate the distribution of lengths in the long run since they are associ- 
ated with the most efficient yielding process. According to this argument, we 
expect 

Z ~ 1.66 pb/fJ^, ~ Pjjj ^ (51) 

The remarkable agreement with the observations of Durham et al. (1977) is 
undoubtedly fortuitous, but nevertheless gives some credence to this simple 
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argument. It is noteworthy that this only differs slightly from the 

results of takeuchi and Argon (1976b) . which involve much more sophisticated 
modeling techniques. 

Observations of subgrain dimensions in olivine compiled by Durham et al . 
(1977) are well fit by a relation of the form 




a o„ 


(52) 


with K in Che range 15 to 50. When combined with (51) and (3), this yields 

a constant in the range 9 to 30. Under mantle conditions we may speculate 

Chat the material has been annealed over a very long period of time, so that 

intragrain dislocation density may be somewhat lower Chan in laboratory samples, 

and Che degree of organization higher. 'Hiese factors conspire to decrease Che 

estimate of K^, so that 10 may be an appropriate value, especially in view 

of the uncertainties involved, Ttius emerges a micros Cructural model involving 

subgrains of dimension L ~ 20 containing rather few mobile dislocations of 

mean length 1 ~ 2f.^. All lengths scale inversely with the lon<s term tectonic 

stress a,„. If most dislocations are to be found in cell walls so that p in 
T w 

equation (4) is Indeed over an order of magnitude greater than p^, then we must 
take the constant K in equation (5) to be greater than 30. We now explore the 
consequences of this highly idealized and simplified microstructure for various 
creep mechanisms. 

2) Creep mechanisms . 

A convenient method of comparison of rheological mechanisms is through 
their characteristic times. For steady state creep we consider the Maxwell time 


T - 0^. / lie 


(53) 
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where e le the attain rate. Reviews of posalble creep nechaniaas are pifovided 
by Weertmann (1970), Stocker and Ashby (1973) and Poirier (1976). For our 
purposes we shall compare three basic mechanisms: diffusional creep, dislocation 

creep controlled by intragrain recovery, and dislocation creep controlled by 
cell wall recovery. 

For diffusional creep we use the Nabarro-Herring creep equation in the fora 

0 o_« 

^NH * ^NH T? la <5A) 


so that 



kTb^ exp [E*/RT] 

p D U 
os 


(55) 


Here we have assumed the cell walls to be efficient sources and sinks of 

it 

vacancies. D , and E are the preexponential diffusivity and activation energy 

OS s 

for self diffusion, and U is the atomic volume. is a constant of order 20 

Nn 

(Stocker and Ashby, 1973). This equation is easily generalized to include Coble 
creep. 

For dislocation creep models we invoke the simple form of Orowan's 
equation discussed by Poirier (1976). At steady state, the hardening rate is 
exactly balanced by the recovery rate and the strain rate is expressed by 


e 


P h 
m 



(56) 


Here A is the mean free path of the dislocation; d is the distance covered at 
the rate controlling speed V. Most models assume that V - V^, the climb velocity, 
and d is the distance covered by climb. 


et 
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Tak«uchl and Argon (1976b) have proposad a modal whara glida and climb 
are equally difficult, due to the drag of a Cottrell atmoaphere by edge Begnents. 
In that case the mean free path A la about equal to the critical length t and 
so is d. An upper bound to the creep rate, and thus a lower bound to the Maxwell 
time is obtained by taking V »• V^, and we may use (e.g., Hirth and Lothe, 1968) 

V . 2, Hg-v V 

O IC Jl ^ V 

Here is the volume of a vacancy, and is a scale length associated with 
vacancy diffusion, of the form 


X^ b in(r 2 /rj) 


(58) 


where r^^ and r 2 are inner and outer cutoff radii (e.g., Hirth and Lothe, 1968; 
Takeuchi and Argon, 1976b; Poirier, 1976). Typical values of the ratio 

4 

for intragrain processes are of the order of 10 , and the Maxwell time is found 
to be 



kT b2 exp (Eg/RT] 
“v 


(59) 


which is of the same form as (55), but where the constant is now of order 

unity. This mechanism is thus more efficient than the generalized Nabarro- 

d ^ 

Herring diffusion. The estimate (59) may be decreased further ~ "4 1-^ 

glide is somewhat easier Chan climb and the mean free path of intracell dislocations 
is longer than £. But in chut case we rapidly reach the point where many, if 
not most dislocations actually collide with cell walls. Then a more appropriate 
model is one where recovery in subgrain boundaries plays a specific role. Such 
models have been reviewed by Weertmun (1972) and by Poirier (1976), and a brief 
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(liBcusBion Is given by Minster and Anderson (1980). For W«ert«an’s (1972) nodel 

A 

with dislocation pilo-ups, the Maxwell time r is given by (59) where the 
proportionality constant is replaced by 

C - 1/0.006K ^ (60) 

pU 8 

Alternative models have been proposed v^hich yield a stronger stress dependence 
of the creep rate. However, as noted by Weertsian (1972) and Poirier (1976) they 

3 

usually stem from ad-hoc modifications of the o creep law which arise naturally 
under minimal assumptions. A model which Involves explicitly the microstructural 
constant K is that of Gittus (1976a). The following is a simplified treatment of 
it, adapted for our purposes. 

Tl»e mean free path of mobile dislocations is bounded from above by the 

subgrain size L. In the limiting case where recovery takes place exclusively 

by climb in cell walls, we actually take A « L. The mean separation between 

- 2 

dislocations in cell walls is r^ « L/K , and Gittus (1976a) gives the following 
expression for the steady-state creep rate 

e « - m - y / (dr /de) 

This is equivalent to (56) if one notes that for each length L of dislocation 
entering the cell walls after sweeping the cell cross section, annihilation 
required to keep a constant density in tite wall takes place by climb over a 
distance 

d • (5r « r ^ /I r.. L/k'"' (61) 

w w 

In addition, climb in the cell walls may take place under conditions of jog under- 
saturation in vacancies (e.g., Friedel, 1964; Hirth and Lothe, 1968); In that case, 
(58) must be replaced by (e.g., Poirier, 1976) 

2^\ b c j ^ « 2 ti b exp [Ej/RT] (62) 
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wh«re Cj 1 b the Jog concentration along the dislccatlon line and is the 
Jog formation energy* the Maxwell time is found to be 

T - C ^IMJL (63) 

where C ~ 3/K^. Cell wall recovery is therefore a more efficient mechanism, 

CW 

and is a high temperature mechanism since It possesses a higher effective 
activation energy* 

In the cell wall model of Glttus (1976a) the force causing climb Is 
associated with the self stress of the boundary network. Climb takes place 
so as to shorten the total dislocation length of the boundary and thus reduce 
the energy of the boundary* In this case, must be replaced in (57) by 
pb/2r^* This mechanism is the more efficient one when 

b/2r^ > o,j,/p (64) 

2 

This constraint Is satisfied for K > 2K^, In that case, the estimate (63) for 

s 

the Maxwell time remains valid, with a new value of the constant 


cw 


2K ^/k\ ^ 
s m 


(65) 


The expression (^7) for the climb velocity must be modified somewhat when the 
contribution of core diffusion is taken into account and when the wall dislocations 
are straight, with few geometrical Jogs* In that case, the appropriate 
expression is (e.g., Hirth and Lothe, 1968, p* 528) 


V . 4ir fliv ■ AE*/2)/RT] 

b (^/b) 


( 66 ) 


Here, AE* is the difference between the activation energies for bulk diffusion 
and core diffusion, and z is the mean free path of vacancies along the dislocation 
line 

z ~ b /T exp (AE*/2RTJ (67) 

- * , 

Expression (66) holds when Jog spacing is larger than z, that is when > AE /2; 
in other words, it assumes undersaturation. It is easy to see that use 
of (66) yields again a Maxwell time of the form (63), although the activation 
energy is somewhat smaller. The constant is not very different from (65); 

All the mechanisms considered so far (eqs. 55,59 ,63 ) lead to Maxwell times 

-2 .3 

proportional to ; in other words the creep rate, e, varies as in each 

case, and the main sources of variation among these models lie in the different 
constants and activation tiergies. A significantly different behavior obtains 
when the Gittus cell wall model is adopted and climb is controlled by double jog 
nucleation. Indeed, In this model, the dislocation segments in the cell wall 
networks have a length r^ which may be shorter than thermal jog separation. 

In that case the climb velocity (66) must be multiplied by c,/b (Hlrth and 

W j 

2 

Lothe, 1968, p. 528), and dependence of the creep rate results (Minster 

and ,'nderson, 1980). Note that this circumstance implies a high activation 

energy, and that t is only proportional to One might therefore expect 

this mechanism to dominate at high temperatures and low stresses, but this in 

turn would imply high thermal jog densities and longer dislocation segments, 

which might lead to a contradiction. Thus the existence of a domain in (T,o^) 

space where this mechanism dominates depends on whether the inequality 

r < b c.”^ can be satisfied without violating micros tructural scaling laws, 
w j 
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'litis can be written 


> 


u -E./RT 


0.r 


( 69 ) 


which places a lower bound on the stress and/or an upper bound on the temperature 
if the other parameters are known. In that case we have 


^sq " "^cw K exp [E./RT] 
® J 

c » — c 

sq cw 

s 


(70) 


For applications to the earth, we roust correct the Maxwell time expression 
given above for the effect of confining pressure. Tills correction term is 
unfortunately quite uncertain. We shall assume that, for practical purposes, 
the only effect is to increase the effective activation energy by PV , where 
P is the pressure and V is an activation volume. In the absence of a simple 

ft 

and reliable theory for V when core diffusion takes place, it seems unjustified 
to attempt a more sophisticated treatment. 


IV. APPLICATIONS AND DISCUSSION ; 

In view of the models presented above, attenuation and creep properties 
may be closely connected via the dislocation raicrostructure. We shall use the 
hypothesis that a single microstructural model explains both aspects of the 
rheology to place constraints on the model parameters. 

1) Constraints on creep models 

Tlte most stringent and useful constraints we can place on the various creep 
models outlined above is that they should satisfy laboratory data on olivine. 
These data have been summarized by Kohlstedt and Goetze (1974), Durham et al. 
(1977) and Goetze (1978). The range of stresses for which experimental data 
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arc available Ih barHt but the claas of moilols whicli we conalder here 

la probfibly only a|iptlcablc for Htrcancs smaller than 2 kbars (Goetsse, 1978; 
Kohlstodt, 1979). Kurtbornwrc, sliioc we are Interested In large regions of 
the mantle where the tectonic stress la probably only a few bars, we must assume 
that extrapolation of laboratory trends to U»w stresses is post lb le, using the 
scaling laws described earlier. 

Based on ilie general aspect ot the microstructure, and on some observa- 
tions on dislocation glide, Durham lU a L ( 1977), and (!oet?.e (1978) suggest 
that Takeuchi and Arg»>n's (19/(jh) creep model (e.g., e«|.!>9 ) may be appropriate 
lor laboratory detormed olivine, in the range ~ 100-2000 bars. The comparison 
is qualitative, however, and no numerical comparison was given. We have already 
shown lAuderson and Minster, (l9M0b) , that a straightforward application of 
ti’lttus' (|9/(>a) creep model could explain observations quantitatively. In this 
paper, we c»n>.slder a more general class of models and comv>are them systematically 
against various observational constraints, starting with direct estimates of 
some phvslc.)! parameters. 

When multi-species coupled diffusion takes place, 
self diffusion is largely controlled by the bulk diffusion of the slowest moving 
species (e.g,, Blocker and Ashbs% |97l). In silicates, this has long been 
thouglit to be oxvgen (e.g., Ashby and Verrall, 1978; Sammis et, »!.• » 1977). A 
recent est imate of the dlffusivity of oxygen in forsterite is (Reddy ^ iU,* » 1980) 

D (01 • 1,‘> It)”' exp I- a9000/RTl cm"/.sec 

Keci'utlv, bowevt'r, a lower dlffuslvitv has been found for silicon (.laoul 
et ad. . 19 79; PounK'llec et al,, 1980) 

D (St) « l.^ lO'*’ exp 1-90000/RTl cm“/sec 
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In this paper we shall consider either alternative to be possible and present 
results for both cases. In Mg2S10^, the appropriate atomic volume Si appearing 
in (55) is — tile molecular volume (7.2 cm /mole) if oxygen is the controlling 

3 

species, and the total molecular volume (28.8 cm /mole) if silicon is the 
controlling species (Stocker and Ashby, 1973). 

Undoubtedly a similar distinction should be made about the vacancy volume 
in the dislocation creep models. Such differences are rather trivial, 
however, compared with uncertainties in the other parameters. We shall make 
the simple choice 

E * » 90 kcal/mole 

s 

/ -3 2 

3.5 10 cm /sec (oxygen) 

D « , „ 

1.5 10 cm /sec (silicon) 

n ^ b^ 

V 

b ^ 7 10 cm (Stocker and Ashby, 1973) 

b) Creep activation energy ; The various models of dislocation creep which 
we have at our disposal entail activation energies which result from combinations 
of poorly known iquantlties, such as jog formation energy, etc.*. . On the other 
hand, the creep activation energy for olivine at high temperature has been 
estimated from laboratory observations. Recent interpretations of experimental 
data converge to 125 ± 3 kcal/mole (Kohls tedt and Goetze, 1974; Durham and 
Goetze, 1977; Ashby and Verrall, 1978; Goetze, 1978). 

Our approach is to assume that the creep activation energy is 
*= 125 kcal/mole and to consider tlie variation of unknown quantities subject 
to this constraint. Thus for cell^wall climb under applied stress (eq. 63) 






(71) 




we liave 


1 


K + K, ■ 125 
8 J 


; tij • 35 kcal/molo 
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(72) 


For coU-wal l cllml) under ael l-streHH , and core diffusion (et|. 65). 


Kj - AK /2 ■ 35 keal/inole 


(73) 


Here AK is the difference between the activation energy for bulk diffusion and 

^ ^ 

that for core cliliusion. For metalsi AE is 1/3 E to 1/2 E in many cases 

s s 

(e.g., Hirth and hotlie, I9b8) , On the other hand, we have argued (Anderson and 
Minster, (198(h) that AK may be close to zero for olivine since activation energies 
for lu)t gressing experiments (Schwenn and Ooetze, 1978) and dislocation loop 
mobility (Ooetze and Kolilstedt, 197'^) do not yield activation energies much lower 

it it 

than K^ and K,, respectively. 

A it 

For our purposes, we shall consider the entire range 0 :£ AB s 1/2 Eg . 

In the ci’.se of the o^ creep law (eq. 70), the constraint (73) must be rewritten as 


2K, - AK /2 “ 35 kcal/mole 


(74) 


Kquations (73) and (74) lead to the following respective ranges for the Jog 
formation energy K ; 


35 s. Kj < 58 kcal/mole 
17 V Kj V 29 kcal/mole 


(75) 

(76) 


TIu’se ranges indicate much higher values than for metals (e.g., Krledel, 1964), 
an iiccoptalvle conclusion in view of the greater cell sizes, and the higher Peierls 
energies Involved in silicates. These values may also be compared with tlie 
estimate of 26 keal/mole for the kink energy suggested hy Stocker and Ashby (1973); 
Jog energies should indeed be greater than kink energies (Hlrtb and Lothe, 1968). 
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llio rmigeB (75) and (76) Involve much Bmuller values than the geometrical 

estimates of order O.l |»b ~ 100 to 300 kcal/mole (Friedel, 1964). Geometrical 

estimates are admittedly very uncertain (Hirth and Lothe, 1968), even In the case 

of the simplest metallic crystal structures. In the absence of a suitable theory 

for silicates, we hypothesize that the discrepancy between E* and E* Is an 

c s 

empirical measure of the Jog formation energy, subject to the constraint equations 
(73) or (74) if core diffusion takes place. 

c) Microstructural constants. The scaling constants K_ and K_. can be 

' ' ‘ -■•■■■ .I..-.- ii.i rr- i.. I, g ||| 

estimated directly from the observations of Durham et al. (1977). Again these 

estimates are subject to fairly large uncertainties, and there is no guarantee 

that they can be used for extrapolation to mantle conditions. We take ■ 15 

and K " 10; these values lie in the lower range of observations, and we have 
m 

already argued that they may be appropriate. The constant K will be treated as 
a model parameter in our discussion. Domains of validity of the various creep 
models can be summarized as follows: 

- For cell-wall climb under self stress instead of applied stress (eq.64 ) 

K" •> 2K (77) 

s 

- For 0* creep law to be permissible (eq. 69) 

K" •> exp (- Ej/UTl (78) 

- An additional constraint must be placed on K for cell-wall recovery 

models with climb under Jog undersaturation in vacancies. It places an upper 

r 

bound on dislocation separation in Che walls: — * n < n , where n is a 

b w c c 

critical number of lattice spacings of the order of exp (Ej/RT], which is large 
In the present applications. It is commonly assumed that saturation occurs if 
the separation is of the order of 10 b — e.g. , in subgrains — (Friedel, 1964). 
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We shall assume, rather arbitrarily, that is at least a factor of 10 smaller 
(n^ < 1000), although this conrept is not very clear, If the transmission electron 
micrographs presented by Durham ^ al. (1977) are assumed to be typical, then n^ 
appears to be of the order of a few hundreds for applied stresses of a few 
hundred bars. 'Hie constraint can be written 

> K / ^\ (79) 

If it fails to be satisfied, then saturation may take place and the climb 
velocity is of the form (-’7); Jog formation is not a controlling process. On 
the other hand, as we mentioned earlier, for straight cell wall dislocations, 
with few geometrical Jogs, thermal Jog pair nucleation Is Important and the 
saturation condition 1s simply (e.g., Hirth and Lothe, 1968) 

Ej AE*/2 (80) 

Tills is always satisried iindor ooiistralnt (73), and imdor constraint (74) it 
implies AE ^ 70 kcal/mole. 

d) M ax w ell tim es! As shown by Ooetze (1978) « creep data on olivine at 
stresses greater titan about 2 kbars are well Fitted by a Dorn creep law. Below 
2 kbars a power law creep appears to give a reasonable interpretation. Figure 4 
shows the creep data of olivine collected for a variety of published sources by 
Kolilstedr and lJoet-/.e (1974), and supplemented with the creep observations of 
Durham and Ooetze (1977) on olivine single crystals with odd orientations (their 
figure 4), All data have been reduced to a common temperature of 1673“K 

using an jictlvation energy of E^^ « 125 kcal/mole. 

The observations of creep rate below 2 kbars can support eitlier a 
c " o' or a e « Interpretation, and Che scatter does not allow us to restrict 
tlu' number of acceptable models on Mtat basis alone. The bulk of the data, and 
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detailed obnervatiuna on individual samples by Durham and Coetze tend to 
support the higher exponent, but the lower stress observations seem to indicate 
a less sensitive stress dependence of the creep rate. As indicated on the 
figure, this can result in several orders of magnitude variations in the extrapo- 
lated creep rate at 10 bars, in addition to the intrinsic scatter. 

Each data point in figure 4 was converted into a preexponential Maxwell 

time using a rigidity of 650 kbars. The resulting stress dependence is 

-2 -1 

shown on figure 5, together with the o and a extrapolations to low stress, 

corresponding to the alternate interpretations of figure 4. Rough ranges for 

T at 10 bars are: 
o 

- For 0^ creep laws t (10 bars) ~ 10 to 5 10 seconds 

-13 -11 

- For 0 ^ creep (10 bars) 2 10 to 10 seconds 

We use these ranges as observational constraints to compare the various models 
discussed in the previous section. 

2) Comparison of creep models 

A convenient procedure to compare the varicus models is to compute the 
predicted Maxwell time preexponential factor calculated at a fixed temperature 
(1673“K) and stress (10 bars). 

a) Dlffiistonal creep (e<i.^^ ) : 

According to the preceding discussion, we do not have any free parameter 
left in this model, if oxygen diffusion is rate controlling, we find 

(1673,10) “3.7 10~^ sec, and if silicon diffusion is rate controlling, we 
find T^“ 2. I sec. Tliese values are many orders of magnitudes too large and 
we can eliminate this possibility. In addition, the predicted activation energy 
is 90 kcal/mole, which is much smaller than the observed value. 

b) Intracell recovery creep (eg. 59 ) ; 


In that case, we find 
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~ 6.2 10 ^ Hec for oxygen diffusion 

T 1.4 10 ^ sec for silicon diffusion 
o 

These eskimaCes arc still much too large to satisfy laboratory observations, 

and the predicted activation energy of 90 is too small. We nay therefore 

eliminate this model us well. 

c) Pile-up model (eg. 60); 

If oxygen controls self-diffusion, constraint (81) requires 130 < K_ < 470 
in that case* This range becomes 1700 < K < 6200 if silicon is rate controlling* 
These values are much larger than the observed range. Furthermore, the predicted 
activation energy is 90 kcal/molc. We must therefore turn to a more flexible 
model . 

d) Cell-wall recovery under applied stress (eg. 03 ) ; 

We now have the parameterB K and Cj at our disposal. As pointed out 

earlier, a Jog formation energy of 35 kcal/mole would reconcile the activation 

energy with the observed value. If oxygen is the controlling species for bulk 

diffusion, then at 10 bars and 1673"K, ~ 1.84 10 ^/K^; when combined with 

the constraints (81), this results in 14 < K < 37. If silicon Is the controlling 

species, we get 96 < K < 256. Thus, if silicon controls self diffusion, the data 

require a very pronounced concentration of tlie dislocation population in the cell 

walls, p being two to three orders of magnitudes greater than p . If oxygen 
w m 

controls self diffusion, then p^ is only required to be one order of magnitude 
or less greater than p . 

In either case, however, K is too large to satisfy constraint (77). In 
other words, the model Is self invalidating if the concept of cell-wall climb 
under self stress can be applied. 

e) Cell-wall r ecovery under self stress (eg. 64- 65) ; 

Testing this modeling hypothesis proceeds in the same fashion as above. We 
must postulate a Jog formation energy of Ej = 35 keal/mole to reconcile the 
activation energies. If oxygen controls self diffusion, then at 10 bars and 
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1673*K, ~ 4.1A 10 The obBervationnl constraints (81) yield in that case 

9 « K 19. If silicon is the controlling speciest we get 35 a. K s: 68. 

Thus If oxygen is the controlling species, the data imply a rather low 

range for K. In fact, it would he legitimate to argue that the micros tructure 
would be so poorly defined that an intracell recovery model should really be 

considered. We have already seen, however, that this model does not appear 
to be successful. 

Tlie hypothesis that silicon-controlled self diffusion limits the recovery 
rate by climb in cell walls under self stress is a more attractive one. Tlic 
density of cell wall dislocations Is predicted to be one order of magnitude 
or more greater than the density of tntracell, mobile dislocatior s in 
accordance with the premises of this model. This model cannot be rejected by 
our testing procedure. 

f) Cell wall recovery controlled by double Ior nucleation and core 
diffusion (eqs. <>6-68) ; 

This model is a generalisation of the previous one. We have now the extra 

it it it 

parameter AE at our disposal. If we let AE vary from zero to 1/2 E , the 

s 

corresponding variation of K, as Implied by the constraints (81), is shown on 
figure 6. The results are shown under the alternate hypotheses that oxygen or 
silicon controls self diffusion. 

Once more, the results for oxygen tend to conflict with the premises of 
the mlcrostructural model, whereas the results for silicon yield a self consistent 
picture, except possibly at very low stresses of the order of one bar, for which 
dislocation spacing in cell walls is large even with sizable values of K. This 
model is thus an acceptable alternative to the previous one, especially if 
further experiments should indicate a low activation energy for core diffusion, 
and a concomitantly high jog formation energy. 
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'I1ic presently avnUnblo dntn, a1 though rather limited and Indirect, tend 
to point to low values of AE (Anderson and Minster, 1980b) ; in that case the 
approximate range is 40 < K < 80, pointing to a sharply defined microstructure. 

g) o| creep law (eq. 70): 

This law applies to a cell-wall model where cell wall dislocation segments 
are shorter than the thermal jog separation. This leads to the constraint 
(78), and interpretation of the observations In terras of this rheology leads to 
the second set of bounds on t^^( 1673,10) In (81). If oxygen diffusion controls 
the self diffusion rate, then this model is unacceptable and must be rejected, 
as can be seen on figure 7; It does not satisfy (78) for admissible core diffusion 
activation energies. If silicon is the controlling species, then once more a 
range of satisfactory values of K is obtained. The model is generally associated 
with quite low values of Hj , however, and Is also incompatible with our conjecture 
that AE is small. We must therefore cmiBlder it to be less satisfactory than 
the previous ones. 

3) Creep in the mantle 

The systematic comparison of available models which we just performed does 
not yield a unique acceptable model for low-stress^ high-temperature creep. 
However, two general models satisfy all of our criteria in a self consistent 
fashion. Both involve cell wall recovery under self stress, both assume that 
silicon controls self diffusion, and both can be accommodated by the same 
microstructure, with 40 ^ K 70; either model constitutes a generalization of 
the creep model originally developed by Gittus (1976a). Their only major 

difference is that core diffusion is ignored in one case, and accounted for in 

* 

the other. Note that both models do not quite coincide for AE » 0, due to 
various approximations made in estimating the climb velocity (e.g., Hirth and 
Lothe, 1968). 
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In order to assess the applicability of our «odel to untie rheology* ve 
have calculated the Maxwell tlus for a range of applied stress Croa 1 to 1000 
bars, as a function of temperature, and for preexponential factors compatible 
with laboratory observations. For simplicity we have chosen 7 10 sec 

at 10 bars, noting however that laboratory evidence can allow for values about 
five times larger or smaller. Hils in turn corresponds to K ~ 50, which means 
that the total length of dislocations in cell walls if 25 times greater than 
the total length of intracell dislocations. 

The results are displayed on figure 8. Also shown is a box outlining a 
plausible range of mantle environments for comparison only (temperatures from 
1000"C to 1800“C; viscosity from 10*^ to 10^^ poises). Thus the model can 
easily accommodate geophysical observations in a quantitative sense. This 
calculation ignores the effect of confining pressure, however, and this brings 
a new complication to the problem. 

In a multispecies diffusion situation, the total activation volume for 
til 

diffusion of the i species is (e.g., Stocker and Ashby, 1973) 


* 


V 


i 


a Va + 8 Vg + Y Vq + 6 Vp + . . . 
a-f8+Y+<5+. . . 


V 

m 


(82) 


where a, 8 ... are stoichiometric coefficients; V^, Vg ... are atomic volumes 

of individual species; is the activation volume for motion of the i*^ species, 

which we assume to be close to an atomic volume. For Ng^SiO^, this leads to 
3 3 

11 cm /mole for oxygon, and about 4 cm /mole for silicon. The preferred ranges 

3 

of values of Sammis e^ (1977) for the mantle decrease from 11 to 6 cm /mole 

3 

for the upper mantle, and from 6 to 3 cm /mole for the lower mantle. Tliey point 
out that this bounds the total variation of viscosity across the transition 
region (300-1000 km) to be less than two orders of magnitude. This argument was 
used in reverse by Anderson and Minster (1980a) , based on the modeling work of 
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3 

Peltiar (1980), to argue for an effective activation voluM between 4 and 9 c« / 
mole. 

If we ignore for the moment the preaeure dependence of the activation 
volume, the estimatee for silicon and oxygen given above have an intereatlng 
and potentially Important implication. Should silicon be the controlling 
species, then the pressure correction to the effective creep activation energy 
is about O.l P kcal/raole, where P is in kilobars; this correction is about 
0.27 P kcal/mole if oxygen is the controlling species. If we now fix the micro- 
structural model, the Maxwell time prcexponentlal factors under both hypotheses 
are practically in the ratio of the diffuslvlties. At zero pressure this ratio 
is essentially independent of temperature because the activation energies are 
almost identical. According to our systematic evaluation of creep models, 
oxygen diffusion is then too rapid to satisfy the laboratory observations, 
and silicon diffusion appears to be rate controlling. The effect of pressure 
may reverse these roles, as shown on figure 9. We calculated Maxwell times, 
for a fixed load of 10 bars, « 7 for silicon diffusion), under the 

alternate hypotheses that silicon or oxygen are the controlling species, and 
for confining pressures up to 300 kbars. Because of the larger activation 
volume assumed for oxygen diffusion, we find that oxygen becomes the slower 
species at high pressure for a fixed temperature. As an indication of mantle 
conditions, we considered a geotherm with a linear temperature gradient of 
0. j "C/km, and a temperature of HOO^K at 400 km depth, and plotted the corres- 
ponding trajectory on figure 9. The switcli over from silicon controlled diffusion 

takes place <nt a pressure between 100 and 200 kbar. (300 to 600 km depth). 

Activation volumes are expected to decrease with increasing pressure based 

on an elastic continuum model (e.g., O'Connell, 1977} Sanunis et al . , 1977). 
in applying (82), tt is clear that the activation volume for silicon diffusion 
is controlled in turn by the much larger atomic volume of oxygen. As a result. 


. , 
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it can be expected that both cHtlmatCK of the activation volume might decreane 
at depth in about the same proportion. Thiu will have two consequences: 

1) The switch over from silicon controlled diffusion to oxygen controlled 
diffusion will happen at greater pressure and thus greater depth, ^nd 

2) The smaller activation volume in the lower mantle will tend to limit the 
variations of apparent viscosity with depth (Sammis et , 1977). 

These comparisons are very crude, and the results can vary considerably 

with model parameters. We have ignored chemical variations, changes in lattice 

structure with pressure, effects of water content on the activation energy... 

etc. Our microscopic models are probably not yet precise enough to support 

such sophistication. However, some general inferences can be made, which are 

robust with respect to the variations in interpretive models. For example, 

15 

at mantle temperatures, very low effective viscosities (as low as 10 poises) 
would be observed under loads of 1 ki lobar. Thus the model precludes the 
persistence of such loads for times longer thin the corresponding Maxwell time 

3 

(<^ 10 sec). Loads of the order of 100 bars might explain the low apparent 
viscosities Inferred by Nur and Mavko (1974) on the basiu of post seismic rebound 

3 

data. On the other hand, persistent loads of the order of 10 bars can be 
sustained by a cold lithosphere over geological times (e.g., Anderson and Minster, 
1980a;Minster and Anderson, I960). Reinterpretation of isostatic rebound data 
in terms of a phenomenological power-law flow compatible with our models was 
performed by Crough (1977). His analysis is a particularly striking illustration 
of the tradeoffs between temperature, confining pressure, and applied load of the 

kind shown in figures 8 and 9. His results, for an applied load of 10 bars, 

8 7 

yield Maxwell times ' 3.5 10 sec for the Canadian shield and 1.7 10 sec for 

Bonneville. As seen on figure 9, this is in good quantitative agreement with 

the predictions from the microscopic model. 
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4) C onBtraintK on otteiiuatloii 

For the purposes of dlscuBHlng attenuation, we may recast equation (9) 
in a form comparable to the expressions for the Maxwell time by introducing the 


critical length (e.g. , Figure 3); 

^2 kT cxpIEq / kT] 


/iL.f 

? (o.j.1 

^ (kJ 




(83) 


(84) 


where is the critical lengtii associated with the long term background 
("tectonic") stress o^. The most interesting implication of this formalism 
is, of course, that in order to insure consistency with the argument leading 
to (31), we must now postulate a very narrow distribution of lengths, clustered 
about the value i == 1.66 1 ^. This has several Important consequences. 

The first and most straightforward is that the relaxation strength is 
given by (16), and if the scaling law (51) holds, A is of the order of 10% 
and is Independent of stress. 

The second consequence is t*'at the absorption band Itself depends rather 
strongly on the background tectonic stress. It shifts to high frequency with 
increasing stress and to low frequency with decreasing stress. Comparison of 
(83) with the Maxv;ell time yields in fact 


D 


' I - bJ)/KTI 


(85) 


cw 


Thus, given microscopic models for creep and attenuation - in particular, 
given the dif fusivities and activation energies we have a formal connection 
between steady state creep, attenuation and tectonic stress, through the 
microstructure. This opens the possibility of formalizing the empirical 
correlations studied for example by Anderson (1966), Meissner and Vetter (1979) 
and Berckhemmer £t aJL, (1979) within the framework of a physical model. 
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'nu> ihlrU ronKtMiut'm'o of rirtK t1u' micros true tuns to be 

comp/it lb U> with the or«M'p imulel l« that a distribution of lengths la no 
longer (iv.iilahle. ’Huis eq. (12) cannot be at the source of a broad absorption 
band. To get a spread of proexponential factors, we must turn to either a 
distribution of or of X. 'ilu' alternative, of course, la to rely on a 
distribution of activation energies to broaden the absorption band, in that 
case, the width and intensity of the band depends mainly on tomporature, but 
its position along the frequency or temperature axis depends on tiH tectonic 
stress 0 .J,, 

Numerous authors have considered diffusion contrelled dislocation 
damping as a possible mechanism for anelastlc attenuation at low frequencies. 

A fairly recent treatim*nt Is that of Simpson and Sosln (1972). Some of 
the imiln al teinat Ives are listed below. 

'*) I Vleri s barrier .ii»di^l, j>.f . tiu^^ Ht>rdonl peak 

A I’clerls harrier model of the HordonI peak gives a relaxation time 
of the correct form. If we use the simplified theory of Hirth and Lothe (1968, 
p, 100) , tia' activation energy Is 2l''j^» where Is tin? kink formation energy, 
estimated to he -v* 2h kcal/mole by Stocker and Ashby (1973). The preexponential 
dlMusIvltv Is V h** where v la an atomic frecpiency of the order of 

10^-^ sec“^, so that 'V* 10 - to I0~‘' cue /sec. The scale length \ is found 
to he \ V (p/o,j,)' ( 0 |/u), where Oj denotes the Internal stress. For 

Intragraln processes, with fi'w dislocations, Oj Is expected to he controlled 
malnlv hv the line tension (Oittns, 197ha), and thus to be comparable to 
o,j,. \ Is therefore of order ^p/o,|, and Is of the order of 10”^^ to 10“^^’ 

seconds, and Is practically Independent of o,j,. We heive therefore a high 
frequency, low temperature mechanism, which does not operate In the seismic 
band and In the mantle. A more detailed review of the HordonI peak is 
given bv Nowlek and llerrv (1 9/:'), 


A-m-39 


I’ ) oou trol lod by .k Ui f ftiai on 

(Jlldo controlled by kink diffusion (Uirtli and Lotlie, 1968, p 697) miglit 
be appropriate In view of tlio blgb I’elerl -nergy of olivine (e.g. Gueguen, 
1979). For an Initially straight segment of lengtii I, bowing out by kink 
pair generation and diffusion we have 


A 




2 b (t 

t -f b exp IE,^/RT) 


( 86 ) 


1)^ -V vb^ exp 1-15^/RT| 


(87) 


Two clreumstanceH iiuiv arise, depending on wlileh term dominates the 
denominator of (Ho). 

Let us assunK>, for tin* sake of arguiwnu, that « 26kcal/roole (Stocker 
and Ashby, 1973), and tiuit o.j, ■ 10 bars. Then, for temperatures lower than 
the transition temperature T^ ~ 1120”K we have A ‘V 2 1 exp (-Ej^/RTJ. Tlte 
effective activation energy is then E^^ 'v 2 Ej^ 52 kcal/mole, and we have 
'V. lO” * sec. Above , we have A 2b, E^^ a. 26 kcal/mole, and 10 
In either case this yieldu an absorption bank at much higher frequency than the 
seismic band. It must bo rocogalxod, i»ow/ver, that the kink formation energy 
Is very uncertain, and that the transition temperature is very sensitive to 
both K|^ and o.j,, as shown in Table 1. Ibus, for * 10 bars and Ej^ “ 35 kcal/mole, 

I ‘ ^ 

we find 1 '• 5 I0~ ■ sec, E,. 70 kcal/mole, and the absorption band approaches 

o'! 

1 Hr. near I500"K and longer periods at lower temi>eratures. Tills mechanism is 
therefore a potential oaiulldate, particularly In regions where tlve stress is 
fairlv hlgj» and the temperature fairly low, such as in the lithosphere. 

Glide of J o gged s crew dislocation s 

Yet another alternative is provided by the glide of jogged screw dislocations, 
where tiu* Jogs can only move by the production or absorption of point defects 
(Hirth and l.othe, l9nH, p 335). In that case the diffuslvlty Is that for 
self diffusion H which in turn Is related to the diffusion of oxygen or 
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silicon according to our discusHlon of creep. For small stresaes we have 

X - Aitb/C| where Cj is the Jog concentration. Two cases arise; either 

the Jogs are geometrical, and X is of order b to within one or two orders of 
* * 

magnitude, F,,. is then •'< « 90 kcal/mole; or the dislocation segments lie in 

'< ” 

Peierls valleys and we must consider thermal Jogs (Gittus, 1976a), In the 

it it 

latter case X is 4vb and the activation energy is "V E + E, ^ 125 kcal/mole 

0 s j 

in accordance with the creep model. In either case, it is easy to verify that 

the retardation time is many orders of magnitude longer than the seismic 

periods. For example, under an applied stress of 10 bars, we get t “ 10^® sec 

at 1600®K if silicon controls self diffusion and thermal Jogs are the dragging 

defects. On that basis, Anderson and Minster (1980b) rejected this mechanism 

as a candidate for seismic attenuation. For the case of geometrical jogs, a 

lower bound for t is obtained by taking X * Awb and using the dlffusivlty of 

oxygen. Then v 'v 5 10® sec at lb00®K and 10 burs. Using silicon dlffusivlty 

and/or, taking the effect of pressure Into account can easily Increase this 

value by several orders of magnitude. We may conclude that this mechanism might 

be a candidate for the attenuation of the Chandler wobble, although it is 

not clear that the anelastlc theory of section 11 still applies when the 

retardation times become comparable to the Maxwell time for steady state creep. 
Dragging of poin t defects 

Dragging of an atimisphere of point defects bound to the dislocation is 
a more versatile and successful model. Original applications to internal 
friction are due to Weertman (1955), and Frledel £t |CL. (1955). A specific 
model was proposed by Schoeck (1963) to explain the cold work peak in b.c.c. 
metals. Further improventents and modifications were subsequently proposed by 
Barrand and Leak (1964), Cibala (1967), Ino and Sugeno (1967), and more 
recently by Miner e_t Jij . (1976). Covalent crystals were considered by 
Southgate and Mendelson (1965). Dragging of a Cottrell atmosphere Is an 
essential Ingredient of Takeuohi and Argon’s (1976) model for steady state 
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creep. Binding of these defects to the dislocation line is usually assumed 
to be elastic, although electrical interaction probably is important in 
ionic crystals at hlRh temperature (Eshelby et 1958; Brown, 1961; 

Menezes and Nix, 1974 a,b). 

This class of models has been demonstrably successful in the case of 
simple crystal structures. We propose to extend it to mantle materials. 

The point defects which contribute to drag may include self Interstitials, 
Impurity interstitials, or substitutional atoms with a slightly different 
ionic radius if the binding mechanism to dislocations is elastic in nature. 

Defects associated with a local charge Imbalance (Smyth and Stocker, 
1975; Stocker and Smyth, 1978; Stocker 1978) might be Involved if electrical 
interaction is important. Let E^ be the binding energy of the point 
defect to the dislocation and c the bulk (lattice) concentration of defects 
away from dislocations. Then a very simple model of equilibrium between 
bound (Mid unbound defects is obtained under the assumption that Ferml-Dirac 
statistics are applicable (Hlrth and Lothe, 1968). The high temperature 
approximation yields a mean distance between defects along the core 
(e.R. Miner iHc aj_. 1976) 


X 'V -- exp (-Ejj/kT) 


( 88 ) 


where binding energy is of tiu' order of a few kcal/mole for many metals 
(e.g. Nowick and Berry, 19711). But the analysis leading to (88) assumes 
implicitly a high level of purity of the material. This might be appropriate 
for intragrain environments if most impurities eventually migrate to grain 
boundaries as we have argued before (Anderson and Minster, I980b). On the 
other hand, if the binding energy is only 1 e.v., then the exponential is 
of order 10”^ at mantle temperatures, and is thus of the same order or 
smaller than Cj . The higli temperature approximation (88) is no longer 
valid in that case and we must consider the case where a very large fraction 
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of core sites are occupied by a draRglng point defect. The scale length X 
Is then comparable to b, and this limit yields an upperbound for t. 

Under these circumstances both 1)^^^ and E* are effective quantities 
for the diffusion of defects bound to the dlslocatlo.i core. At sufficiently 
high temperatures, where (86) holds, E* would be augmented by the binding 
energy. In the absence of a good theory we shall Ignore such complications 
for the present. Dlffuslvity data which are potentially relevant Include 
the cationic diffusion parameters In oxides, silicates and alumlnates compiled 
by Ahrens and Schubert (1975), and a larger data set for metal oxydes compiled 
by Kofstad (1972). If we restrict our attention to measurements which involve 
a broad range of high temperatures, typical preexponential dlffuslvities are 
in the range 10“^ - 10“^cm^/sec, although order of magnitude fluctuations are 
not uncommm. For a tectonic stress of 10 bars, this Implies a preexponential 
retardation time -v 10"®to lO"® seconds. If X > b this estimate should be 
lower. If we choose « 10"® seconds, then In order to bring the absorption 
band Into the seismic range (1 sec to 10 sec) at 1600®K, we require 
45 'vEq^ 65 kcal/mole. An increase of X to 10^ b changes this requirement to 
65 ^ Ep ^ 90 kcal/mole. These values must be lowered if we choose a lower 
reference temperature, or if a lower tectonic stress Is assumed. The 
permissible range encompasses the scatter of activation energies listed by 
Ahrens and Schubert (1975). It Includes the range of estimates of 47-58 kcal/mole 
for interdlf fusion of Mg^ and Fe^ In olivine (Misener, 1974). 

Unfortunately, there Is little experimental evidence available to further 

constrain the class of acceptable models. Jackson (1969) found x "v 4 10"^^ sec 

o 

* 

and E 57 kcal/mole for an absorption peak near 1 Hz In fine grained olivine. 

He Interpreted It as a grain boundary peak. However, we are Interested In 
larger grain sizes, and higher temperatures, and therefore In the high 
temperature background attenuation seen in Jackson’s data, for which estimation 
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of the relaxation parameters is difficult (Anderson and Minster, 1980)# 

Berckhemmer et al_. (1979) used linearity to interpret transient creep 
superposed on a steady state regime, in natural peridot ites, in terms of 
attenuation Q(h>). However, the transient stress and steady state load used 
were both in the 1 - 10^ bars range; it is not clear whether the superposition 
of transient and steady state creep can be handled as we have in our model, 
since the time scales and stress ranges tend not to be so clearly separated. 
The authors report an apparent activation energy of 30 kcal/mole for Q(m), 
which they interpret as uE_ with « *'< 0.3 (e.g. Anderson and Minster, 1980a). 

If instead we choose to interpret this number as Eq, it falls on the low 
activation energy side of the bulk of cationic diffusion data. This however 
would be inconsistent with the observed frequency dependence, i.e. (hit)” 
varlncion of Q In these experiments. 

Although we do not have enough experimental constraints at our disposal 
to propose a better specified model and then test it against geophysical data, 
we can still use seismologlcal observations to limit the range of permissible 
models . 

5) Absorption band in the mantle 

The model affords a large enough number of parameters (including tectonic 

stress o.j.)» which may all vary with position, that there can be no doubt 

that all available geophysical observations could be explained by it. For 

purposes of illustration we •»i.all only consider here the tradeoffs between 

temperature, pressure, and tectonic stress and consider other parameters to 

be fixed. In particular, we assume a low purity for mantle material, 

so that A b and t 10“® sec at o_, = 10 bars and T ''' 1500“K. 

o 1 

At any given point in the mantle, the width of the band is then only 
controlled by the spectrum of activation energies. Since the strength A is 
fixed to 10% in this model, the intensity of the absorption band (in the 
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case a » 0, eq. 37) is given by 

A n 2 «E 

\ irA RT 


(89) 


where 5E* is the range of activation energies* Several authors have argued 
that a low value of is required by seismic observations , at least In the 

regions of the mantle where the absorption band coincides with the seismic band. 

A value 'v> 25 is in good agreement with the seismic data (Solomon, 1972; 
Anderson and Hart, 1978). Sipkin and Jordan (1979) report a Q (ScS) 170 for 
periods of 10 sec to 1000 sec. If ScS should spend only one fourth of its 
travel time In the low Q regions of the earth, then ^ 40. Equation (89) 
allows us to place a limit on the range of activation energies: at mantle 
temperatures 5E* ^ Q„,/2 (kcal/mole). A range 45 ^ E* ^ 60 kcal/mole is 

therefore consistent with this constraint and with the laboratory data on 

cationic diffusion. Note that this yields an absorption band which is only 
about two decades wide at mantle temperatures. One consequence is that the 
value of a in (11) is rather immaterial, so that we will choose a ■■ 0. 

For a temperature of 1500“K the range of retardation times is 3-500 seconds, 
which lies within the seismic band. An Increase of temperature shifts the 
absorption band to higher frequencies; an Increase in pressure has the opposite 

effect and so does a decrease in tectonic stress. Based on their observation 

that Q(ScS) increases with frequency In the range 0.1-1 Hz, Sipkin and Jordan 
(1979) suggest an effective cutoff for the mantle absorption band averaged 
by ScS near one second. Comparable conclusions have been reached by other 
investigators (e.g. Kanamori and Anderson, 1977; Minster, 1978b; Lundquist, 

1979). In terms of our model, a high value for Q(ScS) in the 1 Hz band means 
that this frequency does not lie within the frequency band for an extensive 
range of depths. The much lower values of Q(ScS) reported by Sipkin and 
Jordan in the 10-1000 sec band suggest in fact that the absorption band 
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remains at periods longer than I second in spite of increasing temperature 
with depth. 

This may be due to either n decreaKc in tectonic stress at depth, or 
to the effect of pressure if the activation volume is large enough. For 
example a change of one order of magnitude In shifts the band by two 
decades In frequency. With an activation energy of 50 kcal/mole and at 
mantle temperatures, a similar shift is achieved by a temperature change of 
about 400*C. Thus It may be expected that the stress-temperature tradeoff 
is most important in regions with a high geothermal gradient and a simultaneous 
decrease of stress with depth, such as the upper mantle (e.g. Schubert et 
1976). On the ocher hand, if dynamical transition zones (boundary layers) 
mark the bottom of convection cells In the mantle (Jeanloz and Richter, 1979) 
Chen both temperature and stress Increases would tend to combine and cause 
a rapid shift of the absorption band Coward higher frequencies. This might 
provide a physical interpretation for a low Q zone of the sort hypothesized 
by Anderson and Hart (1978) at the base of the mantle. 

Discussion of the effect of pressure is predicated on the adoption of 

it 

an activation volume V^, for which we have no data. For purposes of discussion, 
let us consider again a geotherm with a linear 0.3"K/km gradient and a 
temperature of 1700®K at 400 km depth. For a pressure independent activation 
volume V*, and at constant tectonic stress the effect of pressure will 
overcome the effect of temperature If 



(90) 


or 


V* ^ 2.5 10-2 


(91) 


where V* is In cm^/irole and is in kcal/mole. Thus, with a spectrum of activation 


I 
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energies between 45 and 60 kcal/mole, and the geotherm described above, 

> the absorption band shifts to longer period with Increasing depth if 

it * 

Vq ^ 111 cm^/mole. In order to Illustrate this phenomenon we show in 
Figure 10 the absorption band and associated dispersion, for a fixed stress 
I <■ 10 bars, at depths of 100, 200..., SOOkm, for which the assumed 

it 

temperatures are shown on the figure. was fixed at 4 cm^/mole, a value 
large enough to illustrate the initial shift of the band to high frequency with 

I 

I rapidly Increasing temperature, followed by a progressive shift toward longer 

I periods, due to increased pressure. 

I More complicated scenarios can of course be envisaged , with a spectrum 

' of activation volumes, phase transitions, etc... Some of these refinements 

I have been considered by Lundqulst (1979, 1980). They go well beyond our 

I 

I present purposes. 

6) High temperature background attenuation (HTB) 

I 

We have suggested before (e.g. Anderson and Minster, 1980) that the 
dominant mechanism of attenuation in Che mantle corresponds Co the HTB in 
laboratory experiments. A description of its properties is found in Nowick 
and Berry (1972). An Important, unanswered, question is whether the HTB is 
anelastlc or viscoelastic in nature. The work of Friedel ^ (1955), and i 

more recently the analysis of Woirgard (1976) clearly indicate the possibility | 

of a broad, high temperature anelastlc peak. Figure 11 shows our model at zero j 

i 

pressure and at frequencies comparable to laboratory experiments. It points | 

Co a testable aspect of Che model, namely, Che development of a peak, at \ 

sufficiently low frequency and/or high temperature. Observation of such 
a peak, in well annealed material deformed under low stress to steady state 
would not only either provide support for our model or infirm it, but it would 
yield needed direct estimates of the model parameters. 

i 


I 
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V. CONCLUSIOMS 

The main conclusion of this work is that a dislocation micro- 
structure model is capable of simultaneously explaining steady state cteep 
and attenuation in the mantle. Our preferred model involves a sharply 
organized mlcrostructure, where most dislocations are found in cell walls. 
Steady state creep is controlled by recovery through climb la the cell walls. 
Good quantitative agreement with laboratory observations of creep in olivine 
is obtained if silicon controls self diffusion. However, because of the 
larger activation volume of oxygen, we expect that oxygen diffusion will 
become the controlling factor at sufficient depth. We hypothesize that 
the discrepancy between oxygen or silicon diffusion activation energies and 
the observed activation energy for creep is a measure of jog formation energy, 
which is then somewhat greater than 35 kcal/mole, depending on the efficiency 
of core diffusion. 

Attenuation of seismic waves takes place through bowing of the few 
intracell dislocations, controlled either by kink diffusion (particularly at 
high frequency) or by the dragging of an atmosphere of point defects bound to 
the dislocations. Due to the paucity of laboratory observations, the success 
of this model is judged on more qualitative and circumstantial evidence. 

The microstructure depends on the long term tectonic stress, which exerts 
a strong influence on both creep and attenuation properties. Within the frame- 
work of this model, stress, viscosity and Q are related, so that estimating 
one of t!»ese parameters provides constraints on the others. The scaling laws 
do not favor a wide spread of length scales, so that the relaxation strength 
is fixed in this model, and is of the order of 10%. A broad absorption 
band can be achieved with a spectrum of activation energies, but selsmological 
requirements ns to tlie intensity of attenuation in the band limit its width 
to 2 or 3 decades. The main parameters controlling the position of the low Q 
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bnnci are tempornturo, pro»Kurc, and tectonic stress. The absorption band 
coincides with the seismic band In the upper mantle, but shifts to longer 
periods at (treater depths due to the combined effects of increased pressure 
and possibly lowering the shear stress. It must be noted, however, that 
sharp shifts to shorter periods are expected In hypothetical boundary l.i\yers, 
due to larger temperature gradients and higher stresses. In the bulk of 
the mantle, however, seismic frequencies appear to be on the high frequency '> 

low temperature side of the absorption band. The absorption band can also be 

* 

expected to be broadened in the mantle because of the distributions In D^,E , 

I and M expected for a polymineral ic aggregate and the anisotropy of the 
grains. Each mineral has its own physical properties but, at least in the 
upper mantle, the total breadth of the band Is constrained by the minimum 
allowable damping. 

Our models adopt a very simple view of dislocations In complex crystal 
structures pertinent to silicates. We have ignored the complications associated 
with partial dislocations, and the differences between dislocations with different 
Burgers vectors. A more detailed discussion Is provided, for example, by 
Gueguen (1979) (see also the comment by Poirier at the end of Goetze*s (1978) 
paper). The dislocation model proposed here may be viewed as a possible 
alternative to grain boundary mechanisms described for example by O'Connell and 
Budlansky (1977). Its validity must be ultimately tested against further 
experimental data. 
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TABLE l ! Transit ton temperature T^. 


Om 

\ 

1 bar 

10 bars 

100 bars 

kcal/mola 

25 

540*K 

1080'K 

1350*K 

35 

1250"K 

1510*K 

1885*K 

45 

1620*K 

1940"K 

2420*K 
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FlRure Captions 


Sketch of a bowinR dislocation of length 1 strongly pinned 
at both ends. Y is the migration velocity, the radius of 
curvature. 

Three examples of absorption band and associated dispersion 

characteristics, for a x 0.25, 0, and 0.25 respectively. 

T « 1 sec, T„ ■“ 10** sec. Relaxation spectra have bera adjusted 
n M 

so as to yield identical high frequency behavior. 

Comparison of relaxation time t and multiplication time t^ 
as function of dislocation length normalized to critical length. 

Steady state creep rate observations in olivine compiled by 
Kohlstedt and Goetze (197A) (circles) supplemented by observations 
of Durham and Goetze (1977) (triangles), as a function of applied 
stress. Data for stresses greater than 2000 bars are shown 
with smaller symbols. 

Same data as figure 4, converted to Maxwell times, assuming a 
rigidity p “ 650 kbars. 

Structural f rnraeter K as a function of difference between bulk 
diffusion I core diffusion, for cell wall recovery model controlled 
by double Jog nucleation. Also indicated is the jog formation energy 
which yields a creep activation energy of 125 kcal/mole (eq. 73), 
and on tlie left, the dislocation separation in the cell walls, for various 
values of ilu» applied stress. Shaded areas indicate range of K values 
which bracket laboratory data shown on figure 5. Forbidden area for 
this model (eq. 77) is also indicated. 
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Same as figure 6, for creep lew. Jog formation energy 
calculated by equation (74). Forbidden area corresponds 
to eq. (78). 

Fig. 8 ; Temperature dependence of Maxwell time for several values of 
applied load, for a creep model which satisfies laboratory 
observation* vljn a o® interpretation. 

Fig. 9 t Effect of confitilng pressure on Maxwell times. At pressures 
greater than 200 kbars, oxygen diffusion becomes slower and 
controls self diffusion. Simple thermal model of the mantle 
is shown, and compared with Crough's (1977) estimates for 
isostatic rebound. 
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Hypothesized mantle absorption band for A » 10%, ■ 10 bars, 

45 £ Eq 5 60 kcal/mole, Vq *» 4 cm^/mole, calculated at pressures 
of 100 to 500 kbars, for which Indicated temperatures were assumed. 
Also shown are the dispersive properties associated with the 
absorption band. 


Fig. 11 ; Temperature dependence of absorption band of fig. 10, for three 


frequencies appropriate for laboratory, Indicating possible interpretation 
in term of high temperature background attenuation. 
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ABSTRACT 
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The basalt, eclogite, and harzburgite that are the differentiation 
products of the Earth appear to be trapped in the upper mantle above the 
670 km seismic discontinuity. It is proposed that the upper mantle 
transition region, 220 to 670 km, is composed of ecloglte which has 
been derived from primitive mantle by about 20% partial melting 
and that this is the source and sink of oceanic crust. The remainder 
of the upper mantle is garnet perldutlte which is the source of 
continental basalts and hotspot magmas. This region is enriched in 
incompatible elements by hydrous and CO 2 rich metasomatlc fluids which 
have depleted the underlying layers in the L.I.L. elements and L.R.E.E. 

The volatiles make this a low-velocity, high attenuation, low viscosity 
region. Tlse eclogite layer is internally heated and it controls the 
convection pattern in the upp'r mantle. Material can only escape from this 
layer by melting. The insulating effect of thick continental lithosphere 
leads to partial melting in both the perldotlte and ecloglte layers. Hotspots 
and ridges mark the former locations of continents. Most of the 
basaltic traction of the oceanic lithosphere returns to the ecloglte layer. 

Plate tectonics is intermittent. The continental the'.mal 
anomaly at a depth of 150-220 km triggers kimberlite and c.arbonatite 
activity, alkali and flood basalt volcanism, vertical tectonics and 
continental breakup. Hot spots remain active after the continents 
leave and build the oceanic islands. Mantle plumes rise from a depth 
of about 220 km. Mldocenn ridge basalts rise from the depleted layer 
below this depth. Material from this layer can also be displaced upwards 
by subducted oceanic lithosphere to form back-arc basins. 
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INTRODUCTION 

Although convection plays an important role In plate tectonics 
and heat transport In the Earth It has not succeeded In homogenizing 
the mantle. Magmas are still being produced from mantle reservoirs 

9 

which have remained separate for the order of 1 to 2 x 10 years 
(e.g., De Paolo, 1979; Sun and Hansen, 1975). Oceanic lithosphere is 
continuously returned to the mantle but the difference in element 
ration in the reservoirs, e.g., Rb/Sr, U/Pb, Th/Pb and Sm/Nd, persists. 

If the depth of earthquakes in subductlon zone^. can be used as a guide, 
oceanic lithosphere is presently being delivered to the region of the 
mantle betv^een about 220 km and 670 km. The isotopic data can be 
satisfied if this is also the source region for mid>oceau ridge basalts 
(MORB). This leaves the upper mantle or the lower mantle as the source 
region for continental flood basalts (CFB) , hotspot magmas and ocean 
island basalts (OIB). The upper mantle low-velocity zone (LVZ), or 
aathenosptiere, is the mote likely source region since temperatures there 
are closest to the melting point. 

37 83 

Ocean floor basalts have comparatively uniform and low Sr/ Sr 
^®^Pb/^^^Pb, and ^^^Nd/^^^Nd ratios whereas continental magmas 
and basalts from ocean islands not associated with island arcs have less 
uniform and higher ratios (De Paolo and Wasserburg, 1979). The latter 
magmas are also enriched in volatiles and the incompatible large-ion 
lithophile (LIL) elements (Wiite and Schilling, 1978; Frey et al., 1978). 
The study of Isotopes has introduced the time constraint that reservoirs 
with different element ratios — Rb/Sr, U/Pb, Th/Pb, and Sm/Nd — 
have existed for the order of 1 to 2 b.y. The source region for MORB 
has been providing uniform composition lavas for 16ng periods of time. 
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It must therefore he Immense In size and global in nature (Schilling^ 1975). 
The reservoir for continental and ocean Island magmas also appears to be 
relatively uniform and global although Its products are often mixed with 
varying amounts of MOEB before It is sampled. 

There are two competing petrological viewpoints regarding the 
nature of the source regions. The common view is that all basalts 
represent various degrees of partial melting of a garnet perldotite* 

The alternative position is that some basalts represent extensive melting 
of a deep ecloglte source. Both ecloglte and garnet perldotite 
inclusions are common in kimberlite pipes. The ecloglte inclusions, 
although not rare, represent only about 20% of the total. This suggests 
that ecloglte is either a less abundant component of the mantle or 
it occurs deeper than the garnet perldotite, as befits its higher density. 
Neither ’of the two types of fragments can represent primitive mantle 
(Allsop er. al., 1969). They must therefore be a result of a previous 
differentiation event. Hie ecloglte minerals are depleted in the 
trace elements which are enriched In perldotite nodules, the plume 
source region and the continental crust. It is therefore desirable 
to test the hypothesis that ecloglte, perldotite and continental crust' 
are the principle products of mantle differentiation and that xenoliths 
in kimberlites may be samples from the mantle source regions. If true, 
this would have considerable impact on our ability to model the composition 
and evolution of the mantle. 

CHEMICAL STRATIFICATION OF THE MANTLE 

The mantle Is also heterogeneous in Its seismic properties. It* has 
not been clear, however. If or how the seismological and geochemical 
heteroteneities are related. The largest lateral variations in seismic 
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velocities occur In the outer 200-250 km of the Earth and are related 
to such surface tectonic features as shields, trenches, rises, and 
volcanic belts. The mantle is also Inhomogeneous radially with the 
lithosphere, asthenosphere, and transition zone being the main sub- 
divisions of the upper mantle. 

A chemically layered upper mantle could provide distinct and 
isolated reservoirs and Is more suitable in many ways than models 
Involving Isolated blobs or regional Inhomogenelties (Hofmann et al., 
1978). It has been proposed that t'bs low-velocity zone Is the depleted 
reservoir and the source of mid-ocean ridge basalts (Schilling, 1973). 
Plume basalts, i.e., magmas from the L.I.L. enriched reservoir, have 
been attributed to deeper sources. If the LVZ Is enriched In volatiles, 
as proposed by Anderson and Sammis (1970) on geophysical grounds, then 
this explanation is untenable. Prey et al. (1978) have discussed other 
objections to this model. They argue that volatl].es should have 
enriched the upper layers of the mantle. 

On the basis of seismic velocities and seismicity patterns, 

Anderson (1979c) proposed that there were chemical discontinuities 
in Che mantle at 220 and 670 km. The former is the base of Che LVZ 
and near the maximum depth of earthquakes In continents] collision zones 
and regions of subductlon of young, <50 Ma, oceanic lithosphere. The 
latter is a sharp seismic discontinuity and is near the maximum depth 
of earthquakes. Only old oceanic lithosphere penetrates this deep. 

The seismic velocities between 220 and 670 km are consistent with 
ecloglte. 

The continental lithosphere extends no deeper than about 180 km 
(Anderson, 1979c). It may terminate at the boundary between granular 
and sheared Iherzolite nodules, 'vl50 km (Boyd and Nixon, 1975). We 
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will assurntt that the sheared ncdules are representative of the 
mantle below the lithosphere and above the Lehmann discontinuity at 
220 km. The shallower granular nodules have apparently been subjected 
to basalt extraction since they contain less CaO and AI 2 O 2 than the 
sheared variety. They may be an important, perhaps major, component 
of the continental lithosphere. Both varieties of nodules are enriched 
in the L.I.L. elements compared to oceanic crust, the MORfi source 
region and the minerals of ecloglte Inclusions. The fertile peridotite 
presumably rises to shallower depths under the oceans, of the order of 
80 km. Thus, the average thickness of the fertile peridotite layer is 
about 120 km. Volume trlcally, this is an adequate source region for 
continental and hotspot magmas but not for the more voluminous MORB. 

We suggest that differentiation of the Earth leads to two layers 
in the upper mantle, a thick basalt crust over residual peridotite. As 
the Earth cools the base of the original crust transforms to ecloglte 
wh? di sinks through the upper mantle. The present upper mantle is 
peridotite overlying a thick (450 km) ecloglte section. Partial melting 
in the ecloglte section allows material to escape and to melt extensively 
upon ascent. This is proposed as the source of oceanic crust. 

The basalt and ecloglte portions of the oceanic lithosphere return 
to the ecloglte layer by subductlon. The fertile peridotite layer 
of the upper mantle can partially melt and provide basalts when the 
upward convection of mantle heat is prevented by the insulation of 
continental lithosphere. Thus, both reservoirs are global and underlie 
both oceans and continents. It appears that the MORB source region can 
also provide magma to back-arc basins, perhaps when material is displaced 
by the descending slab. 
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THE ECLOGUE SOURCE REGION 

Pipe eclogices have a strong resemblance to oceanic thcleiites 
in both the major and trace elements . The Rb/K and other ratios in 
blmlnerallc eclogitc closely resemble the corresponding ratios in 
abyssal tholeilte basalts. The similarities are even more pronounced 
if the eclogites are compared With the average composition of the 
oceanic crust. 

The first column of Table 1 gives a composition which is 
representative of oceanic tholelites. More likely estimates of the 
composition of the primary magma are the total composition of the 
oceanic crust (column 2) and basaltic komatlites (column 3). These 
compositions are remarkably similar and have appreciably >}iore MgO 
and less Al^O^ and Ma 20 than tholelites which are considered 
to be the last crystallizing liquid from a more primary magma. The 
bimerali:.’ eclogites in kimberlite (columns 4, 5 and 8) are 
virtually identical to these estimates of the average composition of 
the oceanic crust. Trace element comparisons between kimberlite 
eclogites and abyssal tholelites are given in Table 2. Again 9 the 
correspondence is remarkable. 

It appears that material similar to eclogltC' inclusions in 
kimberlites is a suitable parent, for the oceanic crust. The inclusions 
themselves may represent cumulates from mantle diapirs that were trapped 
in the continental lithosphere. Diapirs rising from such great depth 
wov:].' melt extensively if their ascent were unimpeded .by the continent. 

THE GARNET-PERI DOmE LAYER 

The K, Rb, and Sr contents of some kimberlite garnet peridotite 
inclusions are given in Table 2. Also given are estimates of CFB and 
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of ths "pliuno" sourcA. Note the agreement between tholelitea and 

eclogites and between peridotltes and the inferred plume source region. 

Another v<xy to estimate the trace element content of a partial melt 

from peridotite is to assume that the difference in composition between 

fertile and barren peridotite is due to basalt removal. The trace 

element content of the resulting liquid is given in Table 3 and compared 

with continental and oceanic basalts. The peridotite compositions are 

from Rhodes and Dawson (1975) and it is assumed that the basalts 

represent 20% melting. These are extremely fresh peridotite xenoliths 

from the Lashame tuff-cone In northern Tanzania that have apparently come 

from a depth of ^.150 km. They are chemically and mlneraloglcally 

similar to peridotite inclusions from kimberlites except that they appear 

87 86 

to be relatively less contaminated. The Sr/ Sr ratios of these samples 
are about 0.705. The inferred melt is much higher in K, Rb, and Sr 
than oceanic tholeiites, a characteristic of continental basalts. The 
K/Rb And Rb/Sr ratios are also much different than abyssal basalts. 

Fertile garnet peridotite therefore seems a suitable, source material 
for continental f.lood basalts but not for MORB. It also has the 
characteristics inferred for the "plume** source region (White and 
Schilling, 1978). This part of the mantle has probably been subjected 
to metasomatlc enrichment of the incompatible trace elements (Lloyd 
and Bailey, 1975). Such enrichment has also been proposed for ?.he source 
region of continental (Boettcher and 0*Neil, 1979) and plume basalts 
(l^hite et al., 1979). The fact that enriched xenoliths are extensively 
sampled by kimberlites argues for the shallowness of the plume reservoir. 

Mid-ocean ridge basalts generally have ®^Sr/^^Sr ratios between 
about 0.702 and 0.704 while continental basalts are usually greater than 
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0,704 and range up to 0.710 (Garfar et al., 197'J; De Paolo, 1979), 

Basalts from oceanic Islands are Intermediate in value and may represent 

mixtures. The data on kimberlite xenollths is sparse and equivocal 

(Allsop ct al., 1969; Barrett, 1975; Simazu, 1975). Pipe perldbtltes have 
87 86 

Sr/ Sr values of 0.7060-0.7075 and other characteristics appropriate 

for the source region of continental basalts. Eclogite xenollths may have 

been brought into t-he continental lithosphere by deeper dlaplrs and evolved 

for some time in an environment different from PEL prior to pipe 

eruption. Whole rock measurements on eclogite xenollths from S. Africa 

87 86 

generally have high Sr/ Sr ratios (0.704-0.711). Allscp et al. (1969) 
estimated the ratio in "ideal" blmlneralic eclogite as 0.702. A 
sample from Tanzania has a value of 0.7004; discrete dlopslde nodules 
give 0.7029-0.704 (Barrett, 1975). 

LOCATION OP THE TWO SOURCE REGIONS 
As discussed earlier, at least part of the oceanic lithosphere seems 
to be returned to depths between 220 and 670 km. The mantle 
discontinuities at these depths are sharp and they are associated with 
changes in seismicity, as if they were acting as barriers to slab penetration. 
This could be due to density jumps caused by changes in mantle chemistry. 

The isotopic data, although useful in finger printing the source regions 
and giving age control does not provide information about major 
element chemistry or intrinsic density. This is \diere the kimberlite 
Inclusions become useful. 

Eclogite is appreciably denser than garnet peridotlte and should 
therefore occur deeper in a gravitationally stable mantle. I have suggested 
that the Lehmann discontinuity at 220 km is the boundary between garnet 
peridotite and eclogite, and the discontinuity at 670 km is the 

C_,. - 
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boundary botwoan acloglte In the garnetlte assemblage and perldotlte 
in the llmenlte plus spinel assemblage (Anderson, 1979b). The eclogite 
layer Is perched (PEL) in the upper mantle anc. forms the transition 
region. In a convecting system composed Of two superposed layers there 
is a thermal boundary layer, l.e., region of high thermal gradient, on 
each side of the interface. This is where temperatures are most 
likely to approach or exceed the solidus and where diaplrs would 
originate. There is also a thermal boundary layer associated with the 
lithosphere-asthenosphere boundary. Temperatures at 670 km and below 
are likely to be well removed from the melting point. 

Differentiation of a silicate planet results in two distinct products, 
basalt and residual perldotlte. The basalt, resulting from low pressure, 
high temperature partial melting of primordial mantle possibly resembling 
peridotiHc komatilte, would be originally concentrated in a thick 
layer at the surface. Subsequent cooling, primarily by convection, 
brings upper mantle temperatures into the stability field of eclogite 
which is denser than residual mantle. This leads to a massive 
overturning of the outer layers of the planet, subduction of the eclogite 
protoplate and destruction of the early geological record. Tills may 
explain the rarity of crustal rocks older than 3.8 Ga. 

Given that planetar' differentiation concentrates basalt in the 
outer layers and that the Earth has generated and subducted massive 
amounts of oceanic lithosphere what is the likely distribution of 
basaltic material in the interior? To answer this I have estimated the 
density as a function of depth for basalt and peridotite (Anderson, 

1979a). Basalt below about 50 km converts to eclogite which is denser 
than normal mantle even after olivine has converted to spinel. Pyroxene 
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and garnet react at higher presaurea to form a garnet solid solution. 

Noriiml mantle also undergoes a series of phase changes but remains less 
dense than garnetlte until ilmenlte and perovskite structures become 
stable below 670 km. The eclogite cannot sink below this levol. 

The addition of AI 2 O 3 CaSiO^ and (Mgi Fe)SiO^ expands the stability 
field of gar.ct and Increases the pressure required for transformation to 
such dense lower mantle phnses as perovskite and ilmenlte. Tills means that 
eclogite cannot sink into the lower mantle. Whole mantle convection, which 
may have been possible prior to the establishment of the eclogite layer 
and the chemical discontinuities in the upper nrnnt iu, would be replaced 
by separate convection systems in the lower mantle, the eclogite layer 
and the upper mantle above 220 km. 

The oceanic uart of the plate tectonic cycle in this scheme Is very 
simple (1'iGuve 1). Heating of the eclogite layer causes partial melting 
and the rise of eclogitic diapirs. The latent heat for complete melting 
is provii. cl hy adiabatic decompression. Oceanic crust forms from 
this eclogite liquid. MORB forms the surface veneer and represents the 
latest I'ceuzing fraction. Subduction causes the crust to reinvert to 
eclogite and it sinks back to the PEL. Tlie hn*zburgite part of 
the lithosphere remains in the upper mantle. The continuous recycling 
and rcmeltlng of the material in the oceanic crust depletes it in the 
Lit and, in particular, the LREE. 

PRIMITIVE MANTLE 

Tl»e isotopic data indicates that the two source regions are the 
results of an early differentiation event. If we accept the 220 and 
670 km discontinuities as its boundaries, the eclogite layer represents 
about 20% of the mass of the mantle. By assuming that the whole mantle 
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was involved In this early differentiation wa can obtain a spectrum of 
estimates of primitive mantle composition. Several of these are given 
in Table 4. 

There are other approaches that have been used for estimating 
primitive bulk Earth chemistry. Ganapathy and Anders (1974) have provided 
a cosmochemlcal mixing estimate which is also given in Table 4. There 
is surprisingly good agreement between these estimates and the resulting 
compositions are distinct from any modern rock type. Peridotite 
komatiites are widespread in early (>3.5 Ga) Precambrlan terrains. Vil,ioen 
and Viljoen (1969) propose that these approximate the composition rt 
primitive mantle. Indeed, these have Mg/Si ratios in the range of whole 
Earth estimates. They may represent primitive mantle that has left some 
garnet in the source region. I propose that a material similar to 
those in Table 4 was the parent from which the current mantle 
reservoirs were derived. These reservoirs are a shallow peridotitic 
layer and a deeper eclogite layer. In this scheme pyrollte would not 
represent primitive mantle but mantle which has already been 
depleted in a basaltic component. 

It has long been recognized that the source region of MORB is 
depleted in HL elements compared to alkali basalts, continental flood 
basalts, and hot-spot magmas. One would expect, however, that the 
original primary differentiation would enrich the basalt/eclogite fraction 
relative to the residual peridotiie. Whole rock analyses of pipe 
eclogite Indeed show such enrichment. Tlte major phasas, omphacite 
and garnet are, however, depleted and the enrichment occurs in the 
intergranular' material (Allsop ot al., 1969). The intergranular material 
LIE content is similar to that of the continental crust (Table 2). The 
eclogite layer may have become depleted ami the peridotite layer enriched 
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by the upward transport of fluids as discussed by Frey et al, (1978) » 
Boettcher et al. (1979) and Mysen (1979). 

Tatsurooto (1978) and Hedge (1978) proposed a model, based on lead 
and strontium Isotopes that Is similar to the present result, l.e., 
the LVL ("asthenosphere") Is undepleted or enriched and supplies "hot- 
spot" magmas; the underlying mantle ("mesosphere") is depleted and provides 
abyssal tholelltes. This Is the opposite of Schilling's (1973) model. 
Continental and hotspot related magmas represent a wide range of partial 
melting, from about 4% to 25% (Frey et al., 1978; White et al., 1979). 

Tltls suggests that they come from a wide range of (shallow) depths. 

MORB's are Invariably tholelltlc, indicating extensive (>25%) melting 

and a consistently deep origin. 

The Incompatible trace. elements in both source regions are 

enriched relative to recent estimates of bulk Earth composition 
(Ganapathy and finders, 1974). For example, if the lower mantle 
is Identical in composition to the peridotitlc upper mantle and If 

the silicate portion of the planet Is 21% eclog.lte and 
0.5% continental crust then the major elements are in agreement with 
bulk Earth estimates such as Ganapathy and Anders (1974) but such 
elements as K, Rb, and Sr are about a factor of 2-1/2 higher, fhls 
can be accounted for if the mantle below 670 km has transferred its 
incompatible trace elements to the crust and upper mantle. This 
presumably occurred during the early differentiation of the Earth and 
accompanied basalt extraction from primitive mantle. The calculated 
Rb/Sr ratio of the continental crust plus upper mantle (peridotlte 
plus eclogite) is 0.028. This Is also the value inferred for the 
bulk Earth (Ganapathy and Anders, 1974; DePaolo and Wasserburg, 1976). 
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This plus the complententary nature of the two source regions suggests 
that the material above 670 km may exhibit bulk Earth patterns of the 
1.1L elements, and have a 2-3 times enrichment. The lower mantle is 

extremely depleted in the LXL elements. 

EVOLUTION OF THE MANTLE 

The evolution of the Earth's mantle according to the present scheme 
is shown in Figure 2. The primitive mantle has roughly the composition 
given by a 1:4 mix of ecloglte and garnet peridotite. Early 
differentiation processes lead to the development of a thick enriched 
basalt crust and a residual peridotite mantle. As the outer layer cools 
it converts to ecloglte which settles through the upper mantle and becomes 
perched near 670 km by a phase bovmdary in the peridotite mantle. Incompatible 
trace elements are removed into the continental crust and the upper mantle 
garnet peridotite regions by raetasoraatic fluids. The eclogite layer 
therefore becomes depleted in those components which are concentrated In 
the overlying layers including the continental crust. 

The highest temperatures in the mantle, relative to melting 
temperatures, are in the thermal boundary layer near the top of the PEL, 

'\'220 km depth. This is where the density contrast between eclogite 
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and peridotite can ba overcome by partial malting and whara ecloglte 
dtaplrs originate. Peridotite dlaplra originate from the top of the 
thermal boundary layer. Their shallower depth and the broad malting 
Interval of peridotite leads to relatively small amounts of partial 
meltings a requirement of alkali basalt petrology (Frey et al.| 1978). 

This is a perauasi\*e but overlooked argument for a shallow location of the 
plume source region relative to the source region for tholelites. 

Continental and alkali basalts are usually emplaced at greater 
elevations than the oceanic tholelltlc basalts. This is sometines taken 
as evidence that the MORB source region is shallowar than the plume 
source region. Alternatively* the alkali basalts are emplaced at higher 
elevations because they are intrinsically less dense than tholelites. 

The pantm»cers of the thermal boundary layer depend on the thermal 
propertLof?, viscosity and heat flow at the interface. The thickness is 
calculated to be about 20 km and the temperature rise 

is 300“(-00"C. This is comparable to the near surface gradient and brings 
the average temperature close to the melting point of mantle silicates 
at 220 km. The thermal perturbation by a stationery continent* or a 
large continent moving slowly* causes a further temperature rise and may 
be. the trigger that Initiates partial melting. 

Depleted peridotite* the refractory product of partial melting of 
the garnet peridotite layer* is lighter than any other component of the 
mantle and becomes part of the continental lithosphere and the suboceanic 
harzburgite layer. The lower mantle need not be involved at all in 
the current magmatic cycle. 
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MANTLE METASOMATISM AND THE REDISTRIBUTION OF TRACE ELEMENTS 

Calculation of the partitioning of the rate Earth eXamenta among the 
varioua regions indicated In Figure 2 for a dry planet give results which 
are contrary to observations (Schilling, personal communication). In 
particular the eclogite layer being the result of partial melting of 
primitive mantle should be LREE enriched even after removal of 
the continental crust. The perldotlte layer should be depleted relative to 
primitive mantle and be LREE deficient. This suggests that there has been 
upward transport of a LREE phase which serves to deplete the PEL and 
enrich the overlying layers, including the peridotite plume source and 
the continental lithosphere. An ll20-rlch vapor phase strongly 
concentrates the REE and, in paticular, the LREE (Mysen, 1979). Evidence 
that the mantle has experienced such metasomatism prior to the transport 
of samples to the Earth’s surface by kimberlites has been provided by Ridley 
and Dawson (1975), Erlank and Shimazu (1977) and summarized by Mysen (1979). 
Parallel evidence from alkali basalts is given by Boettcher and O'Neil 
(1979). Carbonatites, kimberlites, "depleted” granular perldotltes, 
alkali basalts and the interstitial phase in eclogite xenoliths all show 
extreme LREE enrichment. The source region of MORB on the other hand is 
LREE depleted. The degree of LREE enrichment in a water-rich fluid 
increases with garnet content (Mysen, 1979). These observations all 
suggest that a vapor of fluid phase removes the REE from the eclogite 
layer and deposits them in the perldotlte layer and the continental 
lithosphere. Tlie other LIL elements also show a complimentary pattern 
between MORB and plume basalts. 

Since the enrichment of REE in general, and LREE in particular, in 
the fluid Increases rapidly with pressure as well as garnet content, we 
expect an upward increase in LREE enrichment in the rocks as we proceed 
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from the eclogite to the fertile perldotite to the "depleted** perldotite 
(continental lithosphere) and finally to the crustal layers • Mass 
balance calculations suggest that the continental crust is not the only 
repository of LREE enriched material. The remainder we propose is in 
the continental lithosphere and the hVZ. This is supported by the evidence 
from material which has been sampled from these regions. Tlie depleted 
nature of the MORB source region is an argument that it lies below the 
REE enriched regions, l.e., below the depth of generation of 
kimberlites and alkali basalts and the xenollths they contain. 

PLATE TECTONICS 

A major source of mantle heat flow is the Perched Eclogite Layer 

2 

(PEL). With thcleiitic concentrations it would provide 0.7 pcal/cm sec, 
about !/■’ of the global average surface heat flow. The scale length of 
convection in the PEL will he of the order of the layer thickness, n<450 
km. Convection In the overlying perldotite layer would be driven by this 
non-uniform heating from below and would have a similar or smaller scale 
length and be characterized by narrow ascending plumes. The 
high thermal gradient in the boundary layer leads to a large decrease 
in viscosity at the Interface and the two convectlng systems may be 
thermally coupled rather than coupled by viscous drag forces. That 
la, cold descending regions would occur in close proximity in the two 
layers. A descending slab, for example, may trigger detachment of the 
cold boundary layer in the PEL. This is shown schematically in Figure 3. 
Seismic waves will see a continuous cold region with high velocities. 

There will also be other regions of cold descending plumes in the PEL which 
are not directly related to slabs. Likewise, there may be regions 
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of rising curronCB In both the PEL and tho ovarlying layar» which do not 
express themselves in surface features such as mid~ocran ridges. The 
convective pattern, however, may be evident in detailed analyses of 
topography, gravity and seismicity. 

Earthquakes do not extend below about 250 km at most convergent 
plate botindaries. In other regions there is a gap in the seismic zone 

between about 250 and 500 km. Even where the zone is continuous 

it is usually contorted near 250-350 km. In many cases the deeper zone 
is more~or~less spatially continuous with the upper zone but in Chile, 

Peru and New Zealand the two zones appear to be displaced. This is 
suggestive of the type of two-layered convection considered here. 

There are several ways to estimate the lateral extent of the convection 
cells. We assume that convection in the ecloglte layer controls 
convection In the thinner and shallower peridotite layer. Thiessen et 

al. (1979) suggested that the distribution of high spots in Africa 

reflects the underlying convection pattern. By comparison with laboratory 
data they inferred a vertical extent of convection of about 500 km, 
slightly greater than the thickness of the ecloglte layer. They proposed 
that this pattern could only be observed through a stationery continent. 

Jordan (1978) showed that terrain, crustal thickness, and Bouj:;uer 
gravity anomalies have correlation distances of the order of 550 km, 
remarkably similar to African hlghspot distances. This again suggests 
a scale length of convection comparable to the thickness of the transition 
layer. 

Menard (1973) attributod depth anomalies in the eastern Pacific 
and the bobbing motion of drifting islands to convection cells in the 
upper mantle of half-wavelength 250-500 km. The depth anomalies, 
having amplitudes of ±300 meters can be explained by temperature 


A-IV-X8 


dlC;er«ncei in a 200 km -hick layer of about 200*C. The depth anomaUee 
seem to be fixed lelative to hot apotB. Menard believes that motion 
of plates over these bumps explain many aspects of vertical tectonics. 

Hie "bumpy" asthenosphere envisaged by Menard is a natural consequence 
of the convection pattern proposed here. 

The pattern of convection in the PEL is probably fairly complicated. 

The migration of trenches and continents may change the locations of 
the descending plumes in an individual cell or groups of cells in the 
PEL but it seems likely that the cells themselves cannot move far relative 
to one another. This provides a rationale for a fixed hot-spot reference 
frame and a mechanism for allowing the surface expression of a 
hot spot ^0 wander on the order of 5”. The proposed upper mantle 
convection pact 'rn is shown schematically in Figure 4. 

HOT SPOTS AND PLUMES 

Morg.Mi (197?) suggested that island and sea-mount chains are 
produced by plate motion over convective plumes extending from near the 
core-mantle Interface to the base of the lithosphere. Anderson 
(1915) proposed that plumes came from a distinctly different source 
region than midocean ridge basalts and that kimberlites, carbonatites 
and continental flood basalts were all related to hot-spot or plume 
activity. I suggested that plumes were a result of a thermal perturbation 
due cither to continental Insulation or a dlapir rising from the 
deep mantle. In either case the plume source region differed in 
chemistry from mantle providing abyssal tholelltes and was the source region 
for distinctive ocean Island magmas and undersaturated continental magmas. 

It was proposed that this source region was rich in Ti, Ba, Sr, Ys La, 

Zr, Nb, Rb, CO 2 , P 205 » H 2 O, etc., compared to "normal" mantle (the source 
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region of MORB)* Thoes charv,ct«ri»cics have elnce bean found Co be 
cy»"ical of plume chemlaCry (Unni and Schilling, 1978; Schilling ec al., 
1976; BonaCti oC al., 1977; Wtiite and S< hilling, 1978). 

Although 1 initially favored a deep mantle origin for plumea 
it now appeara Chat they originate above the 220 km mantle diecontinuiCy 
in a region of the mantle that has bsvn enriched in incompatible componenta 
by raeCaaomatic proceaaea that have depleted the source region of MORB. 

Since hot spot activity ia restricted both in space '^nd time we 
need a thermal anomaly to initiate partial molting in the peridotite 
layer. One possibility is themal blanketing by tho thick conductive 
continental lithosphere (Anderson, 1975). A large stationary 
continent may thereby cause its own break-up and a subseqnent period 
of rapid plote motion. The temperature anomaly and partial melt 
zone remains after the continent leaves and it becomes an oceanic hot 
spot. Hot spots will have a finite lifetime which appears to be at 

least 200 Ma after the continent starts to move off. 

2 9 

UiCh a mantle heat flow of 0.7 pcal/cm sec there are 2 x 10 

2 g 

cal /cm delivered to the base of subcontinental lithosphere in 10 
years. Even if only 20% of this is trapped this is enough to heat a 50 
km thick section of man.-^. .^y 200*C and to melt it to the extent of 
20%. Continental blanketing therefore seems to be an adequate 
mechanism for turning on melting spots below the continental lithosphere, 
The hot spot tracks in the Atlantic and Indian Oceans can be 
traced back to continental interiors. The timing of Mesozoic and 
Cenozoic continental flood basalts in North and South American, Africa, 
Europe and Siberia is appropriate for their location over hot spots when 
they formed (Morgan, 1979). If our hypothesis concerning the origin of 
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hot spoti undar stationary continants Is corract thsn va would expsct 
that a continent would have been ovar Hawaii soma 200 ai.y* ago. Tha 
Hawallan-Bmparor seamount chain disappears into the Aleutian Trench so 
wa cannot trace it beyond about 70 m.y. By backing up the Pacific plate 
we can infer that some of these continental fragments may have 
been incorporated into northwestern North America. The 210 m.y. old 
flood basalts from central Alaska to Northern Oregon, so-called 
Wrangallia (Jones et al. , 1977)^ originated near the equator (-v-lS*) and 
subsequently moved north to become attached to North America. Greenstones 
in central Japan were formed near the equator in the late Paleozoic 
(Hattori and Hlrooka, 1979). Other fragments are elsewhere around the 
Pacific Margin (Nur and Ben-Avraham, 1979). 

In the pftfuent scheme it is the ther perturbation caused by the 
deep ('v-i: ) km) continental lithosphere that is responsible for the onset 
of hotsp.i-: activity. The hotspots generate kimberlites, carbonatites, 
alkali has nit 3 , and lead to continental breakup. There follows a period 
of rapid continental drift and seafloor spreading. Continental igneous 
activity wanes as the continents drift off their hotspots but ridge and 
ocean Island volcanism Increase, In this scenario hotspots play an 
Important, perhaps dominant, role in breaking up and driving the plates. 

At the end of an interval of rapid spreading they would tend to be 
centr.ally located in the oceans, much as they are today. An Important 
force in plate tectonics may be the "hotspot fleeing force". The thermal 
nerturbatlon by continents may also control locations of mid-ocean 
ridges (Nur and Ben-Avraham, 1979). 
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HOT SPOT PROPULSION 

Elder (1976) has considered the propulsion of continents by a 
horizontal temperature gradient. In the absence of resisting forces 
■ at plate boundaries the velocity is ' 

3 

u ■ ogATh /3vf 

vrhere ot Is the coefficient of thermal expansion (3 x 10 ^/K), g 

3 2 

Is acceleration due to gravity (10 cm/sec )» AT Is the temperature 

anomaly (200K) , h is the layer thickness (100 km), v is the viscosity 
20 

(10 P) and A is the horizontal scale of continent (1000 km). The 
assigned values are just for the purpose of obtaining an order of magnitude 
estimate for u which turns out to be about 6 cm/yr. 

Periods of extensive continental magmatlsm are correlated on a 
global basis and seem to last of the order of 0.3 to 0.4 Ga (Windley, 

1977). They are separated by periods on the order of 0.7 to 1.0 Ga. 

Periods 6? rapid plate motion, or at least of rapid apparent polar wander, 
last for about 30-60 Ma (Gordon et al. , 1979). At a velocity of 10 cm/yr. 
this would lead to total displacements of 3000 to 6000 kw which are of 
continental and Inter-contlnental distances. The continents, on the 
average, then would come to rest far from their own or other continents' 
hot spots. The gestation period for forming a hot spot appears to be 
200 to 400 Ma and the lifetime estimated from the duration of continental 
magmatlsm and the duration of the subsequent hot spot track may be as much 
as .^OO Ma. Since convection is more efficient, through the oceanic mantle 
we expect that hotspots will start to dissipate as soon as their continents 
move off and that regions of high heat flow in the oceans would mark the 
previous locations of stationery or slowly moving continents. 
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Reconstruction of the continents Indicate that most of the 
Atlantic and Indian ocean hotspots were beneath continents from about 
100 to >350 m.y. ago. Tlie present African hotspots were apparently 
beneath Europe at the earlier time. 

SUMMARY AND DISCUSSION 

Isotopic evidence indicates that there are at least two source 
regions of basaltic magma in the mantle which have remained separate for 
the order of 1 to 2 Ga. The sub-continental and hotspot source region 
Is enriched In incompatible components compare'd to the source region for 
ntldocean ridge basalts. Eclogite and garnet perldotlte xenollths 
In kimberlite pipes seem to havf the appropriate characteristics 
to provide mid-ocean ridge basalts, and continental basalv 
respoctiV',' tv- The eclogite layer, the source of mid-ocean ridge 
basalts, is denser, and therefore deeper than the enriched layer. 

Melting iu iioth layers may result from the thermal insulation provided 
by the thick continental lithosphere. 

The latent heat for extensive melting of eclogite dlaplrs is available 
if rapid ascent of 160 km or greater is possible (Yoder, 1976). Tliis 
is leas likely under thick continental lithosphere than in the oceanic 
asthcnosphere. Therefore, oceanic tholeiites occur in the oceans 
and eclogite xenollths occur under continents. The temperature rise 
required for the initiation of garnet peridotite diapirs may require 
a long period of continental insulation. Therefore, hotspots start 
under continents and lead directly to uplift and breakup and 
provide the Initial driving force for continevital drift. 

It is possible that continental insulation is also required to 
initiate melting In the eclogite layer. The difference in geometric 
style between ridges and hotspots would then reflect the difference 
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in convection in thn ecologite and peridotite layerst rising sheets In 
tho internally heated eclogite layer and rising plumes in the overlying 
peridotite layer. 

One would expect the basalts formed in the primary differentiation 
of the earth to be LREB enriched and enriched in the Incompatible elements. 
Present MORB is coming from a source region which h4S been depleted 
In these elements. Continental and hotspot magmas are enriched in 
H 2 O, P> Cl, F, Rb, Sr, Ba, Tl, etc. Ultramaflc mantle xenollths 

are also enriched in these components , even those which have been 
depleted in a basaltic component. The minerals in eclogite xenollths 
are low In the incompatible elements but the intergranular material 
is enriched in a material having abundances similar to the continental 
crust. We suggest that LIL elements were removed from the MORB 
reservoir not only by extraction of the continental crust but also by 
removal of a fluid or vapor phase which has enriched the continental 
lithosphere and the upper mantle peridotite layer. 

A schematic continental geotherm with a thermal boundary layer is 
shown in Figure 5. The ascent path of an eclogite diapir is also shown. 
Voder (1976) has estimated that eclogite must rise about 160 km from 
its source region in order to completely melt if the heat of melting 
is obtained by adiabatic rise. For a 220 km deep source region complete 

melting would be achieved at 60 km. Note that completely molten 
eclogita, i.e., basaltic magma, can be delivered to the surface without 

melting dry garnet peridotite. With the geotherm shown, wet peridotite 

dlaplrs can rise from the thermal boundary layer but they will be 

only partially molten even after ascent to the surface because of the 

broad melting interval. Extensively molten peridotitic 

diapirs are not possible today because of their high liquidus temperature 

and their relatively shallow origin. This suggests that komatiites 
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could only form when peridotito dlnplrs could rise from depths greater 
than some 300 km, i.e., prior to the establishment of the thick ecloglte 
layer. Only such deep diaplrs could extensively melt (>60%) before 
they reach the 8ur*'.f',o. This assumes that the geotherm approaches the 
adlabat only at gri U-'i depths. 

The gradual heating associated with continental Insulation will 
mobilize mantle fluids before extensive melting occurs. Therefore, 
a metasomatlc precursor and a redistribution of LIL can be expected 
prior to continental magmatism. Alkali basalt activity can also be 
expected to precede and accompany tholelJtlc eruption from the deeper 
levels. 

We have suggested that a thermal perturbation may be required to 
initiate the rise of diaplrs from the MOIUJ source region. Tholeiltic 
volcanlsni is also associated with back-arc spreading. This suggests 
that suhil>‘iM: Ion can also trigger the rise of dlapirs from the MORB 
source. Thi'V will not have been preheated to the extent of normal MORB 
magmas and may ho more volatile rich. The spatial relationships of 
marginal basin and island arc volcanlcs, relative to the Benioff zone, 
suggest tliat the latter come from shallower depths than the back-arc 
tholoiites. Since island arcs are generally less than about 150 km above 
the Benioff zone (Rlngwood, 1975) this suggests that marginal basin basalts 
come from somewhat deeper, consistent with our interpretation of the 
location of the MORB source region. 
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Tabl« 1 

POSSIBLE COMPOSITIONS OF THE TRANSITION ZONE ECLOGITB UYBR 



(1) 

( 21 . 

(?) 

(4) 

(5) 

(6) 

(7) , 

(8) 

SiO^ 

50.3 

47.8 

46.2 

45.7 

49,5 

46.6 

45.5 

47.2 i 2.4 

TIO2 

1.2 

0.6 

0.7 

0.4 

0.5 

0.8 

1.9 

0.6 ± 0.3 

AI 2 O 3 

16.5 

12.1 

12.6 

17.9 

8.5 

13.7 

12.4 

13.9 ± 4.5 

FeO 

8.5 

9.0 

11.4 

11.2 

8.8 

9.1 

9.5 

11.0 ± 3.6 

MiO 

0.1 

0.1 

0.2 

0.2 

0.2 

0.2 

0.2 

0.2 ± 0.1 

MgO 

8.3 

17.8 

16.6 

11.9 

16.2 

16.1 

18.8 

14.3 ± 3.0 

CaO 

12.3 

11.2 

10.5 

7 A 

10.6 

11.8 

9.7 

10.1 t 2.2 

NajO 

2.6 

1.3 

1.2 

2.0 

1.7 

1.3 

1.6 

1.6 ± 1.1 

K,0 

0. >. 

0.03 

0.02 

0.4 

1.1 

0.02 

0.1 

0.5 ± 0.4 


\ (I) Oceanic tholelite (Kay et al, 1970, Engel and Engel, 1964) 

I 

(2) Oceanic cruat, calculated from ophlolite section (gxthon, 3979) 

(3) Basaltic "komatiite”, Gorgona Island (Gansaer ed al, 1979) 

( 4 > ) 

jgjj Mantle eclogites (Nixon, 1973) 

(6) Possible ecloglte extract in fractionation of primary magma in 

' upper mantle (O’Hara et al, 1975) 

(7) Picrite (Rlngwood, 1975) 

(8) Average biminerallc eclogite in kimberlite (Tto and Kennedy, 1974) 
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Tabl« 2 

ESTIMATES OP LITHOPHILE ELEMENT CONCENTRATIONS 
(ppm) IN BULK EARTH, THE ECLOGITE AND PERIDOTITE 
SOURCE REGIONS AND VARIOUS PRODUCTS OP THESE 
SOURCE REGIONS 


K Rb Sr Rb/Sr U Ref. 


Abyssal tholeiite 732 

Kimberlite ecloglte 820 

Kimberlite parldotlte 617 

Kimberlite peridotite 483 

'•Plume" source >468 

Continental tlood basalts 6400 
Ocean Island basalts 3160 

Continental crust 13 000 

Intergranular material 

in eclogites 16 000 


0.75 

92 

0.008 

0.16 

(1) 

0.7 

95 

0.007 

0.17 

(2) 

3.4 

55 

0.061 


(2) 

2.0 

59 

0.035 


(2) 

2.5 

60 

0.042 


(1) 

17 

320 

0.053 

0.3 

(3,4) 

5.3 

231 

0.023 


(1) 

33 

370 

0.089 

0.7 

(3) 

48 

550 

0.087 


(2) 
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TABLE 2 
(LEGEND) 

(1) WhiCft and Schilling (1978) 

(2) Allsop at al (1969) 

(3) Jacobsen and Waaaerburg (1979) 

(4) Camlchael at al (1974) 
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Table 3 

TRACE ELEMENTS IN INFERRED MIDOCEAN RIDGE BASALTS, 
ECLOGITE ffiNOLlTIIS, PARTIAL MELT OF "FERTILE” GARNET PERIDOTITE 
AND CONTINENTAL FLOOD BASALTS 


Oceanic Continental 



(1) 

(2) 

(3) 

(A) 

(5) 

(6) 

K 

820 

732 

700 

2000-3600 

4000 

6400 

Rb 

0.7 

0.75 

1.1 

11-20 

20 

17 

Sr 

95 

92 

134 

30-500 

400 

320 

K/Rb 

1170 

976 

640 

100-300 

400 

376 

Rb/Sr 

0.01 

0.008 

0.01 

0.04-0.3 

0.005 

0.053 

(1) 

"Ideal" 

ecloglte 

xenolith 

(Allsop et al, 

1969), 

Potential 


source region for oceanic tholelltes. 

(2) Oceanic tholellte (White and Schilling, 1978) 

(3) Oceanic tholelit'? (Table 1) 

(4) Inferred partial melt (20%) product from fertile garnet 
peridotite xenoliths with sterile peridotite xenollths 
as residual (Rhodes and Dawson, 1975) 

(5) Karoo basalts (Carmichael et al, 1974) 

(6) Continental flood basalts (Jacobsen and Wasserburg, 1979) 
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Table 4 

INCLUSIONS IN KIMBERLITES AND ESTIMATES 
OF PRIMITIVE MANTLE COMPOfITION 



Xonollths 


Mantle 

Compositions 



(1) 

(2) 

<3) 

(4) 

(5) 

(6) 

(7) 

Si'02 

47.2 

47.3 

47.3 

47.3 

48.0 

46.6 

44.8 

TIO 2 

0.6 

0.05 

0.2 

0.2 

0.3 

0.3 

0.2 

^^2°3 

13.9 

1.6 

4.1 

5.3 

5.2 

3.0 

5.3 

FeO 

11.0 

5.8 

6.8 

7.4 

7.9 

10.4 

10.3 

MgO 

14.3 

43.8 

37.9 

35.0 

34.3 

34,2 

34.3 

CaO 

10.1 

1.0 

2.8 

3.7 

4.2 

4.8 

4.4 

Na20 

l.n 

0.2 

0.5 

0.6 

0.3 

0.2 

0.4 

K 2 O 

0.5 

0.2 

0.2 

0.3 

0.2 

0.03 

0.03 


(1) Averaii;L» ccloglto nodules (Ito and Kennedy, 1974) 

(2) Average garnet: Ihzerolite in kimberlite (O'Hara et al, 1975) 

(3) 20% eelogite, 80% garnet Ihzerolite. 

(4) 30% eelogite, 70% garnet Ihzerolite. 

(5) Primitive mantle (Ganapathy and Anders, 1974) 

(6) Average poridotitic komatlite, S. Africa (Vlljoon and Viljoen, 1969) 

(7) Peridotite with quench texture, W. Australia (Nesbitt, 1972) 
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FIGURE CAPTIONS 

Figure 1. Schematic of the primary plate tectonic cycle with the 
transition zone eclogite layer serving as the source 
and sink of the oceanic lithosphere. The harzburglte 
portion of the lithosphere remains in the upper mantle. 

Figure 2. Flow chart of mantle differentiation. The primitive 
differentiation results In basalt and perldotite. 

Partial melting of the basalt layer concentrates the 
LIE elements Into the continental crust. 

Figure 3. Schematic illustration of isotherms for convection In a 
stratified system showing boundary layer detachment. 

The locations of descending plumes in the lower layer 
will be controlled by the locations of the cold iso- 
therms in the upper layer. Material penetrating Into the 
lower layer may Initiate dlapirlc uprise from this region. 

Figure 4. Illustration of flow in superposed convectlng layers. 

The lower layer Is Internally heated and is characterized 
by broad upwelllngs. The upper layer is mainly heated 
from below and is characterized by narrow ascending plumes. 
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Figure 5. Geotherm in a chemically stratified mantle showing the 
thermal boundary layer and the path of ascent of an 
eclogite dlapir that starts to melt at about 250 km* 
Temperatures at the top of the eclogite layer must appreciably 
exceed the solidus before the dlapir has stifflclent buoyancy 
to rise through the perldotlte layer. It may heat further if 
its upward escape is prevented by an overlying continental 
lithosphere. Extrusion temperatures and depth of complete 
melting may therefore be greater than shown here. Dry 
melUng curves are from Wyllie (1971) and Yoder (1976). 
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Abstract 

The high surface temperature of Venus Implies a permanently buoyant 
lithosphere and a thick basaltic crust* Terrestrial style tectonics 
with deep subduction and crustal recycling is not possible* 
Overthickened basaltic crust partially melts instead of converting to 
ecloglte. Because mantle magmas do not have convenient acesss to the 
surface the ^^Ar abundance in the atmosphere should be low* Venus 
provides an analog to Archean tectonics on the Earth* 

Introduction 

The surface temperature of Venus is about 450 K warmer than the 
surface of the Earth* This affects the buoyancy, thermal expansion, 
thermal conductivity and, hence, the thermal evolution and ultimate fate 
of the lithosphere. The pressure in the upper mantle of Venus is at 
least 12% less than at equivalent depths In the Earth's upper mantle. 
Thus, the depth of melting, the locations of upper mantle phase changes 
and the viscosity of the upper mantle will be different for the two 
planets. The buoyancy and thermal properties of the lithosphere control 
the style of plate tectonics and the associated time and length scales. 
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Mott dltcuttiont of th« comptrativt uctonict b«tw««n Earth and Vtnut 
addrttt only tha diffartncaa In vlacotlty* 

Tha ocaanlc lithoaphara of tha Earth coola at it agaa and it 
avantually bacomat dantar than tha undarlylng nantla* Thia Inatablllty 
davalopa aftar about 40 m.y. (Oxburgh and Paramantlar, 1977). Onca 

cooling roachaa a depth of 50 km tha garnat pyroxanlta or acloglta 
transformationa may make a aubatantlal contribution to the negative 
buoyancy of tha lithoaphera. Contlnanta and ocaanlc plateaua reslat 
subductlon because of their thick low-density crusts. The fate of the 
terrestrial lithosphere, therefore, depends on both chemistry and 
temperature. The purpose of this paper is to investigate the 

Implications for Venus tectonics of its high surface temperature. He 
will conclude that the surface thermal boundary layer on Venus Is 
permanently buoyant and that the reversible part of mantle convection 
occurs below about 100 km. 

The terms lithosphere and thermal boundary layer will be used for 

the cold outer layer of a planet. Neither term is strictly correct 

since the former has a strength connotation and the density of tha 
latter Is usually assummed to depend on temperature alone. The 
thickness of the mechanical or rheological lithosphere depends on 
temperature and, possibly, composition and Is proportional to the 
thickness of the thermal lithosphere that is discussed here. 

Lithospheric cooling 

The average temperature of the oceanic lithosphere decreases about 
660®C as It ages. The thickness of the conductively cooled thermal 
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boundary laytr Incraasta at a rate controllad by eht tharmal 

conductivity and tha diffaranca batwaan tha aurfaca and interior 
tanparatureif Since tha cruatal and harzburgita portions of tha 

boundary layer are both lass danse than the underlying mantle a 

gravitational instability <an only occur after the conductive cooling 
penetrates for a sufficient distance, ~30 km, into tha denser upper 
mantle portion of the boundary layer* In some petrological models the 
lower lithosphere is pyroxenltlc. Upon cooling, plagioclase and spinel 
pyroxenite convert to garnet pyroxenite or eclogite with a substantial 
increase in density. These phase changes also require cooling to depths 
in excess of 3C kilometers. 

The equilibrium thickness, 6, of a conductive boundary layer is 

6 - KAT/Q 

where K is the thermal conductivity, AT/ 6 is the average thermal 
gradient and Q is the mantle heat flow. For K of 7 x 10”^ cal/ctr, sec'C, 
AT of 1300“C and Q of 0.6 ycal/cra^sec, values appropriate for the 

terrestrial oceanic lithosphere, A is about 150 km. For a lithosphere 
composed of 6 km of basalt, 24 km of harzburglte and a lot*-" lithosphere 
composed of undepleted mantle, the equilibrium boundary layer for AT ■ 
1300®C is IZ denser than underlying mantle and therefore gravitationally 
unstable. Phase changes in the lower lithosphere may contribute further 
to the negative buoyancy. 

For a given mantle temperature and heat flow the equilibrium 
thickness of the boundary layer on Venus is reduced by two effects; K 
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if about 30% lower at the higher temperature (Schatt and Simmone> 1972) 
and AT if reduced by 460*C« Thie reducee d on Veiiuf to 42 km and Ite 
deneity» aefumming eiroilarity with oceanic lithoephere, if 2% lighter 
than the underlying mantlet Thie if a reeult of the amaller amount of 
cooling, the higher proportion of baealt and harzburgite in the boundary 
layer and the lack of fignificant cooling where it if moft required, 
ItOt the lower lithoephere# Furthermore, phaae tranfformatlone in the 
lower lithoephere cannot contribute at high temperature and low 
presBure# Baealt has about one-half the thermal conductivity of 
ultramaflc rocksi This reduces the thickness of the boundary layer to 
about 2? km if it is mainly basaltic# 

The buoyancy of the Venus lithosphere, even for a relatively thin 
crust, is greater than the buoyancy of young oceanic lithosphere on 
Earth. Since a thicker crust for Venus is probable, it Is certain that 
the surface thermal boundary layer for Venus is permanently buoyant and 
stable against subductlon# 

The temperature gradient in the conductive layer can be written 

dT/dz - (1/K)(q - Az) <2) 

B 

where C\ is the surface heat flow and A is the heat production rate in 
s 

the layer. With basaltic conductivities and radioactivities (A ~ 10**^^ 
cal/cm^ sec) the temperature and thermal gradient at 20 km are 950 and 
23‘C/km respectively and the solidus of dry basalt will be exceeded at 
depths shallower than 40 km. Eclogite is not stable in the mantle of 
Venus at depths shallower than about 100 km (Anderson, 1980). 


A-V-4 


Since crust cannot subduct and the garnet-rich asseoblages are vwll 
below the boundary layer, all partial melt products of mantle 
differentation that make their way into the outer 100 km or so will 
remain there. The crust is therefore likely to be thicker than either 
the oceanic crust or the average crustal thickness on Earth. It may be 
thicker than either the mechanical lithosphere or the thermal boundary 
layer. Crust which is overthickened by compression, buckling, 
thrusting, or collision will melt at its base. The resulting plutonism 
will further Increase the density contrast between crust and mantle and 
Increase the radioactivity and Incompatible trace element content of the 
parts of the crust so affected. 

The characteristic time for thickening of the boundary layer is 

X - (pCp/4 K)s2 (3) 

which for p - 3.,‘ g/cm^, Cp - 0.25 cal/g*C and 6 - 150 km is ~200 m.y. 
for the Earth. For Venus the characteristic time is 20 m.y. or less. 

Therefore, Venus has a thick, low-density and permanently buoyant 
crust-lithosphere that very quickly reaches thermal equilibrium. Upper 
mantle temperatures are high, resulting in low viscosities. 

The basalt-eclogite transformation will occur below about 100 km. 
The driving mechanism for mantle convection may therefore be partly 
chemical. Basaltic melt rising in an upeurrent will freeze as eclogite 
if it remains below about 100 km, and provide the negative buoyancy 
required for overturn. The 670 km discontinuity which halts eclogite 
subduction on Earth (Anderson, 1979a, c) will occur at about 800 km in 
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Venus giving a convectlng layer about 700 km thick. This should 
generate topographic and gravity anomalies of about 1400 km wavelength. 

Equation (1) suggests that plate tectonics in the Archean on Earth 
may have been similar to present day tectonics on Venus. The small 6 in 

s 

early Earth history was primarily due to the high Q from radioactive 

f 

decay, which was at least three times larger 4 x 10^ years ago. The 
thermal boundary layer thickens as Q decreases until eventually it 
becomes unstable and overturns, destroying the early geological record 
and setting the stage for the present style of tectonics (Anderson, 
1979a, b, 1980). This overturn event would allow much of the accumulated 
^®Ar in the upper mantle to escape and the subsequent steady-state 
rldge/trench style of ^lectonlcs allows continuous outgasslng of the 
terrestrial mantle. 

Early Evolution of Earth and Venus 

The initial evolution of Earth and Venus were likely quite similar. 
The high accretlonal energies would result in partial melting of the 
mantle and upward transport of a picritir. melt. Crystallization of this 
upper mantle melt layer, or magma ocean, yields a thin plagioclase-rlch 
crust and a deep ecloglte cumulate layer. An olivine-orthopyroxene 
cumulate layer may form at intermediate depths. This kind of scenerlo 
is well documented for the Moon except that pressures are too low for 
extensive ecloglte fractionation. Ecloglte fractionation removes AI 2 O 3 
from the melt and reduces the thickness of an early anorthositic crust. 
Ecloglte cumulates sink no deeper than 670 km on Earth because of the 
Intervention of the ilraenite and perovskite phase changes in a 
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peridotitlc mantle (Anderson, 1979c). An eclogite cumulate layer will 
be deeper and smaller in Venus because of the effect of temperature and 
pressure on the phase boundaries and the limited stability field of 
eclogite in the upper mantle. Nevertheless, the early tectonics, 
geochemical dlfferentation and outgasslng of the two planets were likely 
to have been similar. 

If the outer layer becomes unstable as the planet cools, the 
situation changes. Subducted slabs affect both the flow and thermal 
regime of the upper mantle. The sinking of the cold boundary layer deep 
into the upper mantle drags cold isotherms to depth and allows hot 
replacement material to rise at ridges. This leads to relatively large 
lateral and vertical temperature differences in the upper mantle and 
rapid convection in which the surface boundary layer participates. 

The mantle of the Earth contains approximately 15-20% of a basaltic 
component. If the mantle were well differentiated this would give a 
crustal layer 350-450 km in thickness. Only the top 50 km would be 
buoyant, because of the eclogite transformation, but this is still more 
than 4 times the present average crustal thickness. The luvtar and 
Martian crusts are much closer to being the appropriate thickness for a 
well differentiated planet in spite of the fact that accretional heating 
must have been less for these bodies. The difference, of course, is due 
to the continuous removal of basaltic crust from the surface implying a 
long-term storage reservoir at depth. The size of this reservoir has 
been estimated to represent about 20% of the mantle (Anderson, 1980). 

Samples thought to be representative of the upper mantle are 
depleted in AI 2 O 3 , CaO and S 102 compared to cosraochemlcal estimates of 
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mantle composition! They arc also depleted in trace elements which are 
retained by the ecloglte minerals, garnet and clonopyroxene. The degree 
of depletion Implies an ecloglte layer about 450 km thick which is 
probably deeper than 200 km. 

The ^®Ar abundance for Venus is about an order of magnitude less 
than for the Earth (Pollack and Black, 1979). This suggests that late 
outgassing has been less efficient for Venus than for the Earth in spite 
of the higher surface and upper mantle temperatures. This is easHy 
understood with the present model. On Venus, mantle melts and their 
volatiles mainly have access to the surface through relatively transient 
rifts. Because of the conservation of near surface material these rifts 
can only form when the crust compresses or thickens somewhere else on 
the planet. The amount of early outgassing on the two planets may have 
been similiar but Venus did not experience the early lithospheric 
overturn event or the outgassing asssoclated with continuous crustal 
renewal at ridge axes. 

Conclusions 

The high surface temperature of Venus has several Important 
tectonic implications. The most significant is the small amount of 
cooling that the lithosphere experiences before it reaches thermal 
equilibrium. This Increases its buoyancy and long-term siability at the 
surface when compared with the terrestrial oceanic lithosphere. Since 
crust cannot be recycled into the mantle the thickness of the basaltic 
crust on Venus is much greater than on Earth. The combination of a 
thick crust and high temperatures decreases the thermal conductivity of 
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the lithosphere. This means that for a given mantle heat flow the 
temperature gradient in the lithosphere is greater than would be the 
case for a cold lithosphere with a thin crust. This leads to high upper 
mantle temperatures, low viscosities and the possibility of partial 
melting at relatively shallow depths. 
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Abstract 


Trace element concentration patterns of continental and ccean 
Island basalts and of mld-*ocean ridge basalts are complementary. The 
relative sizes of the source regions for these fundamentally different 
basalt types can be estimated from the trace element 
enrichnent/depletion patterns < Their combined volume occupies most of 
the mantle above the 670 km discontinuity. The separate source regions 
are the result of early mantle differentiation and crystal fractionation 
from the resulting melt. The MORB source evolved from an eclogite 
cumulate which lost its late stage enriched fluids at various times to 
the shallower mantle and continental crust. The MORB source is 
primarily garnet and cllnopyroxene, while the continental and ocean 
island basalt source is a garnet peridotlte that has experienced 
secondary enrichment. These relationships are consistent with the 
evolution of a terrestrial magma ocean. 
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Introduction 


Hotspots, or plumes, have not yet been fitted satisfactorily Into 
either the tectonic or geochemical framework for the evolution of the 
mantle • From the point of view of mantle processes and Che chemical 
evolution of the mantle the alkali volcanlsm associated with hotspots Is 
a very Important process even if the abundance of such rocks Is small 
compared to the volume of abyssal tholelites (1)* The processes of 
creation and subducClon of oceanic crust and lithosphere clearly account 
for most of the maso transport into and out of Che upper mantle and the 
study of mid-ocean ridge basalts (MORB) has pla-'ed Important constraints 
on the part of the mantle that Is providing these magmas* Trace element 
and Isotopic considerations show, however. Chat, the source region for 
MORB cannot be representative of the average composition of the mantle 
nor can It have existed as a separate entity with its present 
characteristics for the full age of the Earth* The composition, 
location ami volume of this source region, which we designate MORBS, are 
still uncertain. It appears to be very large, homogeneous, global In 
extent, and to be the result of a previous differentiation or 
fractionation process that depleted It, relative to other source 
regions, of most of the Incompatible elements. On the other hand. It Is 
clearly not lacking In a basaltic component. The continental crust is 
complementary to MORB in both trace elements and isotopic ratios (2,3). 
Formation of the continents has presumably been at least partially 
responsible for the depletion of the so-called oceanic mantle. If the 
continental crust is the only enriched reservoir, then mass balance 
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calculations Indicate that only about 30 % of the mantle has been 
processed, the rest remaining undifferentiated and primordial (3,4). A 
large primitive reservoir would be difficult to reconcile with other 
evidence regarding the early thermal evolution of planetary interiors. 
There is, however, evidence for an additional enriched reservoir in the 
mantle. Magmas from continental interiors, rift aones, oceanic Islands, 
anomalous ridge segments and island arcs, mantle xenoliths from 
kimberlites and alkali basalts, and kimberlites themselves, all indicate 
the presence of a mantle reservoir with trace element concentrations 
that are also complementary to trace element concentrations in MORB. 
This reservoir apparently Is also global and can provide magmas to a 
variety of tectonic environments including all types of plate boun'’ .rles 
as well as plate interiors. These magmas, of which alkali basalts are 
one example, are not nearly as voluminous as MORB but they Indicate the 
presen^.^T nf un enriched region or layer which may be substantial 
compared to the continental crust. We shall designate this as the 
hotspot or plume reservoir, or simply PLUME (5), and assume that it is 
the source of alkali basalts, nephellnltes, melilitltes, basanltes, 
kimberlites and continental tholelltes which it provides by varying 
degrees of partial melting (6,7). One of the diagnostic signatures of 
these basalts is the rare-Earth element (REE) pattern which Implies a 
source enriched by about 10 in LREE and 3 in the HREE, relative to 
chondrites (6,7,8). 

PLUME basali.s also have high ratios of Rb/Sr, Rb/K, Ba/K, Nd/Sm, 
Ba/Nd, Ba/Sr, Nb/Zr, ®^Sr/®^Sr, and ^^^Nd/^^®Nd, and high concentrations 


A-VI-4 


of H 2 O, CO 2 , Tl, K, Rb, Sr, Ba, U, Nd, •«», and U ralatlve to MORBS. 
These characterlsttcH are shared by continental and ocean island basalts 
and are also evident, albeit diluted, in back-arc basins, Island nrcs, 
and transitional or anomalous ridge tholeiltes* Other characteristics 
that may be common but which have been studied in only a few hotspot 
locations are high abundances of Cl, F, Br and the primordial rare gases 
(9). The process that led to the enrichment of PLUMB will have depleted 
other, presumably deeper, regions of the mantle> 'flte absolute and 
relative sizes of these enriched and depleted reservoirs bear directly 
on the problem of the evolution of the mantle and, in particular, on the 
question of whether there is a large primitive reservoir In the mantle* 

Basalt s o urce regions 

The continental crust is extremely enriched In the Incompatible 
trace and minor elements. Continental tholeiltes, basanltes, 
nephelenites, alkali basalts, kimberlites and basalts from oceanic 
Islands are also enriched. The complement to the depleted mid-ocean 
ridge source region may, therefore, be much more voluminous than Just 
the continental crust. The trace element and Isotopic affinities of 
these rocks, regardless of their tectonic setting, suggests that they 
arc all derived from a similar source region. This reservoir, PLUME, 
has previously been referred to as continental mantle, the enriched 
source region or primitive mantle. Its location is uncertain but It 
appears to be a global layer. Proposals for its location Include the 
uppermost (6,10) and lowermost (11) mantle. In trace element ratios 
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such (IS Ba/Nd» Nd/l>m, Ba/Sr, Ba/K and Rb/K thara ie a prograsslva 
increasa from ocaan island and Island arc basalts, basalts from 
anomalous rldga segments, continental tholeiites, alkali basalts to 
kimberlites. Mid-ocean ridge basalts have much lower ratxos and 
tholeiites from back-arc basins are generally slightly higher than MORB. 
Trace element concentrations vary systematically through the 
compositional spectrum olivine melilltite, olivine nephelinite, 
basanite, alkali basalt, and olivine tholeiite and can be explained by 
varying degrees of partial melting ranging from 4 to 25% from a common 
source region which is enriched In the strongly incompatible elements 
(Ba, Sr, Til, U, LREE) by lOx chondrltlc and in the moderately 
Incompatible elements (Ti, Zr, Hf, Y, HREE) by a factor of about 3 
(6-8). Mid-ocean ridge basalts are derived from a very different source 
which is depleted in the large-ion lithophlle (LIL) elements. Tlie 
complementary nature of the t.?o source regions is possibly due to the 
mlgrattim from the MORB source region of a melt or fluid with 
Incompatible element concentrations similar to those of kimberlite. 
This fluid depletes the MORB source region and enriches the 
complementary mantle reservoir. The upper mantle low-velocity zone 
(LVZ) is a likely repository of these volatile and trace element 
enriched fluids and, therefore, a possible source region for the 
enriched magmas which are characteristic of the plume or hotspot source. 
Kimberlites come from depths as great as 220 km, consistent with an 
origin near the bottom of the LVZ. The xenollths in kimberlites and 
alkali basalts represent mantle fragments from various shallower depths 
and also generally exhibit enrichment of the type Inferred for the plume 
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Rourc« rsglon* Thers Is now a waalth ot avidanca lor upparmantla 
metasomatism (12, 30) and enrichment events at various timest The 
source of the enriching fluid has not been much discussssed. 

The complementary nature of continental and ocea : tholeiltes is 
illustrated in the center part of Figure 1. Alkalii basalts, 
nephellnites, melilitltes and basanltes have similar but more enriched 
patterns than continental flood basalts, consistent with their 
derivation from the same source region by smaller degees of partial 
melting* Continental and abyssal tholeiltes exhibit reflection symmetry 
about a line corresponding to a 6 to 7 fold enrichment over average 
mantle concentrations* 

Continental and ocean island basalts are relatively depleted in 
those elements that are retained by garnet and clinopyroxene* The 
reverse Is the case for abyssal tholeiltes. Since garnet and cpx are 
reduced or eliminated during the large amounts of partial melting 
usually Inferred for the for ation of tholeiltes (6), the inverse ga-cpx 
pattern of plume basalts must be related to a prior history involving 
ecloglte fractionation or invasion of PLUME by a fluid that was in 
equilibrium with a ga-cpx assemblage* 

The origin of island arc basalts is still controversial but they 
have PLUME affinities in such trace element ratios as Rb/K, Rb/Sr, Ba/K, 
Ba/Sr and Ba/Nd. These ratios are closer to ocean island basalts, 
alkali basalts and continental flood basalts than they are to mid-ocean 
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ridge basalts* Back-arc basin basalts are Intermediate to MORB and 
hotspot magmas, suggesting a mixture from a deeper MORB source and an 
overlying enriched source region. Ocean Island tholelites also appear 
to be mixtures both in trace element and Isotopic ratios (2). The fact 
that both MORB and basalts with continental affinities are available at 
ridge enrivonments, island arcs, oceanic Islands and well develop ' 
continental rifts suggests a composltionally stratified mantle* The 
order of appearance of these basalt types is consistent with the 
shallower layer being enriched lu the Incompatible trace elements* In a 
stratified mantle the deeper layer should be denser and, therefore, more 
garnet rich. The ratio of concentrations in MORB to those in 
continental tholelites indicates that the MORB source rugion, although 
depleted In most of the incompatible trace elements, has selectively 
retained thosu that are most comfortable In the garnet and cllnopyroxene 
lattice!'. Ttils suggests that MORBS is an ecloglte or garnet pyroxenite 
cumulale. 

The coraivosltlon and volumes of the two source regions 

The concentrations of some key incompatible trace elements In 
various magma >ypee, normalized to average mantle concentrations, are 
shovm in Figure 1* Relative to the average mantle, continental 
tholelites, alkali basalts, kimberlites and the continental crust are 
all most enriched In K, Rb, Ba, La and U* These are the elements that 
are most discriminated agaxnst by the major mantle minerals, olivine, 
orthopyroxene, cllnopyroxene and garnet, i.e., they have the lowest 
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mineral/melt partition coefficients. The other elements are also 
strongly rejected by olivine and orthopyroxenu» having partition 
coefficients of lees than 0.03 for these minerals (6). On the other 
hand, garnet and clinopyroxene have partition coefficients greater than 
about 0.1 for Y, Nd, Sm, Sr, and Yb. These are the least enriched 
elements in the magmas with continental and ocean island affinities and 
the most enriched in MORB. It appears titat garnet and clinopyroxene are 
abundant minerals in the MORB source region, and that the fluid which 
enriched PhUME was previously in equilibrium with a ga-cpx rich 
reservoir. 

The high concentration of Y and Yb in MORB relative to other 
elements and relative to concentrations in basalts from PLUME, the 
enrichel source region, is particularly significant. These elements 
have minora^ /'melt partition coefficients greater than unity for garnet. 
Mass balance calculations suggest that most of the terrestrial inventory 
of Y and Yb may be in MORBS. This can only he accomplished if garnet is 
a dominant phase and if most or all of the mantle has experienced 
differentiation. Partial melting of primitive mantle concentrates Y and 
Yb into the melt; crystallization of this melt at moderate presssurc 
would concentrate these elements into an eclogite cumulate. 

The relative volumes of the two source regions can be inferred from 
Figure 1 if it is assumed that they, plus the continental crust, give an 
undifferentiated terrestrial pattern for the very incompatible elements. 
Using 0.56% for the mass of continental crust relative to the mantle (3) 
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we compute that the mass of MORB must exceed the tholeiitlc fraction of 
PLUME by a factor of 12 ,* ’ Continental tholelltes represent about 
20% melting of a perldotitlc source region (6) and MORB must be mixed 
with at least 15% olivine to recover the composition of Its parent magma 
(13,lA)t The residual crystals in both cases are assumed to be depleted 
because of their very low partition coefficients (6,8). The MORB source 
region is, therefore, at least 2.8 to 3.5 times the mass of PLUME and 
the enrichment of the combined source regions, relative to primitive 
mantle, is about 3.3 to 3.8. Note that MORB is depleted only in a 
relative sense. Compared to average mantle abundances it Is enriched. 
This suggests chat MORBS plus PLUME are complementary to the remainder 
of the numtle from which they have presumably been removed by partial 
melting. The enrichment factor implies whole mantle differentiation and 
a residual mantle about 2-1/2 times the size of the combined reservoirs. 

By assuming that the very incompatible elements have been entirely 
fractionated into the two source regions It Is possible to estimate 
their totval mass. For example, if La Is depleted in the lower mantle 
the combined source regions represent 26 to 30% of the mantle. 'us 
corresponds to a thickness of 560 to 640 km in the upper mantle, or a 
region extending upwards from the 670 km discontinuity to a depth of 27 
to 110 km. If we take the transition region, 220-670 km, to be the 
depleted source region (20), then PLUME is 130-160 km thick, about the 
average thickness of the low-velocity zone (LVZ). 


A similar calculation for the other strongly incompatible elements 
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(K, Rb, Ba, U) gives a range of 21-38% for MORBS plus PLUME. These are 
upper bounds since It has been assumed that the whole mantle has been 
processed and depleted. Thls« however, Is a good approximation If 
olivine and orthopyroxene are the main residual phases. It Is 
significant that the above estimates of the relative and absolute sizes 
of tlie principal mantle reservoirs correspond to the main subdivisions 
of the upper mantle, the LVZ and the transition region. The 670 km 
depth corresponds to a major seismic discontinuity and the maximum depth 
of earthquakes. 

Other elements can be used to estimate the sizes of the Individual 
reservoirs. The partition coefficients of Yh and Y are such that they 
strongly preft-r the garnet structure (15). Both are enriched In MORE 
relative to PLUME, suggesting a higher portion of garnet in the former. 
If the f*nfcire mantle comi^lement of Y and Yb reside in )R)RBS, then this 
would represent 14-17% of the mantle or an upper mantle equivalent 
thickness of 300 to 364 km. Sc Is partitioned into garnet and 
clinopyroxene by a factor of 8 relative to olivine and orthopyroxene . 
Using this ratio of enrichment for the MORB source relative to the rest 
of the mantle it would constitute 26% of the mantle. 

K, Rb, Ba, I.a and U have the highest relative concentrations in 
PLUME. These are the elements with the lowest partition coefficients 
for the major mantle minerals and those that are most likely to be 
concentrated Into PLUME by melts or metasomatlc fluids. The PLUME 
abundances of these elements can account for the entire mantle inventory 
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tf I’l.UME Is 9 to 19% of the mantle. 

Usina the above estimates of the relative sizes of the two source 
regions and Ganapathy and Anders (37) values for terrestrial abundances » 
we can now estimate the total abundances of the trace elements in the 
various reservoirs. The continental crust contains more than 26% of the 
mantle plus crust inventory of K, Rb, and Ba and less than 8% of Sr, Y, 
Zr, Nb, Sm, Yb, and Ilf. The "depleted" MORB source contains more than 
50% of the Earth’s inventory of Y, Zr, Nd, Sm, Yb, Hf and Th and more 
than 20% of Nb, La and U. The PLUME source region, although generating 
highly enriched magmas contains only about 10 to 20% of the Earth’s K, 
Rb, Sr, Nd, Sm, Yb, Th, and U. It is the primary repository for Ba and 
La and has small total abundances of Y, Zr, Nb, and Hf (less than 7%). 
The crust, and the two upper mantle reservoirs account for about 75% of 
the hear, produccion from K, Th, and U. Therefore, only 25% of the 
Earth's heat flow comes from the lower mantle. About A0% is due to 
M0R.BS. The continental crust and PLUME each contribute about 18%. 

This distribution of heat sources affects the style of convection 
In the two reservoirs. Assuraraing PLUME to be the shallower reservoir it 
is primarily heated from below and will therefore be characterized by 
xxiwtou ascending plumes. MORBS is primarily heated from within and will 
be characterized by broader ascending regions and narrow descending jets 


or "slabs" 
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Thft mtneraloRv of the two source regions 

Assuming that the two source regions are related • an apparent 
partition coefficient can be dete*^,lned by forming the ratio of the 
concentrations In MORB to those In continental tholelltes (CFB). This 
Is shown In Figure 2 along with mlneral/roelt partition coefficients for 
garnet, cllnopyroxene, orthopyroxene and olivine. Tholelltes are used 
In this comparison since they presumably are products of large degrees 
of partial melting. They will, therefore, have the trace element 
pattern of their source and a relatively uniform enrichment. The ratio 
of the concentrations In the two reservoirs can be explained If the MORB 
source Is composed mainly of garnet and cllnopyroxene and PLUME has been 
enriched with fluids that were In equilibrium with this layer and 
therefore depleted In such elements as Zr, Nb, Y, Yb, and Hf. 

MORB Is enriched with respect to continental tholelltes In Y and Zr 
and only slightly depleted In Hf and Yb. Both MORB and continental 
tholelltes represent rather large degrees of partial melting of their 
respective source region. A garnet perldotlte would be expected to lose 
most of Its garnet and cllnopyroxene under these circumstances. The 
above results Indicate that primary garnet must be a minor mineral In 
PLUME but a major component of the MORB source. MORB is depleted in the 
most Incompatible elements, such as K, Rb, Ba, and U but not as much as 
would be expected if It were composed entirely of residual crystals. 
This suggests that MORBS evolved from a melt, presumably the result of 
an early differentation event, and consists of a mixture of melt with 
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excess garnet which has settled into it during crystallisation, i.e. it 
is an orthocumulate. The upper part of the nantle would consist of 
lighter, cumulate and residual, crystals and late stage fluids which 
impart an Inverse garnet signature to this region. 

This sequence of events, based on trace elements, is precisely what 
et al. (14) proposed on the basis of petrological and major 
element considerations. They concluded that the parent magma for ocean 
island tholeiltes had experienced a previous history of ecloglte 
extraction. The remaining magma evolved to ocean island tholelite by 
olivine fractionation. On the basis of major element (20) and trace 
element chemistry, I suggest that the ecloglte cumulates constitute the 
source region for mid-ocean ridge basalts. The parent magma from which 
both MORBS and PLUME evolved by crystal fractionation would be plcritlc, 
the result of extensive partial melting of a primitive garnet peridotlte 
mantle. Similar considerations have led to the concept of a magma ocean 
on the moon (41). Because of the higher pressures in the Earth's mantle 
a deep ecloglte cumulate is the analog of the floating plagloclase 
cumulate that forms the lunar highland crust. 

In order to see if this holds up quantitatively for the trace 
elements we investtgcte the following model. The primitive mantle is 
split into two reservoirs by 15% partial melting and melt extraction to 
the surface. The melt fraction reservoir is enriched by about a factor 
of 6 in the Incompatible elements but, because of the low partition 
coefficients, has nearly primitive ratios of Rb/Sr and Sm/Nd. 
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Crystnlilzation of the melt yields an eclogite cumulate layer which is 
modelled as a 50 j 50 nix of garnet and melt, proportions appropriate for 
an ovthocumulate. Part of the final 5% tnelt fraction, a malt that is in 
equilibrium with eclogite, is removed to deplete the cumulate layer and 
enrich the overlying mantle* The mixing ratios were adjusted so that 
the former can yield MORB by 15% olivine fractionation and the latter 
can yield CFB by 20% partial melting (38). 

Results of this model are shown as solid points in Figure 2. They 
are consistent with MORB being the result of nearly complete melting of 
a source region which represents high pressure cumulates from the melt 
fraction of the primary differentiation, l.e. eclogite or garnet 
pyroxenite. 

Prior to providing depleted MORB this region lost its late stage 
fluids or, alternatively, crystallized completely and then lost its 
early melt fraction on a subsequent reheating cycle. PLUME was enriched 
by upward migration of this fluid. Thus, two stages of dif ferentatlon 
and fractionation are required in order to concentrate the LIL in the 
t\w upper mantle reservoirs. Isotoplcally, the MORB source evolves as 
nearly primitive mantle until it transfers its incompatible elements to 
PLUME. It subsequently evolves as a depleted reservoir. Basalts from 
PLUME will have time-integrated depleted, enriched or ”p’^i™ibive” 
isotopic ratios, depending on when the enrichment occurred. 
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Locrttions o f the two Bource r egions 

There is variety of evidence that suggests that the plume source 
region is shallow. Anomalously low seismic velocities can he traced to 
a depth of about 250 kra under Yellowstone but deeper velocities appear 
to be normal (17). The variable and small amounts of partial melting 
required to generate melllltitcs, nephelenites, basanltes and alkali 
basalts (6) from a coimnon source region are consistent with dlapirs 
rising from various shallow depths. Tholeiltcs are the result of larger 
degrees of melting, consistent with adiabatic ascent from deeper levels. 
The xenoliths entrained in kimberlite,'* and alkali basalts are volatile 
and trace element rich and these are samples from the upper 200 km of 
the muitle (12). Intviitively, we expect that volatiles will migrate 
upward and he trapped by the cold upper mantle. Volatile and LIL 
enriched magmas occur not only under continents but also at Island arcs 
where tlie downgoing slab perturbs the upper mantle. Plume type basalts 
occur in continental rifts. These rifts evolve to oceanic ridges, with 
Isolated oceanic Island hotspots, when the rifting has led to the 
formation of an ocean basin. 

There are also sonui suggestive geometric constraints. The volume 
of the low-velocity zone (I.VS) is adequate to provide the plume basalts 
but not the voluminous MORB (10). The preferred mode of convection in a 
region with a high temperature gradient and a rapidly varying viscosity 
takes the form of hexagonal cells wltli upwelling centers (18). This 
seems to be consistent with hotspot patterns (19) and the small areal 


A-Vl-16 


t'xtiMVl of hotspots. In rooro homogen>i»n'« rwaions of the mantlet where 
the tampetritiire gradient Is smaller and the viscosity more uniform, 
linear rolls are a possible mode of heat transport. Linear ridge 
systems and the uniformity of oceanic tholelltas therefore suggest a 
deeper source, one that is below the large vertical and lateral 
variations that occur In tie upper 200 km. The sink of oceanic 
lithosphere, judging from the depth distribution of earthquakes, appears 
to be between 200 km and 670 km. The distinctive Isotopic and trace 
element signature of MORBS can be maintained if It is also between these 
depths. Upper mantle temperatures are closest to the melting point 
between about 150 and 250 km depth and this is therefore the depth range 
where it is most likely for diapirs to originate. 

It is difficult to estimate the depth at which partial melting 
first occurs. Seismic data from several hotspots give low velocities 
down to at least 150 km (17,21). Eruption temperatures of 1300“C and 
large degrees of partial melting require initiation of melting below 
150-200 km (22). There is some evidence for crystal fractionation in 
melts as deep as 280 km (23). Mantle fragments brought up by alkali 
basalts and kimberlites do not, in general, have the trace element 
pattern required for the source region of MORB. Therefore, PLUME magmas 
may originate from shallow depths but MORBS must be deeper than 200 km. 
In a gravitationally stratified mantle one would expect the deeper layer 
to be denser and tlierefore rich in garnet. The inability of young 
buoyant lithosphere to subduct below about 220 km suggests a density 
Increase at this depth and this may be the boundary between the two 
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source regions (20). Taking ai.1 evidence Into account it appears that 
PLUME nay be coincident with the LVZ and dlaplrs rising from the 
transition region supply the magmas that evolve to MORE. 

A terrestrial magma ocean 

The complementary LIL element patterns of the major terrestrial 
magma types is reminiscent of lunar data which has led to the widely 
accepted concept of crystal fractionation and cumulate formation in a 
magma ocean or vast lava lakes (Al). The various source regions are 
attributed to cumulate and residual fluid layers that resulted from 
crystallization of a ~300 km thick magma ocean that» in turn, was 
derived from very early melting of at least half of the moon. If a body 
as small aa the moon experienced such extensive differentiation and 
fractionation, then the Earth should have as well. The amount of 
partial melting required, ~15%, to explain the enrichment of the upper 
mantle reservoirs on Earth, is, in fact, relatively modest. Although 
the energy of accretion of the Earth is much greater than that of the 
moon, the greater size of the Earth results in high Rayleigh number 
convection and rapid increase of the melting curve with depth. This 
plus the latent heat buffer and the high melting temperature of olivine 
may prevent more extensive melting on a global scale. A 15 % melt 
implies a magma ocean ~400 km deep. The pressures in the Earth are 
greater than on the moon. Early eclogite fractionation at depth 
therefore preempts the extensive plagloclase fractionation that resulted 
in the early lunar anorthositic crust. 
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CryHtiilllsation of a magma ocean will {trocoed from the baa« because 
of thu rolailva slopes of the adlabat and the liquldna. The near 
liquiduB phases at depthrt greater than 60 km ate garnet and 
cliuopyroxane (14, 16) and they will form an ecloglte or garnet 
pyroxenlte cumulate layer. Although there Is a seismic discontinuity In 
the mantle near 400 km, ecloglte Is denser than residual garnet 
perldotlte to depths of 670 km (20). These cumulates vlll therefore 
sink to this depth, displacing residual mantle upwards* The shallower 
part of the mantle will therefore consist of olivine and pyroxene 
cumulates and residual fluid from the magma ocean mul displaced residual 
lower mantle. The top of the ecloglte layer could be as shallow as 
220 km (20) or as deep as 370 km based on considerations of this paper* 
The mantle discontinuity near 400 km, which is usually attributed to the 
olivine-spine I phase change could therefore be either a chemical 
discout 'nulty or could represent the completion of the garnet pyroxenlte 
(ecloglte) to garnetlte (garnet solid-solution) phase change (42). 

Isotopic ratios of the two source reRlons 

The study of Sr and Nd Isotopes places important time constraints 
on the evolution of mantle reservoirs (24-26). The continental crust 
nnd the depleted reservoir have mean ages of 1.5 b.y. (3). Kimberlites 
and continental flood basalts, although enriched in LIL, have been 
attributed to a primitive reservoir (2, 3, 24, 27). This Interpretation 
is not required by the data as we show later. Some mantle samples liave 
been din’lved from ancient enriched reservoirs (28). Other samples come 
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from rfiBfirvoirs which apparsntly hava baan anrlchad only racantly 

(2n,?.9). 

Consider a prlmltlva mantle that partially malts and separates Into 
twi) reservoirs » as before » by upward reuoval of the melt* The lower 
mantle consists of residual crystals and Is therefore depleted In the 
Incompatible elements. The melt fractionates Into a deep ecloglte 
cumulate layer and a shallow ollvlne-orthopyroxene cumulate layer. A 5% 
melt fraction Is transferred from the deeper to the shallower layer at 
various times (40). The Nd and Sr Isotopic ratios for the two 
reservoirs are shown In Figure 3. The central horizontal scale gives 
the ages of th.-* depletion and enrichment events. The theoretical 
^Nd'^Sr ‘-orrelation lines agree with the data and Indicate that the 
mid-ocean ridge basalt source region was depleted at times between 1.5 
and 2.5 x 10*^ years ago. A redistribution of LIL makes It possible to 
satisfy mantle Isotopic data even If the primary differentiation 
occurred early In the history of the Earth. Thus» the trace element, 
petrological and thermal constraints on the evolution of the mantle and 
the various reservoirs are not contradicted by the Isotopic data. The 
type of model investigated here is similar to previous Ideas of mantle 
ru'tasoraatlsm and trace element redistribution (30). 

Implications f or mantle evoluti on 


Isotopic studies Indicate that the two major mantle reservoirs have 
been isolated for more than 1 or 2 b.y. (2,3). The evidence presented 
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liar* from tr«c« el«mtntO| «nd previously from major alaroent and aalamlc 
coDsidaratlons <20)» auggaat that tli« innntla is chamlcally stratlflad 
and tiuit tha various regions of the mantle are complementary products of 
terrHstrlal differentiation* The present emphasis has been on magmas 
from these reservoirs but similar conclusions result from the study of 
solid fragments from the mantle (14» 20)* O'Hara et al* (14) pointed 
out that the restricted number of phases found in biminerallc eclogites 
from kimberlites indicate that they either crystal accumulates or 
crystalline residue developed In contact with a liquid* This is exactly 
the situation inferred from the trace element patterns for the source 
region of abyssal tholelites. Eclogite and garnet peridotite xenoliths 
from kimberlites may represent samples from MORBS and PLUME, respec** 
tlvely. They give satisfactory average mantle compositions for the 
major oxides when combined in the portions indicated by the trace 
elements (20). 

Although whole mantle convection, in the conventional sense, is 
precluded in a cliemlcally stratified mantle, transfer of material into 
and out of the various reservoirs is possible because of the large 
volume change associated with partial melting and phase changes, e.g., 
basalt-eclogite* The oceanic part of the plate tectonic cycle may be 
summarized as follovrs* Partial melting in the ecloglte layer allows 
dlapirs to rise to the base of the oceanic llthospViere. Nearly complete 
melting occurs during adiabatic ascent. This is possible because of the 
proximity, in temperature, of the liquidus and the solidus in ecloglte* 
Peridotite diaplrs can only partially melt because of their high 
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liquldl. The resulting melt fractionates in near surface magma chambers 
to form thole iitic melts which are light enough to rise to the surface. 
The pyroxenite residu forms the lower oceanic lithosphere. As the 
lithosphere cools and thickens the lower part transforms to garnet 
pyroxenite or eclogite which is denser than the underlying mantle. The 
oceanic lithosphere becomes gravitationally unstable and it returns to 
its source region which lies between about 220 or deeper and 670 km 
deptii. Partial melting in the shallow enriched perldotlte layer 
generates continental and ocean island basalts and a harzburglte 
residue, both of which are lighter than their source region. These 
PL7JME products remain in the crust and upper mantle. Some volatiles are 
returned to PLUNE by subducted sediments and hydrothernally altered 
oceanic crust. The high temperature gradient in the thermal boundary 
layer at the PLUME-MORBS Interface brings temperatures there closo to 
the melting point (20). When mature oceanic lithosphere passes over a 
hotspot it is quite likely that the lower pyroxenite part will 
contribute xenollths to the alkallc magmas and, by remelting, contribute 
to the trace element and isotopic signatures of ocean Island basalts. 
The evolution of the mantle Ig shown schematically in Figure 4. 

It has generally been assumed that basalts of all kinds represent 
partial melts of peridotites. Mantle compositions based on this premise 
have much lower abundances of Si02, AI2O3, and CaO than models based on 
coamochemical considerations. The evidence used to construct the 
petrological models is obtained from the upper 200 km of the mantle and 
the models are therefore btrictly only valid for the uppermost mantle. 
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A thick eclogite layer perveu to increase the abundances of the above 
components and can reconcile the cosmochemical and petrological 
Interpretations. The posssbility of an eclogite layer in the upper 
mantle was also discussed by Press (43). 

The idea that eclogite may be the source for oceanic basalts is an 
old one (31) but has not been in favor in recent years. The main 
objection is that limited partial melting of eclogite does not generate 
a tholellte. Extensive or complete melting is required and this has 
been thought to be unlikely. 

However, eclogite has an extremely small melting Interval (about 
60“C) compared to peridotite (16, 32). Since diapirs cannot rise out of 
an eclogite layer into a less dense peridotite layer until they are 
already extensively molten, it requires only a small additional 
temperature rise, relative to the liquldus, to complete the melting. 
This can be accomplished in a rising diaplr. It is still unknown why 
melting initiates in the first place. One possibility is insulation by 
the thick continental lithosphere. 

Hotspots and ridges in the Atlantic and Indian oceans were beneath 
continental lithosphere prior to 200 m.y. ago. A large number of 
hotspots are presently under the relatively stationary continent of 
Africa. This suggests that continental insulation, which prevents 
mantle heat from being efficiently removed to the surface, may be the 
cause of the partial melting that ultimately results in ridge and 
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hotspot volcanlsm. The thick* 150 km* continental lithosphere (20) 
prevents dlaplrlc uprise and extensive melting may be possible until the 
continent rifts and rapid spreading and heat removal can take place* It 
Is not so obvious that th.n explanation holds for the ridges and 
hotspots in the Pacific. The thick crust under the oceanic plateaus in 
the central Pacific or a thick oceanic lithosphere* due to a period of 
slow spreading* could also serve to insulate the mantle* raise 
temperatures and cause extensive ' in-sltu melting prior to diapirlc 
ascent. Parts of the Pacific rim continents such as Alaska* Mexico, 
Central America and S.E Asia may also have been located in the the 
central Pacific prior to 200 m.y. ago. Once initiated* the rise of 
upper mantle diapirs should be rapid (33), possibly rapid enough to 
avoid crystal fractionation en route. 

Conclusio ns 

Garnet and possibly cllnopyroxene control the complementary trace 
element patterns of MORB and plume basalts. I suggest that the source 
region for MORB was formed from an eclogite cumulate layer resulting 
from crystallization of the picritlc melt fraction of the original 
differentiation of the Earth. The plume source is enriched in those 
elements that would be concentrated in a melt that was removed from a 
garnet rich region. This fluid, the result of an early partial melt or 
a late stage fluid from crystallizing eclogite in the MORB source, 
infiltrated the shallow mantle plume source at various times. This 
gives plume basalts an inverse garnet trace element signature relative 
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to MORR. The MORB and plume source regions together make up about 
25-30% of the mantle. This requires that roost, if not all, of the 
mantle has been processsed to obtain the observed enrichment. There Is 
no need to invoke a large primitive reservoir in the mantle or deep, 
lower mantle plume sources. On the contrary, the plume source region 
appears to be shallow and may be coincident with the LVZ. The upper 
mantle transition region, 220-670 km, represents about 21% of the mantle 
and the LVZ is about one-third as large. This is about the ratio of the 
sizes of the two upper mantle reservoirs which is required to satisfy 
the trace element data. The MORB source is probably an eclogite 
cumulate while the plume source is probably garnet peridotlte. 
Crystallization of a magma ocean would give this kind of upper mantle 
stratigraphy. Nd and Sr isotopic data indicates that the depletion of 
that part of the MORB source region being sampled today occurred over 
the InLerval 1.5 to 2.5 x 10^ years ago. The fluid Involved in the 
depletiou/entlchraent events, and in the formation of the continental 
crust, appears to be a partial melt or late stage residual fluid from an 
eclogite cumulate layer. The lower mantle is residual peridotite and, 
although it may be convectlng, it no longer communicates with the 


shallower reservoirs. 
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Figure Captions 

Figure 1. Normalized trace element concentrations in the continental 
crust (dots), continental basalts and mid-ocean ridge basalts 
<2, 6, 8, 34, 35). Ml concentrations are normalized to 
terrestrial values of Ganapathy and Anders (37), recalculated 
to mantle equivalents. The lower curves for continental 
tholelltes and MORB are normalized to mantle concentrations 
after removal of the continental crust, as are the curves for 
alkali basalts and kimberlites. The complementary nature of 
MORB and the other magma types is evident. The continental 
crust is not the only "enriched" reservoir. 

Figure 2. Crystal/liquid partition coefficients for garnet, 
cllnopyroxene, orthopyroxene and olivine (6). Heavy line is 
ratio of concentrations in primitive mid-ocean ridge basalts 
(36) and continental tholelltes (2, 6). This indicates that 
MORB results from melting of a garnet and cllnopyroxene rich 
source region and that the continental tholeilte (CFB) source 
region has experienced ecloglte extraction. The points are 
for a model (38) in which the proto-MORB source region is an 
ecloglte orthocunulate from the Initial melt fraction of 
primitive mantle differentiation. The continental and ocean 
island basalt, or plume, source region is depleted residue 
enriched by a late stage melt from the proto-MORB cumulate 


layer 
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Figure 3. Present values (39) for an ecloglte cumulate 

(solid line) which has been depleted at various times by 
removal of a melt fraction (AO). This melt fraction 
enriches the complimentary reservoir (dashed line)* Data 
from references 2, 2A| 25| 29, and 30. 

Figure A. A model for the evolution of the mantle. Primitive mantle 
(1) is partially molten either during accretion or by 
subsequent whole mantle convection which brings the entire 
mantle across the solidus at shallow depths. Lit* elements 
are concentrated in the melt. The deep magma ocean (2) 
fractionates into a thin plagioclase-'rlch surface layer and 
deeper olivine-rich and garnet-rich cumulate layers (3). 
tate stage melts in the eclogite cumulate are removed (A) 
to form the continental crust (c.c.), enrich the perldotite 
layer and deplete MORBS, the source region of oceanic crust 
(o.c.) and lower oceanic lithosphere. Partial melting of 
PLUME (5) generates continental flood basalts (CFB), ocean 
island basalts (OIB) and other enriched magmas, leaving a 
depleted residue (harzburgite) that pcays in the upper 


mantle. 
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Although the mantle of the Earth. appears quite rigid on the time scales 
of seismic waves and solid Earth tides It behaves as a fluid for longer term 
processes such as postglacial rebound, plate tectonics and removal of heat 
from the deep Interior. The mantle Is convectlng and will convect until all 
the original heat and heat sources are removed. Although convection itself is 
not questioned, geophysicists are currently debating whether the mantle Is 
convectlng as a whole or is divided into layers which are convectlng separately. 
A major driving force for convection is thermal buoyancy although the plates 
and subdue tlrn, slabs also affect motions in the mantle. Material deep In the 
mantle warms up either due to its own radioactivity or due to heat brought to 
It from below. As it heats up, it becomes buoyant and rises. Similarly, hot 
material brought to the surface cools off and becomes denser. In a homogeneous 
planet the cold surface material eventually sinks back Into the Interior. 
However, thermal buoyancy Is only part of the story, In a chemically stratified 
or dif ferrntiated planet, convection will be confined to the individual layers 
if thermal expansion and contraction cannot overcome the Intrlnlslc density 
differences between layers. Density differences due to chemistry, phase 
changes and partial melting are also Important in mantle convection. For 
example, by partial melting a mantle dlaplr can become much more buoyant than 
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it could by thermal expansion alone and It can therefore rise to higher levels 
In the mantle. Ou the other hand, If a l/"^ge volume of basaltic melt forms 
below about 50 km depth It will free s ecloglte, which Is denser than 
normal mantle, and It will therefore eventually sink. If the Earth had never 
experienced partial melting or differentiation It would be reasonable to 
suppose that the mantle Is convectlng as a whole. Such a mantle would be 
chemically homogeneous, and would not be very Interesting. 

The two mantle source regions 

There Is now abundant geophysical, petrological and geochemical data 
that points toward a chemically Inhomogeneous mantle. Mid-ocean ridge basalts 
come from a region of the mantle that Is depleted In most of the large-ion 
Incompatible elements such as Rb, Sr, Ba, U and the rare-Earth elements. 

These elements are not retained effectively by the major mantle minerals, 
olivine, garnet and pyroxenes. They therefore become concentrated In melts 
Abyssal tholeiices, however, originate as melts or partial melts and this 
means that their source region Is depleted In these elements. The part of 
the mantle that Is providing these voluminous magmas must therefore have 
originated as a crystal cumulate or have experienced a prior episode of 
partial melting and melt extraction. 

Alkali basalts, on the other hand, whether found on continents or 
on oceanic Islands, are enriched In volatiles and the incompatible elements. 
The enrichment is not uniform but Is complementary to the depletion found In 
ocean floor basalts. That Is to say these basalts are most enriched In the 
elements that are most depleted In MORE (Midocean Ridge Basalts), and vice 
versa. The trace element patterns of continental and ocean island basalts 
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are also complementary to garnet. For example these magmas are particularly 
enrich**! in K, Rb» Ba» Th» U and the light rare-Earth elements and are 
relatively depleted in Y» Zr and the heavy REE, the elements favored by 
garnet. This suggests that garnet has either been left behind or was 
previously removed from the source region of the enrichad magmas. Since 
some of these magmas, particularly the continental tholeiites, represent 
large degrees of partial melting, It is unlikely that garnet Is still a 
residual phase. Olivine and orthopyroxene are the probable main residual 
phases and these are essentially free of the large-ion incompatible elements. 
It appears, therefore, that the mantle source region for the enriched basalt 
types, such as nephellnltes, melelltltes, basnnltes, alkali basalts, 
continental flood basalts and kimberlites, had previously experienced an 
episode of garnet or ecloglte removal. The continental crust has a similar 
enrichment pattern and probably came from the same reservoir. 

The MORB characteristics can be explained if the "depleted" reservoir 
is eclogltlc since these basalts are enriched in those elements which are 
concentrated in garnet and, to a lesser extent, cllnopyroxene. The Important 
role of ecloglte removal and crystal fractionation has been emphasized by 
O’Hara (1968) and O’Hara and Yoder (1967). 

A terrestrial magma ocean 

A similar complementarity among ]Lunar basalts has led' to the concept 
of a crystallizing magma ocean (Wood, 1975; Smith et al, 1970; and 
Anderson, 1973). The thick anorthositlc highland crust is thought to 
represent light plagioclase-rich cumulates and the other basalt types have 
been attributed to remeltlng of denser and deeper olivine-pyroxene rich 


A-V11> 


cumulates, The complementarity of the various reservoirs is controlled by 
plagioclase separation. If the moon was extensively molten early in its early 
history then the Earth and the other terrestrial planets should have been as 
well. On a large body* however, garnet crystallization is Important, Garnet 
settling is the terrestrial analog of plagioclase floatation. A thin plaglo- . 
clase-rlch crust will form in the surface chill layer of a terrestrial magma 
ocean but much of the AI2OJ will be removed from the magma by garnet 
fractionation before temperatures in the. bulk of the ocean are low enough 
for plagioclase formation. 

A thick eclogite ctimulate layer will therefore form at depth in an 
Earth-sized planet. Crystallization of a magma ocean will lead to a stratified 
mantle and the source regions will bear the imprint of garnet or garnet 
extraction. A deep garnet-rich cumulate layer would be an appropriate source 
region for MORB. A possible mantle stratigraphy resulting from crystallization 
of a terrestrial magma ocean is shown in Figure 1. 

How does a magma ocean form and what will be Its composition? The 
energy of accretion of the Earth Is sufficient to vaporize most of the mantle. 
Much of the energy of accretion is radiated away but most of the mantle was 
probably raised above the solidus during accretion. Differentiation, 

Including segregation of the core, would be occurring while the Earth was 
being assembled. Silicate melts remain at the surface and the residual 
crystals settle to the interior. The proto-Earth would therefore be 
concentrating Al20^, CaO, U, Th, Ba, Sr, etc. in the upper mantle while It was 
growing. The concentration of Incompatible elements In the upper mantle 
can be explained If the whole mantle experienced about 15% partial melting 
(Anderson, 1980) . 
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Even if one could assemble a cold and homogeneous Earth-size body 
it would experience whole mantle convection that would bring most or all 
of the mantle through the upper mantle melting zone. Separation of melt and 
solid would ensue, Just as is happening at ridges today. 

The higher abundances of the radioactive elements in the past means 
that the cold conductive surface layer was thinner and that the deep mantle 
adiabatic gradient extended closer to the surface. The upper mantle melting 
zone was therefore broader in the Precambrlan and melting was more extensive. 
The requirement that the larger amount of heat be removed results in a thin 
conductive thermal boundary layer with a high temperature gradient. If a 
significant fraction of this layer is basaltic then it will be buoyant and 
will be unable to subduct. Saibduction only occurs today because part of the 
lithosphere is composed of material denser than basalt and because cooling 
occurs over a ritther large depth extent. Areas having thick crusts, such as 
continents, oceanic plateaus and aselsmlc oceanic ridges resist subductlon. 
Early Precambrlan plate tectonics was likely quite different in style than 
modern plate tectonics. It would be dominated by colllslonal tectonics, 
shallow under thrusting and remeltlng of over thickened crust. As the outer 
layer cools the lower basaltic parts eventually convert to eclogite and a 
massive overturning would result. This presumably erases the early geological 
record and sets the stage for the current style of t;>late tectonics which 
Involves deep subductlon and continuous regeneration of ocean lithospheres. 

The Earth is unique among the terrestrial planets in having plate 
tectonics and^ probably, in h- *in;B a very thin crust. These phenomena are 
related since oceanic crust Is continuously being returned to the mantle. 

In principle, the Earth could support a basaltic crust about 50 km thick. 
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At this depth basalt transforms to ecloglte. A widespread ecloglte layer, 
however, Is unstable and would sink Into the mantle. This transformation 
occurs at much greater depths In Mars and the Moon, because of the lower 
pressures, and In Venus, because of the higher surface and upper mantle 
temperatures. Tlierefore, thick and permanently buoyant crusts are expected 
on these planets. An equivalent quantity of material Is presumably buried 
somewhere in the Earth's mantle, Where Is It, how big la It and how long 
has It been there? 

The last question Is the easiest to answer with presently available 
data. Isotopes tell us that the various source regions of the mantle have 

q 9 . 

been Isolated for at least 10 and possibly more than 3 x 10 years (Patterson 
and Tatsumoto, 1964} Tatsumoto, 1978; De Paolo and Wasserburg, 1979; Jacobsen 

g 

and Wasserburg, 1979). This Is much longer than the n#10 year average age of 
the oceanic crust. The Isolation of these reservoirs may have been a result 
of the initial differentiation of the mantle. The size of the depleted 
reservoir is about ‘?. 0 % of the mantle, which Is equivalent to a 450 km thick 
layer in the upper mantle (Anderson, 1980). 

In aiost subduction areas of the world earthquakes occur as deep as 
220 ton. In a few places, particularly where the oceanic lithosphere is old, 
earthquakes occur as deep as 670 km, the depth of a major seismic discontinuity 
in the mantle. Therefore, oceanic lithosphere appears to be disposed of 
between these depths. This is the transition region of the mantle, 

I have shown elsewhere (Anderson, 1979a, b) that ecloglte is denser 
than garnet peridotite to a depth of 670 km. Below this depth normal low 
aluminum mantle minerals transform to ilmenite and perovskite polymorphs. 

These are denser than garnet, the main mineral of ecloglte at these depths, 
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an<! Chey will form a floor to subducting crust or a garnet cumulate layer 
forming In a crystallizing magma ocean. It appears that the 670 km discontinuity 
Is not only a phase boundary, separating garnet from llmenlte/perovskltc, 
but Is also a chemical boundary, separating eclogltlc from ultrabaslc rocks. 

The lower mantle Is presumably the crystalline residue remaining 
after extraction of a basaltic, probably plcritlc, melt. This melt, the 
result of early mantle differentiation, fractionates Into eclogltlc and 
perldotltic cumulate layers, both of which currently reside In the upper 
mantle. Because of the phase change boundary the differentiation of the 
Earth Is Irreversible. Mixing of the upper mantle reservoirs Is Inefficient 
because of the large density contrast between eclogite and garnet perldotlte. 
Normally, the two upper mantle reservoirs convect separately because of this 
density stratification. However, if the eclogite layer partially melts It 
becomes less dense than the overlying mantle. Partially molten dlaplrs can 
then rise out of the eclogite layer and melt completely during adiabatic 
ascent. This is possible because of the small melting interval of eclogite. 

The dlaplrs then provide plcritlc melts to near surface magma chambers where 
they fractionate to form midocean ridge tholelltes and, possibly, a pyroxenitic 
lower lithosphere. The latter evolves to garnet pyroxenlte as the litho- 
sphere thickens and cools. This Is denser than normal mantle perldotlte 
and the oceanic lithosphere eventually becomes unstable and sinks back to 
its source region. 

The alternative point of view Is that ocean floor tholelltes are 
primary magtaas, unaffected by crystal fractionation, which are formed by 
partial melting of garnet perldotlte, or pyrolltc, leaving behind a thick 
depleted harzburglte layer. Since basalt and harzburglte are both buoyant 
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In norma] mantle, a large amount of cooling la required in order to overcome 
this chemical stabilization. A thick layer of depleted perldotlte in the 
upper mantle is one consequence of this model. 

Hotspots and plumes 

The other major reservoir in the mantle provides magmas for the oceanic 
islands and the alkalic provinces on continents. Continental flood basalts 
apparently come from this same reservoir. Current sites of this kind of 
igneous activity are called plumes or hotspots (Wilson, 1963; Morgan, 1971)'. 

In many respects these basalts are complementary to midocean thol elites and 
appear to come from an enriched source region. What is the size and location 
of this reservoir? 

Although It is usually assumed that plumes originate deep in the 
mantle, there Is a variety of evidence that indicates they come from a 
shallower reservoir than MORB. The limited areal extent of hotspots suggests 
that they are caused by narrow upwelllngs in the mantle. This is typical 
behavior of fluids heated from below and Is consistent with an upper mantle 
location. Uniform fluids which are heated from within exhibit broad, often 
roll-like, upwelllngs. This pattern is more appropriate for oceanic ridges. 
Internal heating of a deep upper mantle layer would generate internal rolls 
and plume-like convection In the overlying layer. Kimberlites and their 
entrained xenollths are enriched in volatiles and the incompatible elements 
and these come from the upper 250 km or so of the mantle. Xenollths brought 
to the surface in alkali basalt flows also apparently originate In the shallow 
mantle. The nature of the enrichment in continental and ocean Island basalts 
is such that their source region could be the upper mantle low-velocity zone. 
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Partial malting In tha uppar mantle can explain the low aalamlc valocltlea 
and high electrical conductivity in tha low-velocity zona. Partial malting la 
conalatent with estimates of uppar mantle temperatures If this region Is 
volatile rich. 

The appearance of enriched magmas at developing continental rifts, 

Island arcs, and back-arc basins also suggests a shallow source origin. Trace 
element patterns Indicate that the enriched reservoir Is much smaller than 
the depleted reservoir, perhaps only one^thlrd as large. Flume basalts 
represent varying degrees of partial melting. This Is consistent with a small 
and variable depth of origin. Deep dlaplrs, brought adlabatlcally to the 
surface, would exhibit a large and more uniform degree of partial melting. 

This property is more typical of MORB than the hotspot magmas. 

Theoretical calculations Indicate that the Earth Is getting rid of 
Its heat at about the same rate as It Is generating It by radioactive decay, 
Convection Is an efficient transmitter of heat and the main bottleneck Is 
the cold layer at the surface that must transmit heat by conduction, There 
are also conductive boundary layers at each chemical Interface In the Earth, 
Because of the temperature dependence of viscosity, these boundary layers 
adjust their thickness to accommodate changes In heat flow and temperature. 

The current heat flow through the surface of the Earth Is probably within 
30% of the current rate of heat productlo.i. Below about 250 km the temperature 
gradient ia probably near adiabatic. In addition to the near surface thennal 
boundary layer, sometimes referred to as the thermal lithosphere, there Is 
also a region of high thermal gradient across the lithosphere-asthenosphere 
boundary which Is thought to be at about 150 km under shields (Anderson, 1979c). 
Pressure Increases the melting point of mantle silicates so rapidly that 
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melcing is unlikely to occuv in a convacting mantle at deptha below about 
250 km. Solid material brought adiabatically from depth, however, will 
Intersect the melting curve in the upper mantle. Thus, upper mantle melting 
Is a natural consequence of mantle convection and the relative slopes of the 
adiabat and the solidus. The presence of volatllas also decreases t*he solidus 
In the upper 150 km of the mantle. Water and CO^ are locked In mineral 
phases at greater depths and the melting point Is not much affected. 

There Is another way In which excess temperatures can be generated 
In the upper mantle. The thin lithosphere under oceans allows a balance 
between the rate at which heat Is convected Into the upper mantle and the 
rate at which It can be conducted to the surface. The lithosphere Is much 
thicker under shields and the geothermal gradient Is relatively low. This 
means that the mantle heat flow cannot be rapidly removed through a thick 
continent. A stationery continent will, therefore, heat up at Its base and 
hotspots can be generated In the mantle beneath large continental masses 
that are stationery or slowly moving. These self -generated hotspots weaken 
and uplift the continent and are probably responsible for continental breakup 
and dispersal. 

It is instructive to plot the locations of present day hotspots In 
the Atlantic and Indian oceans relative to the locations of continents prior 
to the breakup of Gondwanda. Figure 2 shows such a display. The currently 
active hotspots in this hemisphere are In a part of the mantle that was 
insulated by a super continent from before 350 to 200 my ago, l.e., for more 
than 150 my prior to breakup and dispersal. The hotspot grouping In the 
tenter of the map is currently under Africa and those under North America 
are currently along the mid-Atlantic ridge. Two hundred million years ago 
the Eiffel plume was under northern Europe, Iceland and Jan Mayan were under 
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Siberia and St. Helena, Trlatan and Bouvet were under South America, In 
many cases continental uplift, alkallc basalt provinces, continental flood 
basalt and kimberlite activity can be associated directly with a hotspot 
track. I suggest, therefore, that continental insulation is responsible 
for melting anomalies in the mantle. These melting anomalies, in turn, lead 
to uplift, intraplate magroatism and the breakup of large continental masses. 
Melting in the shallow perldotite layer gives enriched continental and ocean 
Island magmas and melting in the deeper ecloglte layer provides the depleted 
magmas wliich emerge at well developed ridges. 

Similar mantle insulation is provided by anomalously thick oceanic 
crust such as the high standing oceanic plateaus in the western Pacific. 

These plateaus may have been responsible for turning on the central Pacific 
hotspots, 

On keeping an opan mind 

It is Appropriate at the start of the joint lUGG, lUGS and ICSU program 
on the Lithosphere to review our assumptions and working hypotheses, Geological, 
geophysical and geochemical data are subject to multiple interpretations, and 
models derived by workers in one discipline are often contradictory to those 
derived by those working in a different field. It is often the case, however, 
that the data are not contradictory, only the models. 

For example, the following assumptions are commonly made: 

1) The majority of the mantle Is relatively homogeneous in composition 
and primitive (pyrollte); 

2) Basalts of all types represent varying degrees of partial melting 
of garnet perldotite or pyrollte. 
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3) Oceanic tholelltes are primary liquids formed by partial melting 
and the upper oceanic mantle is the depleted residue. 

4) The oceanic lithosphere subducts as a result of cooling and thermal 
contraction. 

5) The seismic discontinuities in the mantle represent phase changes In 
an ollvlne-rlch mantle rather than chemical boundaries. 

6) The large sizes of the rigid plates and the separation of ridges 
implies large convection cells having dimensions that correspond to the depth 
of the mantle. 

7) The lower mantle Is primitive and provides the material for hotspots, 
oceanic Islands and continental flood basalts. 

8) The enormous quantities of oceanic crust that have been generated 
by sea-floor spreading have sunk to the core-mantle boundary. 

Many of these assumptions have recently been called Into question. 

I suggest that alternatives to each of the above assumptions are viable and they 
must be considered as we strive toward a coherent theory of the evolution of 
the mantle and plate tectonics. The evidence for early melting and differentatlon 
of small solar system objects make It difficult to understand how the Earth 
could have remained primitive, homogeneous and undifferentiated. The processes 
of crystal fractionation and gravitational stratification are important In 
the Moon and In terrestrial magma chambers and were also probably Involved 
in the early evolution of the mantle. Petrological models for the mantle are 
strongly influenced by samples from the upper mantle which, In many cases, 
are clearly residual after melt removal, or crystal cumulates. Thus, there 
is a strong bias toward an olivine rich mantle. Even primitive appearing 
specimens are samples i^om the upper mantle and do not necessarily represent 
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the average composition of the whole mantle« When the large quantities of 
subducted oceanic crust are taken Into account, the average mantle becomes 
more CaO, ^^2^3* closer to cosmochemical estimates 

of mantle composition. Workers who propose tie existence of a large primitive, 
undifferentiated reservoir In the mantle overlook the fact that elements 
such as Rb, Sr, Sm and Nd are so strongly concentrated into melts that by 
the time one melts <^20% of primitive mantle the melt has essentially "primitive" 
ratios of Rb/Sr, Sm/Nd and their Isotopes. The ratios stay nearly primitive 
even after large .mounts of olivine and orthopyroxene fractionation. Garnet 
and cllnopyroxenc fractionate these elements but orthocumulatcs composed of 
these minerals and equilibrium interstitial fluids give nearly primitive 
ratios. The "primitive reservoi *" theories assume that the continental crust 
is the only enriched reservoir that is complementary to the depleted MORB 
reservoir. In che "differentiated mantle" theory outlined in this essay 
the source region of continental and ocean island basalts is the enriched 
complanenC of the depleted reservoir. Isotopic studies of more hotspot 
locations and mantle xenoliths should resolve these two models. 

The study of the oceanic and continental lithospheres is of central 
importance in unraveling the evolution of the mantle. In the fractional 
crystallization model the lower oceanic lithosphere is pyroxenltlc rather than 
harzburgite or pyrolite. The spinel to garnet pyroxenlte phase change would 
be the luuln cause of instability of old lithosphere. "Eclogite" xenoliths 
found in Hawaii may be samples of this lower lithosphere. Study of the 
lithosphere under shields may shed light on early mantle processes and the 
development of the contincntt.x crust. Shields exhibit long term stability in 
spite of the fact that they are cold. This may result from chemical buoyancy. 
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It is possible that chernii ?? effects, i.e., partial melting and phase changes, 
may be as important in driving the terrestrial engine as thermal convection. 

The evolution of a primitive mantle into two upper mantle reservoirs 
and a depleted lower mantle is shown in Figure 3, The formation of hotspot 
magmas by partial melting of garnet peridotlte in the upper mantle probably 
follows conventional views of petrogenesis. Trace elements indicate that this 
region has experienced eclogite extraction and enrichment in the incompatible 
elements. Diapirs from the eclogite cumulate layer are presumed to evolve 
to mid-ocean ridge tholeilres. The long term isolation of the two major source 
regions strongly supports the concepts of chemical stratification in the 
mantle and chemical driving forces, l.e., partial melting and phase changes 
in the lithosphere, in addition to thermal buoyancy. 
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FIGURE CAPTIONS 

Figure 1 Schematic illustration of the sequence of layers in a crystallizing 
terrestrial magma ocean. The magma ocean is the result of partial 
melting of a primitive garnet perldotite mantle. This melt, probably 
plcrltlc, collects In the upper mantle. A garnet-rich cumulate 
layer fok,;s at the base of the molten layer leaving behind an 
Al20^ and S102~poor liquid. Shallower cumulate layers consist of 
olivine, orthopyroxene and interstitial melt. Olivine and ortho- 
pyroxene react with the fluid to form cllnopyroxene and garnet if 
they settle below about 100 km. These cumulate layers evolve to 
the two major reservoirs which are being sampled today, The lower 
mantle is depleted perldotite, the residue after about 15% melt 
removal. The upper boundary of the lower mantle is controlled by 
phase changes in perldotite. Because of the broad stability field 
of gai net, ecloglte cannot sink below 670 km. The 670 km discontinuity 
Is, therefore, a chemical and phase boundary. 

Figure 2 The location of the continents at about 350 my ago. Also shown are 
Che hotspots in this hemisphere, most of which are currently in the 
Atlantic and Indian oceans and under the continent of Africa. The 
continents move slowly north for the next 150 ray and then break up 
and disperse to their present locations. The majority of present 
day hotspots were beneath Gonwandaland for a long period of time 
prior to 200 my ago. They may have been formed by continental 
insulation. 
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Figure 3 Evolution of primitive mantle into a chemically atratifled mantle 
with separate isolated upper mantle reservoirs. Melt is removed 
from primitive mantle, probably during accretion, and concentrated 
at the surface. As the Earth cools, the melt fractionates into 
a deep eclogite. cumulate layer and a shallow perldotlte layer. 

The eclogite cumulate becomes depleted in large-*lon lithopile elements 
(LIL) which enrich the overlying perldotlte layer. Dlapirs from 
the eclogite source region evolve to mid-ocean tholelltes by near 
surface crystal fractionation. The enriched perldotlte source 
region provides, by varying degrees of partial melting, kimberlites, 
melilltes, nephelinites, alkali basalts and continental tholelltes. 
Depleted perldotites are the residue of hotspot magmatlsm. These 
are less dense than fertile perldotites and therefore occur in 
the uppermost mantle. The oceanic lithosphere, at least the Al^O^- 
rlch portions, return to the eclogite layer. 
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